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Introduction: 


The  role  of  Wnt  signaling  in  epithelial  ovarian  cancer  (EOC)  development 
remains  largely  elusive.  We  have  evidence  to  suggest  that  canonical  Wnt  signaling  is 
activated  through  downregulation  of  Wnt5a  in  EOC.  Wnt5a  is  a  non-canonical  Wntligand 
that  acts  as  an  antagonist  of  canonical  Wnt  signaling.  Furthermore,  inhibition  of 
canonical  Wnt  signaling  or  restoration  of  Wnt5a  expression  inhibits  the  growth  of  EOC 
cells.  The  overall  hypothesis  of  this  proposal  is  that  canonical  Wnt  signaling  activated  by 
loss  of  Wnt5a  contributes  to  EOC  development  by  abrogation  of  OIS.  The  objectives  of 
the  proposed  studies  are  to  determine  the  role  of  canonical  Wnt  signaling  in  regulating 
OIS  of  ovarian  epithelial  cells  during  progression  of  benign  ovarian  tumors  into  invasive 
EOCs,  and  to  investigate  the  effects  of  inhibition  of  the  canonical  Wnt  signaling  on 
malignant  behavior  of  EOC  cells.  The  specific  aims  are:  1):  Elucidate  the  mechanisms  by 
which  inhibition  of  canonical  Wnt  signaling  inhibits  the  growth  of  EOC  cells;  2): 
Determine  whether  loss  of  Wnt5a  contributes  to  EOC  development  by  abrogation  of  OIS; 
3):  Investigate  the  effects  of  inhibition  of  canonical  Wnt  signaling  on  malignant  behavior 
of  EOC  cells  in  immunodeficient  mice. 

Body: 

1.  Research  accomplishments  associated  with  each  task  outlined  in  the  approved 
Statement  of  Work. 

Task  1.  Determine  whether  inhibition  of  canonical  Wnt  signaling  inhibits  the 
growth  of  EOC  cells  through  inducing  cell  senescence.  (Months  1-12) 

Specifically,  we  will  detennine  whether  inhibition  of  canonical  Wnt  signaling  induces  the 
expression  of  markers  of  senescence  in  human  EOC  cells. 
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Figure  1.  FJ9  suppresses  the  growth  of  human  EOC  cells  via  Wnt/|3-catenin 
independent  mechanisms.  (A)  Equal  number  of  indicated  EOC  cell  lines  were 
treated  with  lOOuM  FJ9  or  vehicle  control  (dissolved  in  0.01%  DMSO).  After  72 
hours  of  treatment,  the  number  of  cells  were  counted  and  calculated  as  relative 
number  of  cells.  (B)  same  as  (A)  but  examined  for  cell  morphology  changes  using 
phase  contrast  microscopy.  (C)  Same  as  (A)  but  examined  for  expression  of  soluble  b- 
catenin  and  b-actin  in  indicated  human  EOC  cell  lines  treated  with  FJ9  or  vehicle 
control  by  immunoblotting  using  indicated  antibodies. 


Progress  report: FJ9  has  previously  been  demonstrated  as  an  inhibitor  of  canonical  Wnt 
signaling  [1].  We  showed  that  FJ9  inhibits  the  growth  of  human  EOC  cells  in  a  panel  of 
human  EOC  cell  lines  (Figure  1  A).  Consistent  with  the  idea  that  this  is  due  to  senescence 
induction,  cells  treated  with  FJ9  demonstrated  features  of  senescence  such  as  a  large  flat 
cell  morphology  (Figure  IB).  However,  examination  of  markers  of  canonical  Wnt 
signaling  in  FJ9  treated  or  control  cells  showed  that  there  is  no  evidence  to  suggest  the 
observed  effects  are  due  to  inhibition  of  canonical  Wnt/ (beaten  in  pathway.  For  example, 
the  levels  of  soluble  [beaten in,  a  marker  of  active  canonical  Wnt/ (beaten in  pathway  were 
not  decreased  by  FJ9  treatment  (Figure  1C).  This  result  suggests  that  FJ9  inhibits  the 
growth  of  human  EOC  cells  via  canonical  Wnt  signaling  independent  mechanisms. 
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Figure  2.  Wnt5a  restoration  inhibits  the  growth  of  human  EOC  cells  by 
antagonizing  canonical  Wnt/p-catenin  signaling.  (A)  OVCAR5  cells  were 
transduced  with  a  control  or  Wnt5a  encoding  puromycin  resistant  retrovirus.  The 
infected  cells  were  drug-selected  with  3  ug/ml  puromycin.  Expression  of  Wnt5a 
mRNA  in  drug-selected  cells  was  determined  by  qRT-PCR.  (B)  Same  as  (A),  but 
examined  for  expression  of  Wnt5a  and  P-actin  in  control  or  Wnt5a-infected  OVCAR5 
and  primary  HOSE  cells  by  immunoblotting.  Relative  levels  of  Wnt5a  expression  was 
indicated  based  on  densitometric  analysis  using  NIH  Image  J  software.  (C)  Same  as 
(A),  but  equal  number  (5000)  of  drug-selected  control  (open  triangles  and  dotted  line) 
or  Wnt5a-infected  cells  (open  circles  and  solide  line)  were  cultured  on  plastic  plates 
for  4  days,  and  the  number  of  cells  was  counted  (Control  ±  SD  or  Wnt5a  ±  SD  (n=3); 
student  t-test  was  used  for  calculating  p  value)  *p  =  0.0038  and  **p  =  0.001.  Mean  of 
three  independent  experiments  with  SD  and  linear  regression.  (D)  Same  as  (C),  but 
grown  under  anchorage-independent  condition  in  soft-agar.  The  number  of  colonies 
was  counted  2  weeks  after  initial  inoculation.  Mean  of  three  independent  experiments 
with  SD.  (E)  Same  as  (A),  but  examined  for  the  levels  of  soluble  p-catenin  and  p-actin 
expression  by  immunoblotting  (NT  =  non-treated).  (F)  Same  as  (A),  but  examined  for 
expression  of  indicated  p-catenin  target  genes  by  qRT-PCR.  Expression  of  /?-2- 
microglobulin  was  used  to  normalize  the  expression  of  indicated  genes.  *  p  =  0.0095, 
**  p  =  0.0012  and  ***  p  =  0.0286  compared  with  controls. 


We  sought  to  determine  the  effects  of  Wnt5a  reconstitution  in  human  EOC  cells.  Wnt5a 
expression  was  reconstituted  in  the  OVCAR5  EOC  cell  line  via  retroviral  transduction. 
Ectopically  expressed  Wnt5a  was  confirmed  by  both  qRT-PCR  and  immunoblotting  in 
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0VCAR5  cells  stably  expressing  Wnt5a  or  a  vector  control  (Figure  2A-B).  Of  note,  the 
levels  of  ectopically  expressed  Wnt5a  in  OVCAR5  cells  are  comparable  to  the  levels 
observed  in  primary  HOSE  cells  (Figure  2B).  Interestingly,  Wnt5a  reconstitution  in 

OVCAR5  human  EOC  cells 
significantly  inhibited  both 
anchorage-dependent  and 

anchorage-independent  growth 
in  soft-agar  compared  with 
vector  controls  (Figure  2C-D). 
In  addition,  similar  growth 
inhibition  by  Wnt5a 
reconstitution  was  also  observed 
in  the  PEOl  human  EOC  cell 
line  (Figure  3A-C),  suggesting 
that  this  effect  is  not  cell  line 
specific.  Based  on  these  results, 
we  conclude  that  Wnt5a 
reconstitution  inhibits  the 
growth  of  human  EOC  cells  in 
vitro. 

Canonical  Wnt  signaling 
promotes  cell  proliferation  and 
Wnt5a  has  been  demonstrated  to 
antagonize  the  canonical 
Wnt/  (3-catenin  signaling  in 
certain  cell  contexts  [2-5].  We 
hypothesized  that  Wnt5a  would 
suppress  the  growth  of  human 
EOC  cells  by  antagonizing 
canonical  Wnt/fi-catcnin 

signaling.  To  test  our 
hypothesis,  we  examined  the 
effect  of  Wnt5a  reconstitution 
on  expression  of  markers  of 
active  Wnt/fl-catenin  signaling 
in  human  EOC  cells,  namely  the 
levels  of  “active”  soluble  |3- 
catenin[6-8]  and  expression  of 
(3-catenin  target  genes  such  as 
CCND1,  c-MYC  and  FOSL1  [9, 
10].  Indeed,  we  observed  a 
decrease  in  soluble  (3-catenin  in 
Wnt5a  reconstituted  OVCAR5 
cells  compared  with  vector 
controls  (Figure  2E). 
Consistently,  we  also  observed  a 
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Figure  3.  Wnt5a  restoration  inhibits  the  growth 
of  PEOl  EOC  cells  by  antagonizing  canonical 
Wnt  signaling.  (A)  PEOl  cells  were  transduced 
with  a  control-  or  Wnt5a-encoding  puromycin 
resistant  retrovirus.  The  infected  cells  were  drug- 
selected  with  1  ug/ml  puromycin  for  3  days. 
Expression  of  Wnt5a  and  GAPDH  in  drug-selected 
cells  was  determined  by  immunoblotting.  (B)  Same 
as  (A)  but  an  equal  number  of  drug-selected  control 
and  Wnt5a  infected  PEOl  cells  were  cultured  on 
plastic  plate  for  5  days  and  the  number  of  cells 
were  counted  using  trypan  blue  exclusion  assay. 
Mean  of  three  independent  experiments  with  SD. 
(C)  Same  as  (B)  but  grown  under  anchorage- 
independent  conditions  in  soft-agar.  The  number  of 
colonies  were  counted  2  weeks  after  initial 
inoculation.  Mean  of  3  independent  experiments 
with  SD.  (D)  Same  as  (A),  but  examined  for 
soluble  [I -eaten in  expression  by  immunoblotting. 
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significant  decrease  in  the  levels  of  [)-catenin  target  genes  in  these  cells,  namely  CCND1 
(p  =  0.0095),  FOSL1  (p  =  0.0012)  and  c-MYC  (p  =  0.0286)  (Figure  2F).  Similar  effects 
of  Wnt5a  reconstitution  on  expression  of  markers  of  active  Wnt/|3-catenin  signaling  (such 
as  decreased  levels  of  soluble  (3-catenin)  were  also  observed  in  PEOl  human  EOC  cells 
(Figure  3D),  suggesting  that  this  is  not  cell  line  specific.  Based  on  these  results,  we 
conclude  that  Wnt5a  suppresses  the  growth  of  human  EOC  cells  by  antagonizing 
canonical  Wnt/|3-catenin  signaling  in  human  EOC  cells. 

Next,  we  sought  to  determine  the  cellular  mechanism  whereby  Wnt5a  inhibits  the 
growth  of  human  EOC  cells.  We  sought  to  determine  whether  Wnt5a  restoration  might 
induce  senescence  in  human  EOC  cells.  To  do  so,  we  determined  whether  Wnt5a 
restoration  induces  SA-|3-gal  activity,  a  universal  marker  of  cellular  senescence  [11]. 
Indeed,  SA-|3-gal  activity  was  notably  induced  by  Wnt5a  reconstitution  in  both  OVCAR5 
and  PEOl  human  EOC  cells  compared  with  controls  (Figure  4). 
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Figure  4.  Wnt5a  restoration  triggers  cellular  senescence  in  human  EOC  cells.  (A) 

Control  and  Wnt5a  expressing  OVCAR5  EOC  cells  were  stained  for  senescence- 
associated  (3-galactosidase  (SA-|3-gal)  activity.  (B)  Quantitation  of  (A).  Mean  of  three 
independent  experiments  with  SD.  (C-D)  Same  as  (A-B)  but  for  PEOl  human  EOC 
cells. 


Task  2.  Determine  the  effects  of  inhibition  of  Wnt/p-catenin  signaling  on  key 
senescence-regulating  pathways  in  human  EOC  cells.  (Month  13-24) 

Specifically,  we  will  determine  the  effects  of  inhibition  of  Wnt/|)- eaten  in  signaling  by 
FJ9  treatment  or  enforced  Wnt5a  expression  on  status  of  p53,  pRB  and  HIRA/PML 
pathways. 
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Progress  report:  Since  FJ9  failed  to  inhibit  canonical  Wnt  signaling  (Figure  1C),  we 
focused  on  our  study  on  using  Wnt5a  restoration  as  a  way  to  inhibit  canonical  Wnt 
signaling  (Figure  2  and  3).  We  have  previously  shown  that  suppression  of  canonical  Wnt 
signaling  promotes  cellular  senescence  in  primary  human  fibroblasts  by  activating  the 
senescence-promoting  histone  repressor  A  (FIIRA)/  promyelocyticleukemia  (PML) 
pathway  [8].  PML  bodies  are  20-30  dot-like  structures  in  the  nucleus  of  virtually  all 
human  cells.  PML  bodies  are  sites  of  poorly  defined  tumor  suppressor  activity,  and  are 
disrupted  in  acute  promyelocytic  leukemia  [12].  PML  has  been  implicated  in  regulating 
cellular  senescence.  For  example,  the  foci  number  and  size  of  PML  bodies  increase 
during  senescence  [12,  13]  and  inactivation  of  PML  suppresses  senescence  [14]. 
Activation  of  the  FIIRA/PML  pathway  is  reflected  by  the  recruitment  of  FIIRA  into  PML 
bodies  [15]. 


C  D 

Wnt5a  -  +  p.pRB  Wnt5a  -  + 

anti-pRB*-*  anti-pS780-RB^— 

anti-GAPDH  —  — *  PRB  anti-GAPDH 

Figure  5.  Wnt5a  restoration  activates  HIRA/PML  pathway  and  decreases  the 
levels  of  phosphorylated  pRB  in  human  EOC  cells.  (A)  Control  and  Wnt5a 
expressing  OVCAR5  EOC  cells  were  stained  with  antibodies  to  HIRA  and  PML. 
Arrows  point  to  examples  of  co-localized  HIRA  and  PML  bodies.  Bar  =  10  p,m.  (B) 
Quantitation  of  (A).  200  cells  from  control  and  Wnt5a  expressing  cells  were  examined 
for  HIRA  and  PML  co-localization.  Mean  of  three  independent  experiments  with 
SD.(C)  Same  as  (A)  ,  but  examined  for  pRB  and  GAPDH  expression.  (D)  Same  as 
(C),  but  examined  for  pRBpS780  and  GAPDH  expression. 
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As  Wnt5a  antagonizes  canonical  Wnt  signaling  in  human  EOC  cells  (Figure  2E- 
F),  we  sought  to  determine  whether  Wnt5a  restoration  might  activate  the  senescence- 
promoting  HIRA/PML  pathway  in  human  EOC  cells.  Towards  this  goal,  we  examined 
the  localization  of  HIRA  in  OVCAR5  EOC  cells  reconstituted  with  Wnt5a  or  vector 


A  peo-i  B 


Figure  6.  Wnt5a  restoration  triggers  senescence  in  human  EOC  cells.  (A)  PEOl 
EOC  cells  were  infected  with  control  or  Wnt5a  encoding  puromycin-resistant 
retrovirus.  The  infected  cells  were  drug-selected  using  1  mg/ml  of  puromycin.  Drug- 
selected  cells  were  stained  with  antibodies  to  HIRA  and  PML.  Arrows  point  to 
examples  of  co-localized  HIRA  and  PML  bodies.  Bar  =10  mm.  (B)  Quantitation  of 
(A).  Mean  of  three  independent  experiments  with  SD 


control.  Notably,  there  was  a  significant  increase  in  the  percentage  of  cells  with  HIRA 
localized  to  PML  bodies  in  Wnt5a  restored  human  EOC  cells  compared  with  controls 
(Figure  5A-B,  p  =  0.004).  In  addition,  we  also  observed  an  increase  in  the  number  and 
size  of  PML  bodies  in  the  Wnt5a  restored  OVCAR5  EOC  cells  (Figure  5A),  which  are 
also  established  markers  of  cellular  senescence  [14,  16]. 

The  p53  and  pRB  tumor  suppressor  pathways  play  a  key  role  in  regulating 
senescence  [11].  Thus,  we  sought  to  determine  the  effects  of  Wnt5a  on  the  p53  and  pRB 
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pathways.  Interestingly,  pl6INK4a,  the  upstream  repressor  of  pRB,  is  deleted  in  OVCAR5 
human  EOC  cell  line  [17].  In  addition,  the  levels  of  total  phosphorylated  pRB  were  not 
decreased  by  Wnt5a,  while  the  levels  of  cyclin  Dl/CKD4-mediated  Serine  780 
phosphorylation  on  pRB  (pRBpS780)  were  decreased  by  Wnt5a  [18]  (Figure  4C-D). 
Further,  p53  is  null  in  OVCAR5  cells  [19].  We  conclude  that  cellular  senescence  induced 
by  Wnt5a  restoration  is  independent  of  both  p53  and  pl6INK4a. 

Similarly,  we  observed  activation  of  the  HIRA/PML  pathway  by  Wnt5a 
restoration  in  PEOl  human  EOC  cells  (Figure  6A-B),  suggesting  that  the  observed 
effects  are  not  cell  line  specific.  Together,  we  conclude  that  Wnt5a  reconstitution 
activates  the  HIRA/PML  senescence  pathway. 

Based  on  these  results,  we  concluded  that  Wnt5a  restoration  induced  senescence 
of  human  EOC  cells  by  activating  the  HIRA/PML  senescence  pathway. 


Key  Research  Accomplishments: 

•  FJ9  suppresses  the  growth  of  human  EOC  cells  but  likely  through  canonical 
Wntsignaling  independent  mechanisms. 

•  Wnt5a  suppresses  the  growth  of  human  EOC  cells. 

•  Wnt5a  inhibits  canonical  Wnt  signaling. 

•  Wnt5a  induces  cellular  senescence. 

•  Cellular  senescence  induced  by  Wnt5a  restoration  is  independent  of  both  p53  and 
pl6. 

•  Cellular  senescence  induced  by  Wnt5a  restoration  correlates  with  activation  of  the 
HIRA/PML  senescence-promoting  pathway. 

•  Cellular  senescence  induced  by  Wnt5a  restoration  correlates  with  a  decreased 
level  of  serine  780  phosphorylated  pRB. 
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Conclusions: 

Wnt5a  promotes  senescence  of  human  EOC  cells  by  suppressing  the  proliferation 
promoting  canonical  Wnt/|3-catenin  pathway.  We  suggest  that  strategies  to  drive 
senescence  in  EOC  cells  by  reconstituting  Wnt5a  signaling  may  offer  an  effective  new 
strategy  for  EOC  therapy 
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The  definition  of  an  oncogene  is  a 
gene  that  actively  promotes  tumori- 
genesis.  For  example,  activation  of  RAS 
oncogene  promotes  cell  transforma¬ 
tion  and  cancer.  Paradoxically,  in  pri¬ 
mary  mammalian  cells,  oncogenic  RAS 
typically  triggers  cellular  senescence,  a 
state  of  irreversible  cell  growth  arrest. 
Oncogene-induced  senescence  is  an 
important  tumor  suppression  mecha¬ 
nism  in  vivo.  Here,  we  discuss  our  recent 
evidence  that  RAS -induced  suppression 
of  DNA  repair  response  via  dissociation 
of  BRCA1  from  chromatin  promotes 
senescence  while  predisposing  cells  to 
senescence  bypass  and  transformation  by 
allowing  for  secondary  hits.  The  molecu¬ 
lar  mechanism  we  uncovered  helps  rec¬ 
oncile  the  tumor-promoting  nature  of 
oncogenic  RAS  with  the  tumor-suppress¬ 
ing  role  of  oncogene-induced  senescence. 

Oncogenic  mutations  in  the  RAS  gene 
are  present  in  -30%  of  human  cancers.1 
Oncogenic  RAS  proteins  promote  cell 
transformation  through  the  engagement 
of  downstream  pathways  such  as  phos- 
phatidylinositol-3  -kinase  (PI3K)/AKT 
and  the  RAF  family  of  serine/threonine 
kinases  such  as  BRAF.2'4  Paradoxically, 
activation  of  the  RAS  oncogene  in  pri¬ 
mary  mammalian  cells  typically  induces 
a  status  of  irreversible  cell  growth  arrest, 
known  as  cellular  senescence.5  By  driv¬ 
ing  irreversible  growth  arrest  of  cancer 
progenitor  cells  harboring  the  initial 
oncogenic  hit,  oncogene-induced  senes¬ 
cence  is  an  important  tumor  suppression 
mechanism  in  vivo.5  Further  under¬ 
scoring  the  importance  of  senescence 


in  tumor  suppression,  reactivation  of 
tumor  suppressors  such  as  p53  trig¬ 
gers  cellular  senescence  and  associated 
tumor  regression  due  to  activation  of  the 
innate  immune  response  in  mouse  mod¬ 
els.6  Recently,  the  senescence-associated 
secretory  phenotype  has  been  proposed 
to  be  a  cell  non- autonomous  mechanism 
by  which  senescent  cells  promote  trans¬ 
formation  of  neighboring  premalignant 
cells.5,7  Here,  we  discuss  our  recent  dis¬ 
covery  of  a  novel  cell-intrinsic  mechanism 
by  which  oncogenic  RAS  drives  cellular 
senescence  while  predisposing  cells  to  sec¬ 
ondary  hits,  which  ultimately  promotes 
senescence  bypass  in  RAS-expressing  pri¬ 
mary  human  cells.8 

Downstream  effectors  of  RAS  such  as 
BRAF  and  AKT  all  possess  the  ability  to 
induce  senescence  on  their  own  in  pri¬ 
mary  mammalian  cells.9,10  Likewise,  inac¬ 
tivation  of  PTEN,  a  negative  regulator  of 
AKT,  also  triggers  senescence  in  primary 
mammalian  cells.11  Notably,  there  are  dif¬ 
ferences  in  the  senescence  induced  by  dif¬ 
ferent  oncogenes.  For  example,  it  has  been 
demonstrated  that  the  senescence  induced 
by  RAS  or  BRAF  is  associated  with  acti¬ 
vation  of  the  DNA  damage  response.12 
However,  senescence  induced  by  AKT 
or  loss  of  PTEN  is  independent  of  the 
DNA  damage  response.13,14  Accordingly, 
it  has  been  proposed  that  DNA  damage- 
independent  senescence  could  be  utilized 
as  a  novel  mechanism  for  developing  can¬ 
cer  therapeutics.13,15  Aberrant  DNA  repli¬ 
cation  is  thought  to  be  the  trigger  of  the 
DNA  damage  response  during  senescence 
induced  by  oncogenic  RAS  mutants.12 
Consistently,  RAS-induced  senescence  is 
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dependent  upon  S  phase  progression,  and 
inhibition  of  S  phase  progression  by  treat¬ 
ing  cells  with  Aphidicolin  or  contact  inhi¬ 
bition  blocks  RAS-induced  senescence.13 
In  contrast,  senescence  induced  by  loss  of 
PTEN  is  independent  of  S  phase  progres¬ 
sion  and  can  occur  in  quiescent  cells.13 
Thus,  the  difference  in  DNA  damage 
response  observed  in  senescence  induced 
by  different  oncogenes  is  likely  due  to 
their  ability  to  trigger  aberrant  DNA 
replication. 

Aberrant  DNA  replication  induced  by 
oncogenic  RAS  is  transient  and  occurs 
only  at  the  very  early  stages  of  senescence.12 
However,  DNA  damage  accumulates 
in  fully  senescent  cells  days  after  aber¬ 
rant  DNA  replication  in  primary  human 
cells,  which  supposedly  have  intact  DNA 
repair  machinery.12,14  This  evidence  sug¬ 
gests  that  the  impairment  of  DNA  repair 
might  contribute  to  the  accumulation  of 
DNA  damage  during  RAS-induced  senes¬ 
cence.  In  support  of  this  idea,  we  discov¬ 
ered  that  BRCA1  becomes  dissociated 
from  chromatin  in  RAS-infected  primary 
human  cells,  and  BRCA1 -mediated  DNA 
repair  response  is  impaired  in  these  cells.8 
This  observation  is  unique  in  that  other 
markers  of  DNA  damage  that  have  been 
reported  in  the  literature,  such  as  ATM, 
ATR,  Chkl,  Chk2,  yH2A  and  53BP1,  are 
all  activated  during  senescence  induced  by 
oncogenic  RAS.12,16  Importantly,  BRCA1 
chromatin  dissociation  occurs  well  before 
the  RAS-induced  cell  cycle  exit,8  sug¬ 
gesting  that  this  is  not  a  consequence  of 
the  senescence-associated  cell  cycle  exit. 
Moreover,  senescence  induced  by  AKT 
or  loss  of  PTEN  is  not  associated  with 
BRCA1  chromatin  dissociation,  correlat¬ 
ing  with  a  lack  in  DNA  damage  response 
in  these  cells.8  This  observation  further 
supports  the  hypothesis  that  BRCA1 
chromatin  dissociation  is  not  merely  a 
consequence  of  senescence.  Additionally, 
BRCA1  chromatin  dissociation  displays 
the  same  kinetics  as  DNA  damage  accu¬ 
mulation,  and  suppression  of  BRCA1 
chromatin  dissociation  is  sufficient  to 
inhibit  DNA  damage  accumulation.8 
Further,  BRCA1  knockdown  triggers 
DNA  damage  and  senescence  in  primary 
human  cells.8  Together,  these  findings 
further  support  the  premise  that  BRCA1 
chromatin  dissociation  contributes  to  the 


accumulation  of  DNA  damage  during 
RAS-induced  senescence. 

It  is  important  to  note  that  the  levels 
of  activated  RAS  oncogene  determine  the 
outcome  of  RAS  expression.  For  example, 
in  an  inducible  H-RASG12V  transgenic 
mouse  model,  low  levels  of  oncogenic 
RAS  promote  cell  proliferation,  while 
high  levels  of  oncogenic  RAS  drive  cell 
senescence.17  Interestingly,  tumorigenesis 
is  associated  with  enrichment  of  high  RAS 
induction,  senescence  and  senescence 
bypass.17  Notably,  the  levels  of  oncogenic 
H-RASG12V  observed  in  T24,  a  human 
bladder  cancer  cell  line,  is  sufficient  to 
trigger  BRCA1  chromatin  dissociation  in 
primary  human  cells.8  This  suggests  that 
physiological  level  of  RAS  observed  in 
human  cancer  cells  is  sufficient  to  trigger 
BRCA1  chromatin  dissociation. 

As  stated  previously,  oncogene-induced 
BRCA1  chromatin  dissociation  well  pre¬ 
cedes  the  cell  cycle  exit  during  senescence. 
In  addition,  oncogenic  RAS  impairs 
BRCAl-mediated  DNA  repair  response  in 
cycling  RAS-infected  cells  at  a  very  early 
stage  of  senescence.8  Together,  a  large  time 
window  is  created  for  cells  to  accumulate 
secondary  oncogenic  hits  prior  to  the 
senescence-associated  cell  cycle  exit,  which 
ultimately  leads  to  senescence  bypass  in  a 
minority  of  cells  while  the  vast  majority  of 
cells  eventually  exit  from  the  cell  cycle  and 
become  senescent.  Consistent  with  this 
model,  RAS-expressing  cells  accumulate 
significantly  greater  DNA  damage  after 
IR  treatment,  and  IR  treatment  promotes 
senescence  bypass  in  RAS-infected  cells.8 
In  mouse  models,  activated  oncogenes 
(such  as  RAS  or  BRAF)  have  been  shown 
to  initially  induce  senescence  and  forma¬ 
tion  of  benign  lesions.  However,  these 
lesions  ultimately  lead  to  the  development 
of  invasive  cancer.17'19  This  is  consistent 
with  our  model  whereby  oncogenic  RAS 
induces  senescence  in  the  vast  majority  of 
cells  while  promoting  senescence  bypass 
in  a  minority  of  cells,  which  ultimately 
outgrow  the  majority  of  senescent  cells 
and  become  transformed. 

Interestingly,  cells  from  BRCA1  exon 
11  knockout  mice  display  signs  of  cel¬ 
lular  senescence.20  Exon  11  of  the  mouse 
BRCA1  gene  encodes  for  the  two  BRCT 
repeats  of  the  BRCA1  protein.  Recently,  it 
has  been  demonstrated  that  BRCT  repeats 


are  essential  for  the  tumor  suppressive 
function  of  BRCA1.21  Notably,  we  dis¬ 
covered  that  the  C-terminus  of  BRCA1, 
which  contains  two  BRCT  repeats,  is  suf¬ 
ficient  to  be  dissociated  from  chromatin 
in  response  to  oncogenic  RAS.8  BRCT 
repeats  of  BRCA1  bind  to  phosphory- 
lated  partners  such  as  BRIP1,  Abbarax/ 
RAP80  and  CtIP  to  form  three  distinct 
complexes  in  a  mutually  exclusive  man¬ 
ner.22  These  distinct  complexes  regulate 
activation  of  checkpoints  in  response  to 
DNA  damage.22  We  found  that  BRIP1, 
but  not  RAP80  or  CtIP,  is  downregulated 
prior  to  BRCA1  chromatin  dissociation. 
In  addition,  BRIP1  knockdown  dissoci¬ 
ates  BRCA1  from  chromatin,  induces 
the  DNA  damage  response  and  triggers 
cellular  senescence.  Conversely,  ectopic 
BRIP1  suppresses  BRCA1  chromatin  dis¬ 
sociation,  the  DNA  damage  response  and 
cellular  senescence  induced  by  oncogenic 
RAS.8  Together,  these  findings  suggest 
a  key  role  of  BRIP1  in  regulating  RAS- 
induced  senescence. 

BRIP1  was  first  described  as  a  member 
of  the  DEAH  helicase  family  that  binds 
to  the  BRCT  repeats  of  BRCA1.23  The 
interaction  between  BRIP1  and  BRCA1 
depends  upon  the  phosphorylation  of 
BRIP1  at  residue  serine  990,  which  is 
regulated  in  a  cell  cycle-dependent  man¬ 
ner.24  The  BRIPl-containing  BRCA1 
complex  is  implicated  in  regulating  activa¬ 
tion  of  S  phase  checkpoint.25  For  example, 
it  has  been  demonstrated  that  BRIP1  is 
required  for  timely  S  phase  progression.25 
In  addition,  BRIP1  and  BRCA1  facilitate 
DNA  replication  during  the  S  phase  of 
the  cell  cycle  by  mediating  the  loading  of 
CDC45L,  the  replication-licensing  fac¬ 
tor.26  Thus,  RAS-induced  suppression  of 
BRIP1  may  promote  senescence  by  driving 
aberrant  DNA  replication  and  activating 
S  phase  checkpoints  during  RAS-induced 
senescence. 

The  interaction  between  BRIP1  and 
BRCA1  suggests  that  BRIP1  might  be 
linked  to  increased  cancer  risk.  Indeed,  it 
has  been  shown  that  germline  mutations 
that  alter  its  helicase  function  or  disrupt 
its  interaction  with  BRCA1  are  associated 
with  early-onset  breast  cancer.27'28  Counter 
intuitively,  the  levels  of  BRIP1  are  elevated 
in  breast  carcinoma,  and  higher  levels  of 
BRIP1  are  associated  with  higher  tumor 
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Figure  1.  B-Myb  knockdown  suppresses  the  expression  of  BRIP1  and  induces  senescence.  (A)  IMR90  primary  human  fibroblasts  were  infected  with 
lentivirus  encoding  two  individual  short  hairpin  RNAs  to  the  human  B-Myb  gene  (shB-Myb)  or  a  control.  The  sense  sequences  for  the  shB-Myb  are: 
5'-GCT  AAC  AAC  AAA  GTT  CCA  CTT-3'  and  5'-CCC  AGA  TCA  GAA  GTA  CTC  CAT-3'.  Expression  of  B-Myb  mRNA  in  the  drug-selected  cells  was  determined 
by  qRT-PCR.  (B)  Same  as  (A),  but  examined  for  BRIP1  and  p-actin  expression  by  immunoblot.  (C)  Same  as  (A),  but  equal  number  of  cells  (1  x  1 04)  were 
loaded  in  6-well  plates  for  12  d  and  stained  with  0.05%  crystal  violet  for  colony  formation.  (D)  Same  as  (A),  but  labeled  for  1  h  with  10  p,M  BrdU,  and  the 
percentage  of  the  BrdU  positive  cells  was  quantified  by  immunofluorescence  staining.  (E)  Same  as  (A),  but  stained  for  SA-p-gal  activity,  a  marker  of  cel¬ 
lular  senescence.  (F)  Quantitation  of  (D).  100  cells  from  each  of  the  indicated  groups  were  examined  for  SA-(3-gal  activity.  Mean  of  three  independent 
experiments  with  SD. 


grade.29  In  addition,  high  levels  of  BRIP1 
positively  correlate  with  HER2  status  and 
expression  of  Ki67,  a  marker  of  cell  pro¬ 
liferation.29  It  has  been  demonstrated  that 
HER2  overexpression  drives  senescence 
in  primary  human  mammary  epithelial 
cells.30  Thus,  overexpression  of  BRIP1 
may  contribute  to  breast  cancer  by  pro¬ 
moting  proliferation  of  mammary  epithe¬ 
lial  cells  through  suppressing  senescence 
induced  by  HER2  overexpression. 

As  discussed  above,  IR-treatment  in 
RAS-infected  cells  promotes  senescence 
bypass.8  In  IR-treatment-induced  senes¬ 
cence  bypassed  cells,  we  show  that  BRIP1 
remains  downregulated  and  BRCA1 


remains  dissociated  from  chromatin.8 
This  is  likely  due  to  selection  pressure 
imposed  by  exogenous  DNA  damage 
treatment.  In  addition,  we  observed  the 
accumulation  of  DNA  damage  in  these 
senescence-bypassed  cells  compared  with 
controls.8  This  suggests  that  DNA  repair 
remains  impaired  in  these  RAS-expressing 
senescence-bypassed  cells  induced  by 
IR-treatment.  Consistently,  it  has  previ¬ 
ously  been  demonstrated  that  the  DNA 
damage  response  is  attenuated  in  RAS- 
expressing  transformed  cells.31 

Interestingly,  we  observed  rare  colonies 
that  bypassed  R  AS -induced  senescence 
after  extended  culture  without  exogenous 


DNA  damage  treatment.8  The  forma¬ 
tion  of  these  colonies  may  reflect  the  rare 
senescence-bypassed  cells  that  are  associ¬ 
ated  with  endogenous  DNA  damage  trig¬ 
gered  by  oncogenic  RAS.  Thus,  it  will 
be  interesting  to  examine  the  levels  of 
BRIP1  expression  and  chromatin-associ¬ 
ated  BRCA1  levels  in  these  cells.  Because 
restoration  of  BRIP1  expression  and  sup¬ 
pression  of  BRCA1  chromatin  dissocia¬ 
tion  inhibits  RAS-induced  senescence,8 
it  is  plausible  that  the  senescence  bypass 
observed  in  these  rare  colonies  is  due  to 
restoration  of  BRIPI  expression  and,  con¬ 
sequently,  rescuing  BRCA1  chromatin 
dissociation.  Our  model  would  predict 
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that  these  senescence-bypassed  cells  with 
restored  BRIP1  expression  and  chromatin- 
associated  BRCA1  levels  would  be  inde¬ 
pendent  of  the  DNA  damage  response. 
Likewise,  it  will  be  interesting  to  deter¬ 
mine  the  levels  of  DNA  damage  in  HER2 
positive  breast  carcinoma  specimens  with 
elevated  levels  of  BRIP1.29 

RAS-induced  BRIP1  downregulation 
that  triggers  BRCA1  chromatin  disso¬ 
ciation  was  due  to  suppression  of  B-Myb 
expression,8  a  known  regulator  of  RAS- 
induced  senescence.  For  example,  ectopic 
B-Myb  suppresses  senescence  of  primary 
mouse  embryonic  fibroblasts  induced 
by  oncogenic  RAS.32  Similar  to  BRIP1, 
B-Myb  expression  is  required  for  S  phase 
entry  and  acts  as  a  crucial  factor  for  DNA 
replication  during  S  phase. 33  Depletion  of 
B-Myb  triggers  aberrant  DNA  replication 
and  the  DNA  damage  response.33  Given 
its  role  in  S  phase  entry  and  progression,  it 
is  not  surprising  that  gene  amplification  or 
overexpression  of  B-Myb  occurs  in  several 
types  of  cancer.34,35 

B-Myb  directly  regulates  BRIP1  by 
binding  to  the  proximal  promoter  region 
of  the  human  BRIP1  gene.8  However,  it 
remains  to  be  determined  whether  B-Myb 
downregualtion  is  sufficient  to  induce 
BRCA1  downregulation  and  senescence. 
To  address  this  question,  we  developed 
two  individual  short  hairpins  RNA  to 
the  human  B-Myb  gene  (shB-Myb).  The 
efficacy  of  B-Myb  knockdown  by  the 
shB-Mybs  was  confirmed  by  qRT-PCR 
(Fig.  1A).  Indeed,  BRIP1  expression 
was  downregulated  in  B-Myb  knock¬ 
down  cells  (Fig.  IB).  Further,  B-Myb 
knockdown  suppressed  the  proliferation 
of  primary  human  fibroblasts  as  deter¬ 
mined  by  decreased  colony  formation 
and  BrdU  incorporation  (Fig.  1C  and 
D).  Consistent  with  the  idea  that  growth 
inhibition  induced  by  B-Myb  knockdown 
is  due  to  senescence,  expression  of  SA-(3- 
gal,  a  marker  of  cellular  senescence,36  was 
induced  in  B-Myb  knockdown  cells  (Fig. 
IE  and  F).  Together,  these  data  further 
support  that  premise  that  B-Myb  is  a  key 
regulator  of  BRIP1  during  senescence. 

In  summary,  a  very  early  stage  of 
RAS-induced  cellular  senescence  involves 
suppression  of  BRIP1  expression  to  disso¬ 
ciate  BRCA1  from  chromatin  and  thereby 
impair  the  BRCAl-mediated  DNA  repair. 


This  promotes  cellular  senescence  and 
senescence-associated  accumulation  of 
DNA  damage  while  also  allowing  for 
subsequent  secondary  hits  that  may  pre¬ 
dispose  rare  clones  of  cells  to  escape  senes¬ 
cence  and  ultimately  contribute  to  cell 
transformation.  Therefore,  this  newly  dis¬ 
covered  cell-intrinsic  pathway  reconciles 
the  tumor-promoting  nature  of  oncogenic 
RAS  with  the  tumor-suppressing  role  of 
RAS-induced  senescence. 
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SUMMARY 

Here,  we  report  a  cell-intrinsic  mechanism  by  which 
oncogenic  RAS  promotes  senescence  while  predis¬ 
posing  cells  to  senescence  bypass  by  allowing 
for  secondary  hits.  We  show  that  oncogenic  RAS 
inactivates  the  BRCA1  DNA  repair  complex  by 
dissociating  BRCA1  from  chromatin.  This  event 
precedes  senescence-associated  cell  cycle  exit 
and  coincides  with  the  accumulation  of  DNA 
damage.  Downregulation  of  BRIP1,  a  physiological 
partner  of  BRCA1  in  the  DNA  repair  pathway,  triggers 
BRCA1  chromatin  dissociation.  Conversely,  ectopic 
BRIP1  rescues  BRCA1  chromatin  dissociation  and 
suppresses  RAS-induced  senescence  and  the  DNA 
damage  response.  Significantly,  cells  undergoing 
senescence  do  not  exhibit  a  BRCA1 -dependent 
DNA  repair  response  when  exposed  to  DNA  damage. 
Overall,  our  study  provides  a  molecular  basis  by 
which  oncogenic  RAS  promotes  senescence. 
Because  DNA  damage  has  the  potential  to  produce 
additional  “hits”  that  promote  senescence  bypass, 
our  findings  may  also  suggest  one  way  a  small 
minority  of  cells  might  bypass  senescence  and 
contribute  to  cancer  development. 


INTRODUCTION 

Activation  of  oncogenes  (such  as  RAS)  in  primary  mammalian 
cells  typically  triggers  a  cascade  of  molecular  and  cellular 
events,  which  ultimately  culminates  in  a  state  of  irreversible 
cell  growth  arrest  (Campisi,  2005).  This  process  is  termed 
oncogene-induced  senescence  and  is  an  important  tumor 
suppression  mechanism  in  vivo  (Campisi,  2005).  Paradoxically, 
the  definition  of  an  oncogene  is  a  gene  that  actively  promotes 
tumorigenesis.  The  mechanism  underlying  this  paradox  remains 
poorly  understood. 


Senescent  cells  display  several  hallmark  morphological  and 
molecular  characteristics.  These  cells  are  positive  for  senes¬ 
cence-associated  (3-galactosidase  (SA-p-gal)  activity  (Dimri 
et  al.,  1995).  In  addition,  chromatin  in  the  nuclei  of  senescent 
human  cells  often  reorganizes  to  form  specialized  domains  of 
facultative  heterochromatin  called  senescence-associated 
heterochromatin  foci  (SAHF)  (Braig  et  al.,  2005;  Narita  et  al., 
2003,  2006;  Zhang  et  al.,  2005,  2007a).  SAHF  contain  markers 
of  heterochromatin,  including  di-  and  tri-methylated  lysine  9 
histone  H3  (H3K9Me2/H3K9Me3),  histone  H2A  variant  mH2A, 
and  HMGA  (Narita  et  al.,  2003,  2006;  Zhang  et  al.,  2005).  SAHF 
formation  contributes  to  senescence-induced  cell  cycle  exit  by 
directly  sequestering  and  silencing  proliferation-promoting 
genes  (Narita  et  al.,  2003;  Zhang  et  al.,  2007a). 

Oncogene-induced  senescence  is  often  characterized  by  the 
accumulation  of  DNA  damage;  in  particular,  DNA  double-strand 
breaks  (DSBs)  (Bartkova  et  al.,  2006;  Di  Micco  et  al.,  2006).  For 
example,  oncogenic  RAS  mutants  induce  DNA  damage  by 
triggering  aberrant  DNA  replication  (Di  Micco  et  al.,  2006). 
However,  it  remains  to  be  determined  whether  impaired  DNA 
repair  contributes  to  the  accumulation  of  DNA  damage  observed 
during  oncogene-induced  senescence. 

BRCA1  plays  an  important  role  in  DNA  DSB  repair  (Scully  and 
Livingston,  2000).  Germline  mutations  in  the  BRCA1  gene 
predispose  women  to  breast  and  ovarian  cancer  (Scully  and 
Livingston,  2000),  and  inactivation  of  BRCA1  contributes  to 
cancer  development  by  causing  genomic  instability  (Turner 
et  al.,  2004).  BRCA1  interacts  with  various  DNA  damage  repair 
proteins  through  its  two  C  terminus  BRCA1  C-terminal  (BRCT) 
repeats.  The  BRCT  repeats  of  BRCA1  recognize  cognate  part¬ 
ners  by  binding  to  their  phosphoserine  residues  (Manke  et  al., 
2003;  Yu  et  al.,  2003),  and  their  binding  partners  include 
BRCA1 -interacting  protein  1  (BRIP1),  CtIP,  and  RAP80/Abraxas 
(Wang  et  al.,  2007;  Yu  et  al.,  1998,  2003).  In  addition,  BRCA1 
interacts  with  partner  and  localizer  of  BRCA2  (PALB2),  which  is 
necessary  for  localization  of  BRCA2  to  DNA  DSBs  (Xia  et  al., 
2006).  Functional  BRCA1  is  required  for  localizing/sustaining 
PALB2  at  sites  of  DNA  DSBs  and  error-free  homologous  recom¬ 
bination  repair  (Livingston,  2009;  Sy  et  al.,  2009;  Zhang  et  al., 
2009).  A  role  for  BRCA1  in  senescence  is  implied  by  findings 
from  the  BRCA1  exon  11  knockout  mouse  whose  cells  exhibit 
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Figure  1.  Oncogene-Induced  Dissociation  of  BRCA1  from  Chromatin  Occurs  prior  to  the  Oncogene-Induced  Cell  Cycle  Exit  and  Coincides 
with  the  Accumulation  of  DNA  Damage  during  Senescence 

(A)  IMR90  cells  were  infected  with  retrovirus-encoding  control  or  oncogenic  RAS  to  induce  senescence.  Drug-selected  cells  were  stained  with  DAPI  to  visualize 
SAHF  and  an  anti-BRCAl  antibody  at  day  5.  Arrows  point  to  an  example  of  BRCA1 -excluded  SAHF  in  senescent  IMR90  cells. 

(B)  Same  as  (A)  but  examined  for  expression  of  BRCA1 ,  p-actin,  and  histone  H3  in  soluble  (cytoplasmic  soluble  fraction  SN1  and  nuclear  soluble  fraction  SN2)  and 
chromatin  fractions  of  control  (C)  or  RAS  (R)-infected  IMR90  cells  by  western  blot.  Note  that  equal  amounts  of  total  proteins  were  loaded  for  total  cell  lysates,  SN1 , 
SN2,  and  chromatin  fractions. 

(C)  Same  as  (A).  At  days  1  and  6,  control  and  RAS-infected  cells  were  labeled  with  BrdU  for  1  hr  (left)  or  stained  with  DAPI  to  visualize  SAHF  and  antibodies  against 
mH2A  and  H3K9Me2  that  form  foci  that  colocalize  with  SAHF  in  senescent  human  cells  (right). 

(D  and  E)  Time  course  of  BrdU  labeling  (D)  and  formation  of  SAHF,  mH2A,  and  H3K9Me2  foci  (E)  in  control  and  RAS-infected  IMR90  cells.  Mean  of  three 
independent  experiments  with  SD. 

(F)  Same  as  (D)  but  examined  for  expression  of  BRCA1 ,  yH2AX,  and  histone  H3  in  the  chromatin  fractions  of  control  and  RAS-infected  cells  at  the  indicated  time 
points  by  western  blot. 
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signs  of  senescence  (Cao  et  al.,  2003).  These  observations 
suggest  that  senescence  and  tumorigenesis  pathways  may 
converge  on  BRCA1  -associated  DNA  damage  responses. 

Here,  we  report  a  cell-intrinsic  mechanism  by  which  onco¬ 
genic  RAS  promotes  senescence  but  at  the  same  time  predis¬ 
poses  cells  to  secondary  hits,  which  ultimately  leads  to  senes¬ 
cence  bypass. 

RESULTS 

BRCA1  Becomes  Dissociated  from  Chromatin  during 
Oncogenic  RAS-Induced  Senescence 

Senescent  cells  are  characterized  by  the  accumulation  of  DNA 
DSBs  (Bartkova  et  al.,  2006;  Di  Micco  et  al.,  2006;  Halazonetis 
et  al.,  2008),  and  one  of  the  critical  players  in  DSB  repair  is 
BRCA1  (Scully  and  Livingston,  2000).  To  test  the  hypothesis 
that  changes  in  BRCA1  function  occur  during  oncogene- 
induced  senescence,  we  first  examined  changes  in  the  subcel- 
lular  distribution  of  BRCA1  during  RAS-induced  senescence  of 
IMR90  primary  human  fibroblasts  (see  Figure  S1A  available  on¬ 
line).  BRCA1  immunofluorescence  (IF)  staining  was  performed 
in  proliferating  (control)  and  senescent  IMR90  primary  human 
fibroblasts  induced  by  RAS.  Notably,  BRCA1  was  excluded 
from  SAHF  in  senescent  cells  (Figure  1A).  In  addition,  similar 
results  were  obtained  using  multiple  anti-BRCAl  antibodies 
(one  rabbit  polyclonal  and  two  individual  mouse  monoclonal 
antibodies)  (data  not  shown).  We  next  fractionated  protein 
from  proliferating  (control)  and  senescent  IMR90  cells  into 
soluble  and  chromatin  fractions  (Mendez  and  Stillman,  2000; 
Narita  et  al.,  2006)  and  tested  each  for  the  presence  of 
BRCA1.  Compared  with  control  cells,  BRCA1  levels  were 
dramatically  decreased  in  the  chromatin  fractions  of  senescent 
IMR90  cells  (Figure  IB;  Figure  SIB). 

We  next  sought  to  exclude  the  possibility  that  BRCA1 
chromatin  dissociation  was  caused  by  supra-physiological 
levels  of  RAS.  To  do  so,  we  titrated  RAS  expression  in  IMR90 
cells  to  levels  lower  than  those  in  the  classic  bladder  cancer 
cell  line  T24  that  harbors  an  oncogenic  RAS  mutation  (H- 
RASG12V,  the  same  RAS  mutant  used  in  the  current  study)  (Hurlin 
et  al.,  1989)  (Figure  SIC).  Importantly,  we  observed  a  decrease 
in  BRCA1  levels  in  the  chromatin  fractions  of  RAS-infected 
IMR90  cells  with  lower  RAS  expression  than  the  T24  cells  (Fig¬ 
ure  SIC).  Interestingly,  BRCA1  levels  in  the  chromatin  fractions 
of  T24  cells  were  also  notably  lower  when  compared  to  IMR90 
cells,  albeit  similar  levels  of  total  BRCA1  levels  were  observed 
(Figure  SIC).  Together,  our  data  show  that  the  physiological 
levels  of  oncogenic  RAS  observed  in  cancerous  cells  are  suffi¬ 
cient  to  dissociate  BRCA1  from  chromatin  in  primary  cells. 

Next,  we  asked  whether  BRCA1  chromatin  dissociation  is 
unique  to  IMR90  cells  induced  to  senesce  by  RAS.  To  answer 


this  question,  primary  WI38  and  BJ  human  fibroblasts,  which 
both  senesce  after  oncogenic  RAS  is  expressed  (Ye  et  al., 
2007;  Zhang  et  al.,  2005)  (data  not  shown),  were  infected  with 
control  or  RAS-encoding  retrovirus.  Compared  with  controls, 
BRCA1  levels  were  dramatically  decreased  in  the  chromatin 
fractions  of  RAS-infected  WI38  and  BJ  cells  (Figure  SI  D).  Taken 
together,  these  data  demonstrate  that  BRCA1  becomes  dis¬ 
sociated  from  chromatin  during  RAS-induced  senescence. 

BRCA1  Chromatin  Dissociation  Precedes  the  Cell  Cycle 
Exit  during  RAS-induced  Senescence 

We  next  sought  to  determine  whether  BRCA1  chromatin  disso¬ 
ciation  occurs  early  or  late  during  RAS-induced  senescence. 
Toward  this  goal,  we  conducted  a  detailed  time  course  analysis 
of  chromatin-associated  BRCA1,  senescence-associated  cell 
cycle  exit  (determined  by  BrdU  incorporation  or  expression  of 
cyclin  A  or  serine  10  phosphorylated  histone  H3  [pH3S10]), 
and  other  markers  of  senescence  (such  as  formation  of  SAHF, 
mH2A  foci,  and  H3K9Me2  foci)  in  control  and  RAS-infected 
IMR90  cells  (Figures  1C-1H).  Strikingly,  as  early  as  day  1  (Fig¬ 
ure  SI  A),  BRCA1  was  largely  dissociated  from  chromatin  (Fig¬ 
ure  IF;  Figure  S1K),  which  is  well  before  the  senescence- 
associated  cell  cycle  exit  and  accumulation  of  markers  of 
senescence  (Figures  1C-1 H). 

BRCA1  forms  discrete  nuclear  foci  during  the  S/G2  phases  of 
the  cell  cycle  in  normal  cycling  cells  (Durant  and  Nickoloff, 
2005;  Scully  et  al.,  1997a;  Xu  et  al.,  2001).  Therefore,  we  sought 
to  determine  whether  RAS  expression  impairs  BRCA1  foci 
formation  in  cycling  cells.  Toward  this  goal,  at  day  2,  control 
and  RAS-infected  IMR90  cells  were  labeled  with  BrdU  to  iden¬ 
tify  S  phase  cells.  We  next  pre-extracted  soluble  proteins  from 
control  and  RAS-infected  cells  and  stained  these  cells  with 
antibodies  against  BRCA1  and  BrdU.  Indeed,  the  intensity  of 
BRCA1  foci  in  BrdU-positive  cells  was  significantly  weaker  in 
RAS-infected  IMR90  cells  compared  with  controls  (p  =  0.009) 
(Figures  II  and  1 J).  This  is  not  simply  a  consequence  of  DNA 
damage  because  this  did  not  occur  in  ionizing  radiation 
(IR)-treated  cells  (Figures  S1E  and  S1F).  Furthermore,  control 
and  RAS-infected  IMR90  cells  were  stained  with  antibodies 
against  BRCA1  and  cyclin  A,  a  marker  of  the  S/G2  phases  of 
the  cell  cycle  (Erlandsson  et  al.,  2000;  Sartori  et  at,  2007). 
Consistently,  BRCA1  foci  were  either  negative  or  notably 
weaker  in  cyclin  A-positive  RAS-infected  cells  compared  with 
controls  (Figures  S1G  and  S1H).  Finally,  FACS  analysis  re¬ 
vealed  that  RAS-infected  cells  accumulated  at  the  S  and  G2/ 
M  phases  of  the  cell  cycle  compared  with  controls  at  this  stage 
(day  2)  (Figure  IK;  Figures  SI  I  and  S1J).  From  these  results, 
we  conclude  that  oncogene-induced  dissociation  of  BRCA1 
from  chromatin  precedes  the  oncogene-induced  cell  cycle 
exit  during  senescence. 


(G  and  H)  Same  as  (D)  but  stained  for  cyclin  A  expression  (using  a  mouse  anti-cyclin  A  antibody)  (G)  or  pH3S10  (H)  by  IF  at  the  indicated  time  points. 

(I  and  J)  At  day  2,  control  and  RAS-infected  cells  were  labeled  with  BrdU  for  1  hr,  and  soluble  proteins  were  pre-extracted  (see  Experimental  Procedures  for 
details).  Pre-extracted  cells  were  stained  with  antibodies  against  BRCA1  and  BrdU.  The  intensity  of  BRCA1  staining  in  BrdU-positive  control  and  RAS-infected 
cells  was  quantified  using  MetaMorph  software  (J)  (n  =  30). 

(K)  Same  as  (D)  but  examined  for  cell  cycle  distribution  by  flow  cytometry. 

(L  and  M)  Same  as  (D)  but  examined  for  expression  of  yH2AX  (L)  or  53BP1  (M)  expression  by  IF  staining.  A  total  of  200  cells  from  each  of  the  indicated  groups 
were  examined  for  formation  of  yH2AX  or  53BP1  foci.  Mean  of  three  independent  experiments  with  SD. 

See  also  Figure  SI . 
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BRCA1  Chromatin  Dissociation  Coincides  with  the 
Accumulation  of  DNA  Damage 

We  next  asked  whether  BRCA1  chromatin  dissociation  plays 
a  role  in  the  DNA  damage  accumulation  observed  during  RAS- 
induced  senescence.  Strikingly,  accumulation  of  markers  of 
DNA  damage,  including  formation  of  yH2AX  and  53BP1  foci, 
accumulation  of  yH2AX  in  chromatin  fractions,  and  upregulation 
of  p53  displayed  the  identical  kinetics  as  dissociation  of  BRCA1 
from  chromatin  in  RAS-infected  IMR90  cells  (Figures  1 F,  1 L,  1 M; 
Figures  S1K-S10).  For  example,  as  early  as  day  1,  yFI2AX 
had  already  accumulated  in  chromatin  fractions  (Figure  IF;  Fig¬ 
ure  SI  K),  yH2AX  and  53BP1  foci  were  increased,  and  the  p53 
expression  levels  were  upregulated  in  RAS-infected  IMR90  cells 
compared  with  controls  (Figures  1L  and  1 M;  Figures  SI  L-SIO). 
We  conclude  that  RAS  expression  results  in  concomitant 
BRCA1  chromatin  dissociation  and  DNA  damage  accumulation. 

To  determine  whether  BRCA1  chromatin  dissociation  is  a 
RAS-specific  effect,  IMR90  cells  were  infected  with  control, 
RAS,  BRAF,  or  myristylated  AKT1  (myr-AKTI)-encoding  retro¬ 
virus.  Ectopic  expression  of  RAS,  BRAF,  and  myr-AKTI  was 
confirmed  by  western  blot  (data  not  shown).  Notably,  expression 
of  RAS,  BRAF,  and  myr-AKTI  all  induced  expression  of  markers 
of  senescence  in  IMR90  cells  (Figure  S2A)  (Krizhanovsky  et  al., 
2008;  Michaloglou  et  al.,  2005;  Xue  et  al.,  2007;  Zhang  et  al., 
2005).  Next,  we  examined  BRCA1  protein  levels  in  total  cell 
lysates  and  chromatin  fractions  of  control,  RAS,  BRAF,  and 
myr-AKTI -infected  IMR90  cells  by  western  blot.  Strikingly, 
compared  with  controls,  BRCA1  levels  decreased  dramatically 
in  the  chromatin  fractions  of  RAS  and  BRAF-infected  IMR90 
cells,  but  not  in  myr-AKTI -infected  IMR90  cells  (Figure  2A). 
These  results  suggest  that  BRCA1  chromatin  dissociation  is 
dependent  upon  specific  oncogenic  pathways  and  is  not  simply 
a  consequence  of  senescence. 

We  next  sought  to  determine  whether  BRCA1  chromatin 
dissociation  coincides  with  DNA  damage  accumulation  or  is 
associated  with  SAFIF  formation.  Toward  this  goal,  control, 
RAS,  BRAF,  and  myr-AKTI -infected  IMR90  cells  were  stained 
with  DAPI  to  visualize  SAFIF  and  with  an  antibody  to  mFI2A, 
which  is  a  component  of  SAFIF  (Zhang  et  al.,  2005).  Both  RAS 
and  BRAF,  but  not  myr-AKTI,  induced  formation  of  SAFIF  and 
mH2A  foci  (Figure  2B).  Additionally,  compared  with  controls, 
both  RAS  and  BRAF  expression  induced  formation  of  yH2AX 
and  53BP1  foci  and  increased  levels  of  chromatin-associated 
yFI2AX  and  53BP1  (Figures  2C-2E).  In  contrast,  myr-AKTI, 
which  did  not  induce  BRCA1  chromatin  dissociation,  failed  to 
induce  formation  of  yFI2AX  and  53BP1  foci  or  increase  the  levels 
of  chromatin-associated  yFI2AX  and  53BP1  (Figures  2C-2E). 
Consistently,  it  has  been  recently  reported  that  cell  senescence 
induced  by  myr-AKTI  is  not  associated  with  SAFIF  formation  or 
the  DNA  damage  response  (Kennedy  et  all.,  2011).  Similarly, 
PTEN  knockdown  induced  expression  of  markers  of  senescence 
but  had  no  effects  on  BRCA1  chromatin  association  and  also 
failed  to  trigger  formation  of  yH2AX  foci  or  SAFIF  (Figures 
S2B-S2G).  Likewise,  it  has  been  previously  shown  that  senes¬ 
cence  induced  by  loss  of  PTEN  is  not  associated  with  the  accu¬ 
mulation  of  DNA  damage  (Alimonti  et  al.,  201 0).  The  lack  of  both 
DNA  damage  and  BRCA1  chromatin  dissociation  in  the  PTEN 
knockdown  cells  is  consistent  with  the  idea  that  BRCA1  chro¬ 
matin  dissociation  contributes  to  DNA  damage  accumulation. 


We  next  directly  measured  the  extent  of  DNA  damage  in 
control,  RAS,  BRAF,  and  myr-AKTI -infected  IMR90  cells  using 
the  comet  assay.  Compared  with  controls,  there  was  a  significant 
increase  in  DNA  damage  in  RAS  and  BRAF-infected  IMR90 
cells  (p  <  0.05),  whereas  the  difference  between  control  and 
myr-AKTI -infected  IMR90  cells  was  not  significant  (p  >  0.05) 
(Figures  2F  and  2G).  The  combined  data  suggest  that  dissocia¬ 
tion  of  BRCA1  from  chromatin  coincides  with  DNA  damage 
accumulation,  correlates  with  SAFIF  formation,  and  is  indepen¬ 
dent  of  PI3K/AKT  signaling. 

BRCA1  Knockdown  Induces  Senescence 

We  next  sought  to  test  whether  BRCA1  knockdown  is  sufficient 
to  drive  SAFIF  formation  and  senescence  in  IMR90  cells.  Toward 
this  goal,  three  individual  short  hairpin  RNA  to  the  human 
BRCA1  (shBRCAl)  genes  with  different  degrees  of  BRCA1 
knockdown  efficacy  were  utilized.  BRCA1  knockdown  efficiency 
was  confirmed  by  IF  staining  and  western  blot  (Figures  3A  and 
3B).  Two  individual  shBRCAl  s  (#2  and  #3),  which  efficiently 
knocked  down  BRCA1 ,  induced  SAFIF  formation  and  expression 
of  SA-p-gal  activity  (Figures  3C-3F).  Notably,  an  shBRCAl  (#1) 
that  knocked  down  BRCA1  with  ~40%  efficacy  at  the  total 
protein  level  had  no  effect  on  SAFIF  formation  or  expression  of 
SA-p-gal  activity  (Figures  3A-3F).  Interestingly,  similar  to  the  hy¬ 
perproliferation  observed  in  RAS-infected  IMR90  cells  prior  to 
cell  cycle  exit  (e.g.,  Figures  1G  and  1 H;  Figure  S5)  (Di  Micco 
et  al.,  2006),  knockdown  of  BRCA1  in  IMR90  cells  also  triggers 
a  minor  but  statistically  significant  hyperproliferation  prior  to 
the  cell  cycle  exit  as  demonstrated  by  increased  BrdU  incorpo¬ 
ration  (Figure  S3A).  Consistent  with  a  previous  report  (Krum 
et  al.,  2010),  we  observed  an  increase  in  expression  of  yFI2AX 
as  well  as  formation  of  yFI2AX  foci  following  BRCA1  knockdown 
(Figures  S3B  and  S3C).  We  conclude  that  knockdown  of  BRCA1 
is  sufficient  to  drive  senescence. 

Downregulation  of  BRIP1  Expression  Triggers  BRCA1 
Chromatin  Dissociation 

Next,  we  investigated  the  molecular  mechanism  underlying 
BRCA1  chromatin  dissociation.  Notably,  a  fragment  of  the  C 
terminus  of  BRCA1  (aa  1314-1863)  containing  two  BRCT 
repeats  was  sufficient  to  be  dissociated  from  chromatin  in 
IMR90  cells  infected  with  RAS  (Figure  S4A),  suggesting  that 
factors  that  interact  with  the  BRCT  repeats  of  BRCA1  may  trigger 
BRCA1  chromatin  dissociation. 

BRCT  repeats  of  BRCA1  bind  to  BRIP1,  CtIP,  and  RAP80/ 
Abraxas  (Wang  et  al.,  2007;  Yu  et  al.,  1998,  2003).  Compared 
with  controls,  BRIP1  levels  in  both  total  cell  lysates  and  the  chro¬ 
matin  fractions  of  RAS-infected  cells  were  dramatically 
decreased  (Figures  4Aand  4B;  Figure  S4B).  In  contrast,  expres¬ 
sion  levels  of  CtIP  and  RAP80  did  not  overtly  change  in  RAS-in¬ 
fected  cells  at  the  same  time  (Figure  4A).  Notably,  RAS-induced 
downregulation  of  BRIP1  occurred  prior  to  the  cell  cycle  exit,  as 
reflected  by  the  kinetics  of  cyclin  A  and  pH3S1 0  expression  (Fig¬ 
ure  4B).  Together,  these  results  suggest  that  downregulation  of 
BRIP1  may  trigger  BRCA1  chromatin  dissociation  in  RAS-in¬ 
fected  cells. 

BRIP1  was  first  identified  as  a  BRCA1  physiological 
binding  partner  (Cantor  et  al.,  2001).  We  next  sought  to  deter¬ 
mine  whether  BRIP1  is  downregulated  at  the  mRNA  level  in 
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Figure  2.  Oncogene-Induced  BRCA1  Chromatin  Dissociation  Correlates  with  DNA  Damage  Accumulation  and  SAHF  Formation 

(A)  IMR90  cells  were  infected  with  control  or  the  indicated  activated  oncogene-encoding  retrovirus.  Expression  of  BRCA1,  histone  H3,  and  p-actin  in  the 
chromatin  fractions  and  total  cell  lysates  of  indicated  cells  was  assayed  by  western  blot  at  day  4. 

(B)  Same  as  (A)  but  stained  with  DAPI  and  an  antibody  to  mH2A,  which  is  a  component  of  SAHF.  Percentage  of  SAHF-positive  cells  is  indicated.  Mean  of  three 
independent  experiments  with  SD.  Arrows  point  to  mH2A-stained  inactivated  X  chromosome  (Costanzi  and  Pehrson,  1998). 

(C)  Same  as  (A)  but  stained  with  DAPI  or  antibodies  against  yH2AX  or  53BP1 . 

(D)  Quantitation  of  (C).  A  total  of  200  cells  were  examined  for  formation  of  yH2AX  foci  and  53BP1  foci.  Cells  with  more  than  five  nuclear  foci  for  yH2AX  or  53BP1 
were  counted  as  positive.  Mean  of  three  independent  experiments  with  SD.  *p  <  0.05  and  #p  >  0.05  as  compared  to  controls. 

(E)  Expression  of  yH2AX,  53BP1,  and  histone  H3  in  the  chromatin  fractions  of  IMR90  cells  infected  with  vector  control  or  indicated  activated  oncogenes  as 
determined  by  western  blot. 

(F)  Same  as  (A).  At  day  2,  drug-selected  indicated  cells  were  assayed  for  DNA  damage  by  the  comet  assay.  White  bars  indicate  the  examples  of  comet  tails  that 
resulted  from  damaged  DNA. 

(G)  Quantitation  of  (F).  DNA  damage  from  the  comet  assay  is  calculated  as  artificial  Olive  Moment  unit  as  described  in  Experimental  Procedures.  The  Olive 
Moment  in  100  cells  from  each  group  was  measured.  Mean  of  Olive  Moment  with  SEM  is  shown  in  blue.  *p  <  0.05  and  #p  >  0.05  as  compared  to  control. 

See  also  Figure  S2. 


RAS-infected  IMR90  cells.  Expression  of  BRIP1  and  BRCA1 
mRNA  in  control  and  RAS-infected  IMR90  cells  was  examined 
by  quantitative  reverse  transcriptase-polymerase  chain  reaction 
(qRT-PCR).  Compared  with  controls,  BRIP1  mRNA  levels  were 
significantly  decreased  as  early  as  day  1  in  RAS-infected  cells 
(Figure  4C).  However,  the  BRCA1  mRNA  expression  level  did 
not  statistically  change  at  the  same  time  (Figure  4C).  This  is 
consistent  with  the  idea  that  RAS-induced  BRIP1  repression 
triggers  BRCA1  chromatin  dissociation  prior  to  RAS-induced 
cell  cycle  exit. 

We  next  sought  to  determine  the  mechanism  underlying 
BRIP1  downregulation  in  RAS-infected  cells.  BRIP1  mRNA 
expression  is  downregulated  in  RAS-infected  cells  (Figure  4C), 


but  the  stability  of  BRIP1  mRNA  is  not  decreased  by  RAS 
expression  (data  not  shown),  suggesting  that  BRIP1  may  be 
regulated  at  the  transcriptional  level.  Consistently,  the  activity 
of  a  400  bp  (-300  bp  -  +1 00  bp)  fragment  of  the  proximal  human 
BRIP1  gene  promoter  was  significantly  suppressed  in  RAS- 
infected  cells  (Figure  S4C).  Notably,  deletion  of  a  critical 
B-Myb  binding  site  from  the  BRIP1  promoter  blocked  the  RAS- 
mediated  suppression  of  promoter  activity,  suggesting  that 
B-Myb  plays  a  critical  role  in  suppressing  BRIP1  expression  in 
response  to  RAS  (Figure  S4C).  Interestingly,  it  has  previously 
been  demonstrated  that  B-Myb  suppresses  oncogenic  RAS 
induced  senescence  (Masselink  et  al.,  2001).  Consistently,  we 
observed  downregulation  of  B-Myb  in  RAS-infected  cells  prior 
24 
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Figure  3.  BRCA1  Knockdown  Induces 
Senescence 

(A)  IMR90  cells  were  infected  with  lentivirus  en¬ 
coding  three  individual  shRNAs  to  the  human 
BRCA1  gene  or  a  vector  control  and  stained  with 
an  anti-BRCAl  antibody.  Note  that  the  BRCA1 
nuclear  foci  in  control  cells  are  associated  with  the 
S/G2  phases  of  the  cell  cycle  (Scully  et  al.,  1 997b), 
indicating  that  the  BRCA1  staining  is  specific. 

(B)  Same  as  (A)  but  assayed  for  BRCA1  and  (3-actin 
expression  by  western  blot.  The  levels  of  BRCA1 
expression  were  quantified  using  the  LI-COR 
Odyssey  imaging  system  (normalized  using 
p-actin  as  a  loading  control).  Cont,  control. 

(C)  Same  as  (A)  but  stained  for  SA-p-gal  activity  at 
day  5. 

(D)  Quantitation  of  (C).  A  total  of  100  cells  were 
examined  for  expression  of  SA-p-gal  activity. 
Mean  of  three  independent  experiments  with  SD. 
*p  <  0.03. 

(E)  Same  as  (C)  but  stained  with  DAPI  to  visualize 
SAHF  and  an  antibody  against  mH2A. 

(F)  Quantitation  of  (E).  A  total  of  100  cells  were 
examined  for  SAHF  and  mH2A  foci  formation. 
Mean  of  three  independent  experiments  with  SD. 

*p  <  0.02. 

See  also  Figure  S3. 
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to  the  RAS-induced  cells  cycle  exit  (i.e.,  as  early  as  12  hr) 
compared  with  controls  (Figure  S4D).  Significantly,  using  an 
anti-B-Myb  antibody,  chromatin  immunoprecipitation  studies 
demonstrated  that  the  binding  of  B-Myb  to  the  promoter  of  the 
human  BRIP1  gene  was  significantly  reduced  in  RAS-infected 
cells  compared  with  controls  (Figure  4D).  Additionally,  it  has 
recently  been  demonstrated  that  B-Myb  expression  is  sup¬ 
pressed  by  upregulation  of  microRNA  29  (mir29)  during  senes¬ 
cence  (Lafferty-Whyte  et  al.,  2009;  Martinez  et  al.,  201 1).  Consis¬ 
tently,  expression  of  mir29  was  upregulated  in  RAS-infected 
cells  as  early  as  6  hr  compared  with  controls  (Figure  S4E).  We 
conclude  that  RAS-mediated  inhibition  of  B-Myb  contributes  to 
the  downregulation  of  BRIP1  during  RAS-induced  senescence. 

Next,  we  asked  whether  knockdown  of  BRIP1  drives  BRCA1 
chromatin  dissociation  and  senescence.  Three  individual 
shRNAs  to  the  human  BRIP1  gene  (shBRIPI)  were  developed, 
and  the  efficacy  of  BRIP1  knockdown  was  confirmed  by  western 
blot  (Figure  4E).  BRCA1  levels  in  the  chromatin  fractions  of 
shBRIPI  -expressing  cells  were  greatly  reduced  when  compared 
to  controls  (Figure  4E).  In  addition,  we  observed  an  increase  in 
expression  of  yFI2AX  in  shBRIPI -expressing  cells  (Figure  S4F). 
Notably,  BRIP1  knockdown  induced  expression  of  markers  of 
senescence,  including  SA-p-gal  activity  and  SAFiF  formation 
(Figures  4F  and  4G).  Consistent  with  this,  BRIP1  knockdown 


increased  levels  of  pi  6,  p21 ,  p53,  and  hy- 
pophosphorylated  pRB,  demonstrating 
activation  of  both  the  p53  and  the  pRB 
pathways  (Figure  S4G).  To  limit  off-target 
effects,  we  sought  to  determine  whether 
BRIP1  restoration  in  shBRIPI -expressing 
cells  could  prevent  senescence  induced 
by  BRIP1  knockdown.  Indeed,  ectopic 
expression  of  an  shBRIPI -resistant  wild-type  BRIP1  significantly 
reduced  the  expression  of  markers  of  senescence  induced  by 
shBRIPI  (Figures  4H-4J),  indicating  the  specificity  of  the  senes¬ 
cence  phenotype  induced  by  BRIP1  knockdown.  We  conclude 
that  BRIP1  knockdown  is  sufficient  to  dissociate  BRCA1  from 
chromatin  and  induce  senescence. 

The  interaction  between  BRIP1  and  the  BRCT  repeats  of 
BRCA1  depends  on  phosphorylation  of  BRIP1  at  serine  990 
(S990)  (Yu  et  al.,  2003).  Thus,  we  determined  whether  rescuing 
of  senescence  by  BRIP1  depends  on  the  phosphorylation  of 
BRIP1  at  S990.  For  this  purpose,  we  made  a  serine  to  alanine 
mutant  (S990A)  that  mimics  the  nonphosphorylated  state  of 
BRIP1  (Yu  et  al.,  2003).  Compared  with  wild-type  BRIP1, 
BRIP1  S990A  failed  to  rescue  the  senescence  induced  by 
BRIP1  knockdown  (Figures  4PI-4J).  This  result  suggests  that 
the  interaction  between  BRCA1  and  BRIP1  plays  a  critical  role 
in  regulating  senescence. 

Ectopic  BRIP1  Rescues  BRCA1  Chromatin  Dissociation 
and  Suppresses  RAS-induced  Senescence 

We  next  asked  whether  ectopic  BRIP1  might  rescue  BRCA1 
chromatin  dissociation  and  suppress  oncogene-induced  senes¬ 
cence.  Toward  this  goal,  IMR90  cells  were  cotransduced  with 
a  retrovirus  encoding  RAS  to  induce  senescence  and  a  retrovirus 
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Figure  4.  BRIP1  Repression  Triggers  BRCA1  Chromatin  Dissociation 

(A)  IMR90  cells  were  infected  with  control  (Cont)  or  RAS-encoding  retrovirus.  Expression  of  BRIP1,  BRCA1,  RAP80,  CtIP,  and  histone  H3  in  the  chromatin 
fractions  of  control  and  RAS-infected  IMR90  cells  was  determined  by  western  blot  at  the  indicated  time  points. 

(B)  Same  as  (A).  Expression  of  BRIP1 ,  BRCA1 ,  pH3S10,  cyclin  A,  and  p-actin  was  determined  by  western  blot  at  the  indicated  time  points  in  total  cell  lysates  of 
control  and  RAS-infected  IMR90  cells. 

(C)  Same  as  (A).  Expression  of  BRCA1  and  BRIP1  mRNA  was  determined  by  qRT-PCR  at  indicated  time  points  in  control  and  RAS-infected  IMR90  cells.  Mean  of 
three  independent  experiments  with  SD.  #p  >  0.05  and  *p  <  0.01  versus  controls. 

(D)  Same  as  (A).  Control  and  RAS-infected  cells  were  subjected  to  chromatin  immunoprecipitation  using  an  anti-B-Myb  antibody  as  detailed  in  Experimental 
Procedures.  The  immunoprecipitated  DNA  was  subjected  to  quantitative  PCR  analysis  using  primers  that  cover  the  B-Myb  binding  site  in  the  promoter  of  human 
BRIP1  gene. 

(E)  IMR90  cells  were  infected  with  lentivirus-encoding  shBRIPI  or  control.  Expression  of  BRIP1,  BRCA1,  histone  H3,  and  p-actin  in  total  cell  lysates  and  in 
chromatin  fractions  was  determined  by  western  blot. 

(F)  Same  as  (E)  but  stained  for  SA-p-gal  activity  or  with  DAPI  to  visualize  SAHF. 

(G)  Quantitation  of  (F).  Mean  of  four  independent  experiments  with  SEM.  *p  <  0.05. 

(H)  IMR90  cells  were  infected  with  a  lentivirus-encoding  shBRIPI  (#3)  together  with  control,  an  shBRIPI -resistant  wild-type  BRIP1  (BRIP1  WT),  or  an  shBRIPI - 
resistant  S990A  mutant  as  detailed  in  the  Experimental  Procedures.  Expression  of  BRIP1  and  |3-actin  in  drug-selected  cells  was  determined  by  western  blot  (top). 
In  addition,  expression  of  BRIP1  mRNA  was  determined  by  RT-PCR  using  primers  designed  to  its  5'  UTR  region  (only  amplifies  endogenous  but  not  ectopic 
BRIP1  mRNA)  or  its  open  reading  frame  (amplifies  both  endogenous  and  ectopic  BRIP1  mRNA)  (bottom). 

(I)  Same  as  (H)  but  stained  for  SA-p-gal  activity  or  with  DAPI  to  visualize  SAHF. 

(J)  Quantitation  of  (I).  Mean  of  three  independent  experiments  with  SD.  *p  <  0.05  and  #p  >  0.05  compared  with  only  shBRIPI -expressing  cells. 

See  also  Figure  S4. 
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Figure  5.  Ectopic  BRIP1  Rescues  BRCA1  Chromatin  Dissociation  and  Suppresses  Oncogene-Induced  Senescence  and  DNA  Damage 
Response 

(A)  IMR90  cells  were  infected  with  retroviruses  encoding  RAS  together  with  control  (Cont),  Myc -tagged  wild  type  (WT)  BRIP1 ,  or  BRIP1  S990A  mutant.  At  day  4, 
drug-selected  cells  were  examined  for  expression  of  BRCA1 ,  BRIP1 ,  and  histone  H3  in  the  chromatin  fractions  and  for  expression  of  Myc-tagged  ectopic  BRIP1 , 
BRCA1 ,  RAS,  and  p-actin  in  total  cell  lysates. 

(B)  Same  as  (A)  but  stained  with  antibodies  against  BRCA1  and  cyclin  A. 

(C)  Quantitation  of  (B).  A  total  of  1 00  cyclin  A-positive  cells  were  examined  for  BRCA1  foci  formation.  Mean  of  three  independent  experiments  with  SD.  *p  <  0.05 
and  #p  >  0.05  compared  with  RAS-infected  controls. 

(D)  Same  as  (A)  but  stained  for  SA-p-gal  activity  or  DAPI  to  visualize  SAHF  or  an  antibody  against  yH2AX  to  detect  DNA  damage. 

(E)  Quantitation  of  (D).  Mean  of  four  independent  experiments  with  SD.  *p  <  0.015  and  #p  >  0.12  versus  RAS-infected  control  cells. 

(F)  Same  as  (A)  but  equal  numbers  of  cells  (3  x  103  cells/well)  were  plated  in  6-well  plates  in  triplicate  for  focus  formation  assays.  After  2  weeks  in  culture,  the 
plates  were  stained  with  0.05%  crystal  violet  in  PBS  to  visualize  foci.  Shown  are  representative  images  of  four  independent  experiments.  Number  of  foci  is 
indicated  as  mean  with  SD. 

(G)  Same  as  (F),  but  the  number  of  cells  was  counted  at  indicated  time  points.  Mean  of  three  independent  experiments  with  SD. 

(H)  Control  and  the  senescence-bypassed  cells  isolated  from  (G)  were  examined  for  RAS  and  p-actin  expression  by  western  blot. 

See  also  Figure  S5. 


encoding  a  Myc-tagged  wild-type  BRIP1  or  control.  Compared 
with  controls,  BRIP1  expression  notably  rescued  the  levels  of 
BRCA1  in  the  chromatin  fractions  of  RAS-infected  cells  (Figure  5A, 
lane  2  versus  lane  3).  This  was  not  due  to  a  lower  RAS  expression 
level  because  RAS  was  expressed  at  a  higher  level  in  ectopic 


BRIP1 -expressing  cells  compared  to  controls  (Figure  5A,  lane  2 
versus  lane  3).  In  addition,  the  total  BRCA1  protein  level  was 
not  increased  by  ectopic  BRIP1  (Figure  5A,  lane  2  versus  lane 
3),  implying  that  the  increased  level  of  chromatin-associated 
BRCA1  was  not  due  to  increased  levels  of  total  BRCA1 . 
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Notably,  ectopic  BRIP1  restored  BRCA1  foci  formation  in 
cyclin  A-positive  RAS-infected  cells  (Figures  5B  and  5C).  Inter¬ 
estingly,  ectopic  BRIP1  did  not  affect  RAS-induced  hyperprolif¬ 
eration  (Figure  S5),  a  trigger  of  the  DNA  damage  response 
(Di  Micco  et  al.,  2006),  suggesting  that  ectopic  BRIP1  may 
instead  affect  DNA  repair.  In  addition,  expression  of  ectopic 
BRIP1  in  RAS-infected  cells  suppressed  expression  of  markers 
of  senescence  including  SA-p-gal  activity,  SAHF  formation, 
and  senescence-induced  cell  growth  arrest,  and  inhibited 
DNA  damage  response  revealed  by  decreased  yFI2AX  foci 
formation  when  compared  to  controls  (Figures  5D-5FI).  These 
data  further  support  the  conclusion  that  BRIP1  plays  a  major 
role  in  BRCA1  chromatin  association  during  RAS-induced 
senescence. 

We  next  sought  to  determine  whether  suppression  of  BRCA1 
chromatin  dissociation  by  BRIP1  depends  on  the  phosphoryla¬ 
tion  of  BRIP1  at  S990.  Compared  with  wild-type  BRIP1,  the 
BRIP1  S990A  mutant  failed  to  rescue  BRCA1  chromatin  dissoci¬ 
ation  in  RAS-infected  cells  (Figure  5A,  lane  3  versus  lane  4)  and 
was  also  impaired  in  restoring  BRCA1  foci  in  cyclin  A-positive 
RAS-infected  cells  (Figures  5B  and  5C).  These  results  imply 
that  the  interaction  between  BRIP1  and  BRCA1  is  necessary 
for  suppression  of  BRCA1  chromatin  dissociation  by  BRIP1. 
Notably,  the  BRIP1  S990A  mutant  failed  to  suppress  senes¬ 
cence  and  its  associated  DNA  damage  accumulation  in  RAS-in¬ 
fected  IMR90  cells  (Figures  5D-5G).  We  conclude  that  suppres¬ 
sion  of  senescence  by  ectopic  BRIP1  is  dependent  on  its 
interaction  with  BRCA1 . 

Oncogenic  RAS  Impairs  the  BRCA1 -Mediated  DNA 
Repair  Response  prior  to  RAS-induced  Cell  Cycle  Exit 
during  Senescence 

We  next  asked  whether  oncogene-induced  BRCA1  chromatin 
dissociation  inactivates  the  BRCA1 -mediated  DNA  repair 
response.  Upon  DNA  damage  BRCA1  foci  largely  disappear  in 
normal  cycling  cells,  although  new  damage-induced  foci  form 
hours  after  DNA  damage  during  the  S/G2  phase  of  the  cell  cycle 
(Chen  et  al.,  1998;  Durant  and  Nickoloff,  2005;  Scully  et  al., 
1997a;  Xu  et  al.,  2001).  To  determine  the  effects  of  oncogenic 
RAS  expression  on  the  BRCA1 -mediated  DNA  repair  response, 
at  day  2,  control  and  RAS-infected  IMR90  cells  were  treated  with 
2  Gy  of  IR  to  induce  DNA  DSBs.  Notably,  IR  did  not  prevent 
BRCA1  chromatin  dissociation  in  RAS-infected  cells  (Figure  6A). 
This  was  not  due  to  a  lack  of  DNA  damage  induction  by  IR  in 
RAS-infected  cells  because  nearly  100%  of  both  control  and 
RAS-infected  cells  were  positive  for  yFI2AX  foci  (Figure  6B). 
We  next  examined  formation  of  damage-induced  BRCA1  foci 
in  the  S/G2  phases  of  cycling  cells  by  costaining  cells  with  anti¬ 
bodies  against  BRCA1  and  cyclin  A  5  hr  after  IR  treatment  (Peng 
et  al.,  2006;  Zhang  et  al.,  2009).  As  expected,  BRCA1  foci  were 
significantly  induced  upon  IR  treatment  in  controls  (Figures  6C 
and  6D).  Flowever,  formation  of  damage-induced  BRCA1  foci 
was  severely  impaired  in  RAS-infected  IMR90  cells  when 
exposed  to  IR  (Figures  6C  and  6D). 

Functional  BRCA1  is  required  for  relocating/sustaining  BRCA2 
and  its  partner  protein  PALB2  at  damage-induced  foci  (Sy  et  al., 
2009;  Zhang  et  al.,  2009),  which  are  critical  for  BRCA1 -mediated 
DNA  DSB  repair  (Fluen  et  al.,  2010).  Consistently,  formation  of 
damage-induced  BRCA2  and  PALB2  foci  was  also  significantly 


impaired  in  RAS-infected  cells  upon  IR  treatment  compared  to 
controls  (Figures  6E-6G).  Based  on  these  results,  we  conclude 
that  oncogenic  RAS  impairs  the  BRCA1 -mediated  DNA  repair 
response  prior  to  the  RAS-induced  cell  cycle  exit  during 
senescence. 

DNA  Damage  Promotes  Senescence  Bypass 
in  RAS-infected  Cells 

The  impaired  ability  of  BRCA1  to  repair  DNA  damage  promotes 
genomic  instability,  facilitates  acquisition  of  oncogenic  alter¬ 
ations,  and  ultimately  drives  tumorigenesis  (Turner  et  al.,  2004). 
Thus,  we  anticipated  that  impaired  BRCA1  -mediated  DNA  repair 
might  lead  to  DNA  damage  accumulation  and  allow  for  accumu¬ 
lation  of  secondary  hits  that  might  promote  senescence  bypass. 
Consistently,  we  reproducibly  observed  rare  foci  of  senescence- 
bypassed  cells  in  RAS-infected  IMR90  cells  (Figure  S6A).  RAS 
remains  overexpressed  in  those  cells  (Figure  S6B),  suggesting 
that  senescence  bypass  is  not  due  to  loss  of  ectopic  RAS 
expression.  To  directly  test  our  hypothesis,  control  and  RAS- 
infected  IMR90  cells  were  treated  with  2  Gy  IR  at  day  2,  and 
the  extent  of  DNA  damage  was  measured  by  the  comet  assay. 
At  this  time  point,  BRCA1  was  largely  dissociated  from 
chromatin  in  RAS-infected  proliferative  cells  (Figure  IF).  There 
was  increased  DNA  damage  in  IR-treated  RAS-infected  cells, 
which  was  significantly  greater  than  either  IR-treated  control 
cells  or  RAS-infected  cells  without  IR  treatment  (p  <  0.05) 
(Figures  7A  and  7B).  Notably,  ectopic  expression  of  wild-type 
BRIP1,  but  not  the  BRIP1  S990A  mutant,  suppressed  the  DNA 
damage  accumulation  in  IR-treated  RAS-infected  cells  (Figures 
7C  and  7D). 

We  next  sought  to  determine  whether  this  decreased  ability 
to  repair  DNA  might  lead  to  senescence  bypass.  We  tested 
this  by  focus  formation  and  cell  growth  assays.  Indeed,  IR 
treatment  consistently  induced  senescence  bypass  as  evi¬ 
denced  by  both  focus  formation  and  apparent  cell  growth  in 
IR-treated  cells  compared  to  controls  (Figures  7E  and  7F). 
Notably,  senescence-bypassed  cells  formed  colonies  under 
anchorage-independent  growth  condition  in  soft  agar  (Figures 
7G  and  7H).  As  a  negative  control,  IR  did  not  promote  the 
proliferation  of  control  cells  (Figure  7E;  Figures  S6D  and  S6E). 
To  eliminate  the  possibility  that  senescence  bypass  observed 
in  IR-treated  cells  was  due  to  loss  of  ectopic  RAS,  these  cells 
were  isolated  and  analyzed  by  western  blot  for  exogenous  RAS 
expression.  Compared  with  controls,  RAS  remained  greatly 
overexpressed  in  the  senescence-bypassed  cells  (Figure  71). 
In  addition,  pRB  was  hyperphosphorylated,  and  p53  expres¬ 
sion  was  reduced  in  the  senescence-bypassed  cells  (Figure  7J), 
suggesting  that  inactivation  of  the  key  senescence-promoting 
pRB  and  p53  pathways  contributes  to  the  senescence  bypass 
induced  by  IR  treatment.  Notably,  IR  treatment  has  no  effect  on 
senescence-associated  cell  growth  arrest  once  RAS-infected 
cells  have  exited  from  cell  cycle  (e.g.,  at  day  7)  (Figures  S6C- 
S6E).  This  result  suggests  that  the  senescence  bypass 
observed  in  IR-treated  RAS-infected  cells  is  not  due  to  the 
preexistence  of  senescence-resistant  cells.  Ectopic  BRIP1 
suppresses  senescence  (Figures  5D-5FI),  which  prevented  us 
from  determining  whether  rescuing  BRCA1  chromatin  dissocia¬ 
tion  by  ectopic  BRIP1  inhibits  IR  treatment-induced  senes¬ 
cence  bypass  in  RAS-infected  cells. 
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Figure  6.  BRCA1 -Mediated  DNA  Repair  Response  Is  Impaired  prior  to  the  Oncogene-Induced  Cell  Cycle  Exit 

(A)  IMR90  cells  were  infected  with  control  (C)  or  RAS  (R)-encoding  retrovirus.  At  day  2,  drug-selected  cells  were  treated  with  or  without  2  Gy  IR.  After  5  hr  of 
recovery,  expression  of  BRCA1 ,  histone  H3,  and  (3-actin  in  total  cell  lysates  and  in  the  chromatin  fractions  was  determined  by  western  blot. 

(B)  Same  as  (A)  but  quantified  for  yH2AX  foci  formation  in  the  indicated  groups.  Mean  of  three  independent  experiments  with  SD. 

(C)  Same  as  (A)  but  stained  with  antibodies  against  BRCA1  and  cyclin  A  (using  a  rabbit  anti-cyclin  A  antibody).  Percentage  of  cyclin  A-positive  cells  is  indicated  as 
mean  of  three  independent  experiments  with  SD. 

(D)  Quantitation  of  (C).  A  total  of  1 00  cyclin  A-positive  cells  from  each  of  the  indicated  groups  were  examined  for  BRCA1  foci  formation.  Mean  of  three  independent 
experiments  with  SD. 
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If  senescence  bypass  in  IR-treated  cells  is  achieved  through 
accumulation  of  DNA  damage  and  acquisition  of  secondary 
hits  due  to  BRCA1  chromatin  dissociation  induced  by  BRIP1 
repression,  we  anticipated  that  BRIP1  might  remain  downregu- 
lated,  and  BRCA1  might  remain  dissociated  from  chromatin 
in  the  senescence-bypassed  cells.  Indeed,  compared  with 
controls,  BRIP1  remained  downregulated,  and  BRCA1  remained 
dissociated  from  chromatin  in  the  senescence-bypassed  cells 
(Figure  7K).  These  results  further  support  the  idea  that  BRCA1 
chromatin  dissociation  is  not  merely  a  consequence  of  senes¬ 
cence-associated  cell  cycle  exit  because  the  senescence- 
bypassed  cells  are  highly  proliferative.  In  addition,  markers  of 
DNA  damage  (such  as  yH2AX  and  53BP1  foci)  were  expressed 
at  higher  levels  in  the  senescence-bypassed  cells  compared  to 
controls  (Figures  7L-7N).  This  is  consistent  with  the  idea  that 
BRCA1  chromatin  dissociation  contributes  to  DNA  damage 
accumulation  in  those  cells.  Together,  these  results  suggest 
that  the  loss  of  BRCA1 -mediated  DNA  repair  may  also  allow 
for  subsequent  hits  that  ultimately  enable  a  small  fraction  of 
Ras-induced  cells  to  bypass  senescence. 

DISCUSSION 

Our  study  reveals  that  oncogenic  RAS  induces  dissociation  of 
BRCA1  from  chromatin  prior  to  the  cell  cycle  exit  during  RAS- 
induced  senescence.  This  dissociation  of  BRCA1  from  chro¬ 
matin  coincides  with  DNA  damage  accumulation.  Downregula- 
tion  of  BRIP1,  a  BRCA1 -binding  partner,  contributes  to  BRCA1 
chromatin  dissociation.  Conversely,  ectopic  BRIP1  rescues 
BRCA1  chromatin  dissociation  and  suppresses  RAS-induced 
senescence.  Significantly,  the  BRCA1 -mediated  DNA  repair 
response  is  impaired  prior  to  the  RAS-induced  cell  cycle  exit, 
which  renders  cells  susceptible  to  the  accumulation  of 
secondary  hits.  In  some  instances  this  may  ultimately  allow  for 
senescence  bypass. 

The  Role  of  BRCA1  Chromatin  Dissociation  during 
Oncogene-Induced  Senescence 

DNA  damage  persists  in  senescent  cells  (Rodier  et  al.,  2009), 
suggesting  that  defects  in  DNA  repair  may  contribute  to  the 
accumulation  of  DNA  damage  observed  during  senescence. 
Our  data  reveal  that  BRCA1  chromatin  dissociation  coincides 
with  DNA  damage  accumulation,  which  occurs  as  early  as  the 
hyperproliferation  phase  in  RAS-infected  cells  (Figures  1F-1H, 
1L,  and  1M).  A  previous  report  by  Di  Micco  et  al.  (2006)  showed 
that  the  accumulation  of  DNA  damage  closely  follows  the  hyper¬ 
proliferation  phase  in  RAS-infected  cells  during  senescence.  The 
basis  for  this  minor  discrepancy  between  these  two  reports 
remains  to  be  determined.  It  could  be  due  to  quantitative 
approaches  used  in  the  current  study,  which  make  it  easier  to 
reveal  subtle  differences  at  early  stages  (Figures  1L  and  1M). 
We  showed  that  BRCA1  chromatin  dissociation  correlates  with 
the  DNA  damage  response  in  both  RAS  and  BRAF-infected  cells 


(Figure  2).  Furthermore,  BRCA1  knockdown  in  primary  human 
cells  triggers  the  DNA  damage  response  and  induces  senes¬ 
cence  (Figure  3;  Figure  S3B  and  S3C).  Finally,  suppression  of 
BRCA1  chromatin  dissociation  by  ectopic  BRIP1  suppresses 
the  DNA  damage  induced  by  oncogenic  RAS  (Figures  5B-5E). 
Together,  these  results  support  the  notion  that  BRCA1  chro¬ 
matin  dissociation  contributes  to  the  accumulation  of  DNA 
damage  during  oncogene-induced  senescence. 

Flerein,  we  demonstrated  that  RAS-induced  dissociation  of 
BRCA1  from  chromatin  precedes  SAFIF  formation  (Figures  1C- 
1 F),  and  BRCA1  knockdown  drives  SAFIF  formation  (Figure  3). 
Furthermore,  BRCA1  chromatin  dissociation  correlates  with 
SAFIF  formation  in  IMR90  cells  expressing  RAS  and  BRAF  onco¬ 
genes  (Figures  2A  and  2B).  Consistent  with  this  idea  that  BRCA1 
antagonizes  heterochromatin  formation  and/or  maintenance,  it 
has  been  previously  demonstrated  that  targeting  BRCA1  to  an 
amplified  lac  operator-containing  chromosome  region  in  the 
mammalian  genome  results  in  large-scale  chromatin  unfolding 
(Ye  et  al.,  2001).  We  were  unable  to  ectopically  express  wild- 
type  BRCA1  in  primary  human  cells  (data  not  shown),  which 
prevented  us  from  determining  whether  ectopically  expressed 
BRCA1  might  suppress  SAFIF  formation.  Flowever,  ectopic 
BRIP1  was  able  to  rescue  BRCA1  chromatin  dissociation  and 
suppress  SAHF  formation  (Figures  5D  and  5E).  Together,  these 
findings  support  the  idea  that  BRCA1  chromatin  dissociation 
promotes  senescence  by  contributing  to  SAHF  formation. 

Suppression  of  the  DNA  damage  response  by  ectopic  BRIP1 
inhibits  SAHF  formation  (Figures  5B-5E).  Conversely,  knock¬ 
down  of  BRCA1,  which  induces  DNA  damage,  drives  SAHF 
formation  (Figure  3;  Figure  S3).  In  addition,  AKT  or  shPTEN, 
neither  of  which  dissociates  BRCA1  from  chromatin,  also  fails 
to  induce  a  DNA  damage  response  or  SAHF  formation  (Figure  2; 
Figure  S2).  Together,  these  data  suggest  that  DNA  damage 
response  triggered  by  BRCA1  chromatin  dissociation  is  required 
for  SAHF  formation.  Indeed,  there  is  evidence  to  suggest  that 
formation  of  SAHF  limits  the  degree  of  DNA  damage  response 
during  oncogene-induced  senescence  (Di  Micco  et  al.,  2011). 
However,  the  DNA  damage  response  is  not  sufficient  for  SAHF 
formation,  which  also  requires  activation  of  p16/pRB  and 
HIRA/PML  pathways  (Narita  et  al.,  2003;  Ye  et  al.,  2007;  Zhang 
et  al.,  2005,  2007a).  Overall,  these  results  support  the  notion 
that  DNA  damage  is  necessary  but  not  sufficient  for  SAHF 
formation. 

The  Role  of  BRIP1  Repression  during 
Oncogene-Induced  Senescence 

Stable  BRIP1  knockdown  reduces  damage-induced  BRCA1  foci 
formation,  and  BRIP1 -deficient  cells  demonstrate  defects  in  the 
number  and  intensity  of  BRCA1  foci  (Peng  et  al.,  2006).  This 
suggests  that  BRIP1  plays  an  important  role  in  BRCA1 -mediated 
DNA  repair.  BRIP1  has  been  shown  to  be  significantly  downre¬ 
gulated  in  response  to  RAS  expression  in  primary  human  fibro¬ 
blasts  based  on  gene  expression  microarray  analysis,  while 


(E)  Same  as  (C)  but  stained  with  antibodies  against  BRCA1  and  BRCA2. 

(F)  Same  as  (C)  but  stained  with  antibodies  against  BRCA1  and  PALB2. 

(G)  Quantitation  of  results  from  (E)  and  (F).  A  total  of  1 00  cells  from  each  of  the  indicated  group  were  examined  for  BRCA1 ,  BRCA2,  and  PALB2  foci  formation. 
Mean  of  three  independent  experiments  with  SD.  *p  <  0.02  versus  control  IR-treated  cells.  Cont,  control. 
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Figure  7.  DNA  Damage  Promotes  Senescence  Bypass 

(A)  IMR90  cells  were  infected  with  control  or  RAS-encoding  retrovirus,  and  drug-selected  cells  were  treated  with  or  without  2  Gy  IR  at  day  2.  After  5  hr  of  recovery, 
comet  assays  were  performed  to  detect  DNA  damage.  Shown  are  representative  images  of  six  independent  experiments.  The  white  bar  indicates  an  example  of 
the  comet  tails  that  reflect  the  extent  of  DNA  damage.  Note  that  a  shorter  exposure  time  was  used  here  compared  with  Figure  2F  to  avoid  saturation  of  comet  tail 
signal  in  IR-treated  RAS-infected  cells. 

(B)  Quantitation  of  results  from  (A).  A  total  of  1 20  cells  from  each  of  the  indicated  groups  were  measured  for  Olive  Moment.  Means  of  Olive  Moment  with  SEM  are 
indicated  in  blue.  Note  that  there  was  a  significant  increase  in  DNA  damage  in  RAS-infected  cells  without  IR  treatment  compared  with  controls  (p  <  0.05). 

(C)  Same  as  (A).  Comet  assays  were  performed  to  detect  DNA  damage  in  the  indicated  cells.  Shown  are  representative  images  of  three  independent  experiments. 
Cont,  control. 

(D)  Quantitation  of  results  from  (C).  A  total  of  1 00  cells  from  each  of  the  indicated  groups  were  measured  for  Olive  Moment.  Means  of  Olive  Moment  with  SEM  are 
indicated  in  blue. 

(E)  At  day  2,  control  and  RAS-infected  BJ-hTERT  cells  were  treated  with  or  without  2  Gy  IR.  Equal  numbers  of  the  indicated  cells  (3  x  103  cells/well)  were  plated  in 
6-well  plates  for  focus  formation  assays.  After  2  weeks  of  culture,  the  plates  were  stained  with  0.05%  crystal  violet  in  PBS  to  visualize  foci.  Shown  are  repre¬ 
sentative  images  of  four  independent  experiments.  The  number  of  foci  is  indicated  as  mean  with  SD. 
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CtIP  and  RAP80  are  not  among  the  list  of  genes  whose  expres¬ 
sion  is  significantly  changed  in  the  same  analysis  (Mason  et  al., 
2004).  Interestingly,  BRIP1,  but  not  CtIP  or  RAP80/Abraxas, 
regulates  DNA  replication  and  is  required  for  timely  S  phase 
progression  (Huen  et  al.,  2010;  Kumaraswamy  and  Shiekhattar, 
2007),  implying  that  the  BRCA1/BRIP1  complex  may  regulate 
senescence  by  altering  DNA  replication.  In  support  of  this  idea, 
RAS-induced  senescence  is  characterized  as  a  DNA  damage 
response  triggered  by  aberrant  DNA  replication  (Di  Micco 
et  al.,  2006).  Furthermore,  B-Myb,  a  key  regulator  of  BRIP1 
expression  (Figure  4D;  Figures  S4C  and  S4D),  is  critical  for 
proper  DNA  replication  and  regulates  the  DNA  damage  response 
(Ahlbory  et  al.,  2005;  Lorvellec  et  al.,  201 0).  Notably,  BRIP1  is  ex¬ 
pressed  at  higher  levels  in  advanced  breast  carcinomas,  and  its 
overexpression  correlates  with  a  higher  cell  proliferation  index 
(Eelen  et  al.,  2008),  supporting  the  notion  that  high  levels  of 
BRIP1  may  contribute  to  cell  proliferation  by  suppressing  senes¬ 
cence  (Figure  5).  Taken  together,  these  data  imply  a  key  role  for 
BRIP1  in  accumulation  of  DNA  damage  during  oncogene- 
induced  senescence. 

The  Role  of  BRCA1  Chromatin  Dissociation  in 
Senescence  Bypass 

Oncogene-induced  BRCA1  chromatin  dissociation  precedes 
the  cell  cycle  exit  during  senescence  (Figure  1).  Significantly, 
oncogenic  RAS  impairs  BRCA1 -mediated  DNA  repair  response 
prior  to  the  cell  cycle  exit  during  senescence  (Figures  6  and  7). 
This  allows  for  the  creation  of  a  large  time  window  for  cells  to 
accumulate  secondary  oncogenic  hits  prior  to  the  senescence- 
associated  cell  cycle  exit,  which  ultimately  leads  to  senescence 
bypass  in  a  minority  of  cells  while  the  vast  majority  of  cells 
eventually  exit  from  the  cell  cycle  and  become  senescent. 
Consistent  with  this  model,  RAS-expressing  cells  accumulate 
significantly  greater  DNA  damage  after  IR  treatment,  and  IR 
treatment  promotes  senescence  bypass  in  RAS-infected  cells 
(Figure  7).  The  mechanism  we  uncovered  here  may  help  explain 
the  paradox  of  why  activation  of  oncogenes  (such  as  RAS) 
promotes  senescence  but  at  the  same  time  predisposes  cells 
to  transformation. 

EXPERIMENTAL  PROCEDURES 

Chromatin  Isolation  and  Chromatin  Immunoprecipitation 

Chromatin  was  prepared  according  to  published  methods  (Mendez  and 
Stillman,  2000;  Narita  et  al.,  2006).  Soluble  proteins  in  supernatant  1  (SN1, 
cytoplasmic)  and  2  (SN2,  nuclear)  and  chromatin-bound  proteins  in  the  chro¬ 
matin  fraction  were  detected  by  western  blot. 


Chromatin  immunoprecipitation  in  control  and  RAS-infected  IMR90  cells 
was  performed  at  day  6  as  previously  described  (Zhang  et  al.,  2007b)  using 
a  monoclonal  anti-B-Myb  antibody  (Santa  Cruz  Biotechnology)  or  an  iso- 
type-matched  IgG  control.  Immunoprecipitated  DNA  was  analyzed  using 
SYBR  Green  quantitative  PCR  (SABiosciences)  against  the  human  BRIP1 
gene  promoter  region  containing  the  B-Myb  binding  site  using  the  following 
primers:  forward,  5'-ATAAAGCGGAGCCCTGGAAGAGAA-3';  reverse,  5'- 
ATT  CGT  CT  CGGGTT  GT  GTG-GTT  G  A-3' . 

Comet  Assay 

The  comet  assay  was  performed  with  the  CometAssay  (Trevigen)  kit  following 
the  manufacturer’s  instructions.  DNA  damage  was  measured  as  the  artificial 
Olive  Moment  using  CometScore  software  downloaded  from  http://www. 
tritekcorp.com.  To  determine  significance,  the  t  test  was  performed  using 
GraphPad  Prism  software  (http://www.graphpad.com). 

Anchorage-Independent  Growth  in  Soft  Agar  and  Focus  Formation 
Assay 

Anchorage-independent  growth  in  soft  agar  was  performed  as  previously 
described  (Li  et  al.,  2010).  For  focus  formation,  4  days  after  initial  infection, 
control  or  oncogenic  H-RASG12V-encoding  retrovirus-infected  BJ-hTERT  cells 
were  treated  with  2  Gy  IR.  Cells  were  cultured  for  6  days  to  eliminate  the 
apoptotic  cells  induced  by  IR.  Then,  control  and  H-RASG12V-expressing  cells 
with  or  without  IR  treatment  were  seeded  into  6-well  plates  at  a  density  of 
3000  cells/well  in  triplicate.  Two  weeks  later,  the  plates  were  stained  with 
0.05%  crystal  violet  in  PBS  (Kuilman  et  al.,  2008). 

Retrovirus  and  Lentivirus  Infections 

Retrovirus  production  and  transduction  were  performed  as  described  previ¬ 
ously  (Ye  et  al.,  2007;  Zhang  et  al.,  2005)  using  Phoenix  cells  to  package  the 
infection  viruses  (Dr.  Gary  Nolan,  Stanford  University).  Lentivirus  was  pack¬ 
aged  using  the  ViraPower  Kit  from  Invitrogen  following  the  manufacturer’s 
instructions  and  as  described  previously  (Li  et  al.,  2010;  Ye  et  al.,  2007).  Cells 
infected  with  viruses  encoding  a  drug-resistant  gene  to  puromycin  or  neomycin 
were  selected  in  1  and  500  ng/ml,  respectively,  of  the  corresponding  agent. 

Immunofluorescence,  BrdU  Labeling,  FACS,  and  SA-p-gal  Staining 

Immunofluorescence  staining  and  BrdU  labeling  for  cultured  cells  were  per¬ 
formed  as  described  previously  using  antibodies  described  above  (Zhang 
et  al.,  2005,  2007a,  2007b).  Soluble  protein  pre-extraction  with  detergent 
was  carried  out  as  described  previously  (Taddei  et  al.,  2001).  Briefly,  cells 
were  incubated  with  PBS  supplemented  with  0.5%  Triton  X-100  for  5  min  at 
room  temperature,  followed  by  fixation  in  4%  paraformaldehyde  (Sigma- 
Aldrich)  for  10  min.  Fixed  cells  were  incubated  with  a  rabbit  anti-BRCAl 
antibody  for  1  hr  at  room  temperature  and  visualized  by  incubating  the  cells 
with  goat  anti-rabbit  Cy3  (Jackson  ImmunoResearch;  1:5000)  secondary 
antibody  followed  by  detection  of  BrdU  using  a  FITC-labeled  anti-BrdU  anti¬ 
body  (BD  Biosciences).  FACS  was  performed  as  previously  described  (Ye 
et  al.,  2007),  and  FlowJo  software  was  used  to  analyze  cell  cycle  distribution. 
SA-p-gal  staining  was  performed  as  previously  described  (Dimri  et  al.,  1995). 

For  additional  information  about  cell  culture  methods,  RT-PCR,  and  lucif- 
erase  assays,  as  well  as  antibodies  and  plasmids  used  in  this  study,  please 
see  the  Supplemental  Experimental  Procedures. 


(F)  Same  as  (E),  but  an  equal  number  (1x1 04)  of  indicated  cells  were  inoculated  and  cultured.  The  number  of  cells  from  each  group  was  counted  at  indicated  time 
points.  Mean  of  three  independent  experiments  with  SD. 

(G)  Equal  numbers  of  parental  and  senescence-bypassed  cells  were  grown  under  anchorage-independent  conditions  in  soft  agar.  Arrow  points  to  an  example  of 
colonies  formed  by  the  senescence-bypassed  cells. 

(H)  Quantitation  of  (G).  Mean  of  three  independent  experiments  with  SD. 

(I)  Expression  of  RAS  and  p-actin  in  control  (C)  and  the  senescence-bypassed  (SB)  cells  isolated  from  (F)  was  determined  by  western  blot. 

(J)  Same  as  (I)  but  assayed  for  pRB,  p53,  and  p-actin  expression  by  western  blot. 

(K)  Same  as  (I)  but  examined  for  expression  of  BRIP1 ,  BRCA1 ,  and  histone  H3  in  total  cell  lysates  and  in  the  chromatin  fractions  by  western  blot. 

(L)  Same  as  (I)  but  stained  with  an  anti-yH2AX  antibody. 

(M)  Same  as  (L)  but  stained  with  an  anti-53BP1  antibody.  Arrows  point  to  examples  of  53BP1  foci-positive  cells. 

(N)  Quantitation  of  (L)  and  (M).  A  total  of  100  cells  from  each  indicated  group  were  examined  for  yH2AX  and  53BP1  foci  formation.  Mean  of  three  independent 
experiments  with  SD. 

See  also  Figure  S6. 
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SUPPLEMENTAL  INFORMATION 

Supplemental  Information  includes  six  figures  and  Supplemental  Experimental 
Procedures  and  can  be  found  with  this  article  online  at  doi:10.1016/j.devcel. 
2011.10.010. 
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Abstract 

Epithelial  ovarian  cancer  (EOC)  remains  the  most  lethal  gynecologic  malignancy  in  the  United  States. 
EZH2  silences  gene  expression  through  trimethylating  lysine  27  on  histone  H3  (I  I3K27Me3).  EZH2  is  often 
overexpressed  in  EOC  and  has  been  suggested  as  a  target  for  EOC  intervention.  However,  EZH2  target  genes 
in  EOC  remain  poorly  understood.  Here,  we  mapped  the  genomic  loci  occupied  by  EZH2/EI3K27Me3  using 
chromatin  immunoprecipitation  followed  by  next-generation  sequencing  (ChIP-seq)  and  globally  profiled 
gene  expression  in  EZH2-knockdown  EOC  cells.  Cross-examination  of  gene  expression  and  ChIP-seq 
revealed  a  list  of  60  EZH2  direct  target  genes  whose  expression  was  upregulated  more  than  1.5-fold  upon 
EZH2  knockdown.  For  three  selected  genes  ( ALDH1A1 ,  SSTR1,  and  DACT3),  we  validated  their  upregula- 
tion  upon  EZH2  knockdown  and  confirmed  the  binding  of  EZH2/H3K27Me3  to  their  genomic  loci. 
Furthermore,  the  presence  of  H3K27Me3  at  the  genomic  loci  of  these  EZEI2  target  genes  was  dependent 
upon  EZH2.  Interestingly,  expression  of  ALDH1A1,  a  putative  marker  for  EOC  stem  cells,  was  significantly 
downregulated  in  high-grade  serous  EOC  (n  =  53)  compared  with  ovarian  surface  epithelial  cells  (n  =  10, 
P  <  0.001).  Notably,  expression  of  ALDII1A1  negatively  correlated  with  expression  of  EZH2  ( n  =  63, 
Spearman  r  =  —0.41,  P  <  0.001).  Thus,  we  identified  a  list  of  60  EZH2  target  genes  and  established  that 
ALDH1A1  is  a  novel  EZIT2  target  gene  in  EOC  cells.  Our  results  suggest  a  role  for  EZEI2  in  regulating  EOC 
stem  cell  equilibrium  via  regulation  of  ALDEI1A1  expression.  Cancer  Prev  Res;  5(3);  484-91.  ©201 1  AACR. 


Introduction 

Epithelial  ovarian  cancer  (EOC)  accounts  for  more  deaths 
than  any  other  gynecologic  malignancy  in  the  United  States. 
EOCs  are  classified  into  distinct  histologic  types  including 
serous,  mucinous,  endometrioid,  and  clear  cell  (1).  The 
most  common  histology  of  EOC  is  serous  (~60%  of  all 
cancers)  and  less  common  histologies  include  endome¬ 
trioid,  clear  cell,  and  mucinous  (1).  Recently,  an  alternative 
classification  has  been  proposed,  in  which  EOC  is  broadly 
divided  into  2  types  (2).  Type  I  EOC  includes  mucinous, 
low-grade  serous,  low-grade  endometrioid,  and  clear  cell 
carcinomas,  and  type  II  EOC  includes  high-grade  serous 
carcinomas,  which  is  the  most  lethal  histosubtype  (2). 
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Enhancer  of  zeste  homology  2  (EZH2)  is  a  histone 
methyltransferase  that  mediates  gene  silencing  by  catalyz¬ 
ing  the  trimethylation  on  lysine  27  of  histone  H3 
(IT3K27Me3;  ref.  3).  EZEI2  is  often  expressed  at  higher 
levels  in  human  EOC  cells,  and  its  expression  positively 
correlates  with  cell  proliferation  in  these  cells  (4).  Further 
underscoring  the  importance  of  EZH2  in  EOC,  EZH2 
knockdown  triggers  apoptosis  and  inhibits  the  invasion  of 
human  EOC  cells  (4).  In  addition,  EZH2  is  overexpressed  in 
ovarian  tumor-associated  endothelial  cells,  which  pro¬ 
motes  angiogenesis  (5).  Finally,  there  is  evidence  to  suggest 
that  EZH2  is  overexpressed  in  ovarian  cancer  stem  cell-like 
populations  enriched  by  chemotherapy  (6).  Accordingly, 
EZTI2  has  been  suggested  as  a  putative  target  for  developing 
EOC  therapeutics.  Thus,  it  is  important  to  identify  EZH2 
target  genes  in  EOC  to  gain  insights  into  the  biology  of  the 
disease  and  to  facilitate  translational  EOC  research  related 
to  EZH2.  Although  a  number  of  EZEI2  target  genes  have 
been  characterized  in  a  few  cancer  types,  including  prostate 
and  breast,  using  chromatin  immunoprecipitation  (ChlP)- 
on-chip  analysis  (7,  8),  studies  that  aim  to  globally  identify 
EZTI2  target  genes  in  EOC  cells  have  yet  to  be  conducted. 

Here,  we  report  the  identification  of  direct  EZTI2  target 
genes  in  human  EOC  cells  using  a  combination  of  genome¬ 
wide  approaches.  Specifically,  we  identified  the  genomic 
loci  occupied  by  EZH2/H3K27Me3  using  ChIP  followed  by 
next-generation  sequencing  (ChIP-seq).  In  addition,  we 
discovered  a  list  of  genes  whose  expression  was  upregulated 


484  Cancer  Prev  Res;  5(3)  March  201 2 


^^C^American  Association  for  Cancer  Research 


36 

Downloaded  from  cancerpreventionresearch.aacrjournals.org  on  October  4,  2012 
Copyright  ©  201 2  American  Association  for  Cancer  Research 


Published  OnlineFirst  December  5,  201 1 ;  DOI:1 0.1 158/1 940-6207.CAPR-1 1-0414 


EZH2  Targets  ALDH1A1  in  EOC 


more  than  1. 5-fold  in  EZH2-knockdown  EOC  cells  com¬ 
pared  with  controls  using  gene  expression  microarray  anal¬ 
ysis.  Cross-examination  of  gene  expression  profiling  and 
ChIP-seq  analysis  revealed  a  list  of  60  genes  that  are  direct 
EZH2/H3K27Me3  target  genes,  including  56  novel  putative 
EZH2  target  genes.  For  validation,  we  selected  3  genes  that 
are  implicated  in  regulating  stem  cells,  apoptosis,  cell 
growth,  or  invasion.  We  validated  their  upregulation  upon 
EZH2  knockdown  in  EOC  cells  and  confirmed  the  binding 
ofEZH2/H3K27Me3  by  ChIP  analysis.  Interestingly,  expres¬ 
sion  ofALDHIAl,  a  putative  marker  for  EOC  stem  cells  (Si¬ 
ll),  was  expressed  at  significantly  lower  levels  in  high-grade 
serous  EOC  than  in  normal  human  ovarian  surface  epithe¬ 
lial  (HOSE)  cells  and  negatively  correlated  with  expression 
of  EZH2. 

Materials  and  Methods 

Cell  culture,  short  hairpin  RNA,  lentivirus  packaging 
and  infection 

The  SKOV3  human  EOC  cell  line  was  cultured  according 
to  the  American  Type  Culture  Collection  and  as  previously 
described  (4, 12).  SKOV3  cell  line  identification  was  further 
confirmed  by  DNA  Diagnostic  Center.  The  sense  sequences 
of  2  individual  short  hairpin  RNAs  (shRNA)  to  the  human 
EZH2  genes  are  as  we  have  previously  published  (4).  Len¬ 
tivirus  packaging  was  conducted  using  ViraPower  system 
(Invitrogen)  according  to  the  manufacturer's  instruction 
and  as  previously  described  (4).  Briefly,  SKOV3  cells  at 
40%  to  50%  confluency  were  infected  with  lentivirus 
expressing  shRNA  to  EZH2  or  vector  control.  The  infected 
cells  were  drug  selected  with  3  pg/mL  of  puromycin  to 
eliminate  noninfected  cells. 

Antibodies,  Western  blot  analysis,  RNA  isolation,  and 
quantitative  reverse-transcriptase  PCR 

The  following  antibodies  were  used  for  Western  blot 
analysis:  mouse  anti-EZH2  (1:2,500;  BD  Bioscience),  rabbit 
anti-H3K27Me3  (1:1,000;  Cell  Signaling),  and  mouse  anti- 
GAPDII  (1:10,000;  Millipore).  RNA  from  cultured  human 
SKOV3  EOC  cells  was  isolated  using  TRIzol  (Invitrogen) 
according  to  the  manufacturer's  instruction.  For  quantita¬ 
tive  reverse-transcriptase  PCR  (qRT-PCR),  TRIzol-isolated 
RNA  was  further  purified  using  an  RNeasy  kit  (Qiagen) 
following  the  manufacturer's  instruction.  The  primers  for 
ALDH1A1,  SSTR1,  and  DACT3  genes  used  for  qRT-PCR 
were  purchased  from  Applied  Biosystems.  Expression  of 
the  housekeeping  gene  ji-2-microglobulin  was  used  to  nor¬ 
malize  mRNA  expression. 

ChIP-seq  analysis  and  ChIP  validation  for  selected 
EZH2  target  genes 

Briefly,  SKOV3  cells  were  fixed  with  1%  formaldehyde  for 
1 5  minutes  and  quenched  with  0.125  mol/L  glycine.  Chro¬ 
matin  was  isolated  by  adding  lysis  buffer  (1%  SDS,  10 
mmol/L  EDTA,  50  mmol/L  Tris-HCl,  pH  8.1,  1  mmol/L 
phenylmethylsulfonylfluoride)  followed  by  disruption 
with  a  Dounce  homogenizer.  Lysates  were  sonicated  using 


a  Misonix  Sonicator  3000  to  shear  the  DNA  to  an  average 
length  of  300  to  500  bp.  Lysates  were  cleared  by  centrifu¬ 
gation  to  collect  chromatin  suspensions.  Prior  to  their  use  in 
the  ChIP  protocol,  protein  A  agarose  beads  (Invitrogen) 
were  preblocked  using  blocking  proteins  and  nucleic  acids 
for  3  hours.  For  each  ChIP  reaction,  an  aliquot  of  chromatin 
(20-30  pg)  was  precleared  with  30  pL  preblocked  protein  A 
agarose  beads  for  1  to  2  hours.  ChIP  reactions  were  set  up 
using  precleared  chromatin  and  antibody  (anti- 
H3K27Me3,  Millipore  07-449;  anti-EZH2,  Millipore  07- 
689)  and  incubated  overnight  at  4°C.  Preblocked  protein 
A  agarose  beads  were  added  and  incubation  at  4°C  was 
continued  for  another  3  hours.  Agarose  beads  containing 
the  immune  complexes  were  washed,  and  the  immune 
complexes  eluted  from  the  beads  were  subjected  to  RNase 
treatment  at  37°C  for  20  minutes  and  proteinase  K  treat¬ 
ment  at  37°C  for  3  hours.  Cross-links  were  reversed,  and 
ChIP  DNAs  were  purified  by  phenol-chloroform  extraction 
and  ethanol  precipitation. 

ChIP  DNA  was  amplified  using  the  Ulumina  ChIP-Seq 
DNA  Sample  Prep  Kit.  In  brief,  DNA  was  resonicated  and 
ends  were  polished  and  5'-phosphorylated  using  T4  DNA 
polymerase,  Klenow  polymerase,  and  T4  polynucleotide 
kinase.  After  addition  of  3'-A  to  the  ends  using  Klenow 
fragment  (3'-5'  exo  minus),  Illumina  genomic  adapters 
were  ligated  and  the  sample  was  size-fractionated  (300- 
400  bp)  on  a  2%  agarose  gel.  After  a  final  PCR  amplification 
step  (18  cycles,  Phusion  polymerase),  the  resulting  DNA 
libraries  were  quantified  and  tested  by  qPCR  at  the  same 
specific  genomic  regions  as  the  original  ChIP  DNA  to  assess 
quality  of  the  amplification  reactions.  DNA  libraries  were 
sequenced  on  a  Genome  Analyzer  II.  Sequences  (36-nucle¬ 
otide  reads)  were  aligned  to  the  human  genome  (NCBI 
Build  37.1/hgl9)  using  Eland  software  (Illumina).  Aligned 
sequences  were  extended  in  silico  at  their  3'-ends  to  a  length 
of  240  bp,  which  is  the  average  genomic  fragment  length  in 
the  size-selected  library,  and  assigned  to  32-nucleotide  bins 
along  the  genome.  The  resulting  histograms  were  stored  in 
BAR  (Binary  Analysis  Results)  files.  Peak  locations  were 
determined  using  the  MACS  algorithm. 

For  validation  of  binding  of  EZH2/H3K27Me3  to  the 
genomic  loci  of  the  selected  EZFI2/FI3K27Me3  target  genes, 
SKOV3  EOC  cells  were  transduced  with  lentivirus  encoding 
control  or  shEZH2.  Drug-selected  cells  were  subjected  to 
ChIP  analysis  as  previously  described  (13,  14).  The  follow¬ 
ing  antibodies  were  used  to  conduct  ChIP:  anti-EZIT2  (Cl  1, 
BD  Biosciences),  anti-H3K27Me3  (C36B11,  Cell  Signal¬ 
ling),  and  anti-histone  H3  (05-928,  Millipore).  An  iso¬ 
type-matched  IgG  was  used  a  negative  control.  Immuno- 
precipitated  DNA  was  analyzed  with  PCR  against  the  geno¬ 
mic  regions  ofALDHIAl  (forward:  5'-TGGCACTGGTTATT- 
CAACGTGGTC-3'  and  reverse:  5'-GAGGGTGGAAGCT- 
CTTGTAGGTTT-3'),  DACT3  (forward:  5' -  C AC  AC  AC  AC A- 
CACAAACAGTGCCT-3'  and  reverse:  5'-TTCCTCCAAC- 
TAGGCTGGCAGTTT-3')  and  SSTR1  (forward:  5'-TAGCC- 
TAAGCTGCCTGCTGTGTTA-3'  and  reverse:  5'-AAAGTG- 
CATGTGCGGTCTGTTAGC-3').  PCR  products  were  visual¬ 
ized  on  a  2%  agarose  gel. 
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Gene  expression  microarray  analysis 

For  gene  expression  microarray  analysis  in  SKOV3  cells, 
500  ng  of  total  RNAwas  amplified  and  labeled  using  Agilent 
QuickAmp  labeling  kit  following  the  manufacturer's  pro¬ 
tocol.  A  total  of  1.65  |ig  of  Cy-3-labeled  cRNA  targets  were 
hybridized  onto  Agilent  4  x  44  k  whole  genome  arrays  for 
17  hours  at  65°C  and  washed  according  to  procedure 
described  by  Agilent.  The  hybridized  slides  were  scanned 
at  5 -pm  resolution  on  an  Agilent  scanner  (Agilent),  and 
fluorescent  intensities  of  hybridization  signals  are  extracted 
using  Agilent  Feature  Extraction  software. 

Data  sets 

Gene  expression  microarray  data  sets  for  53  cases  of  laser 
capture  and  microdissected  (LCM)  high-grade  serous  EOC 
and  10  individual  isolations  of  normal  HOSE  cells  were 
obtained  from  Gene  Expression  Omnibus  (GEO;  http:// 
www.ncbi.nhn.nih.gov/geo/;  GEO  accession  number: 
GSE18521). 


Statistical  analysis 

Quantitative  data  are  expressed  as  mean  ±  SD  unless 
otherwise  stated.  AN OVA  with  Fisher  least  significant  dif¬ 
ference  (LSD)  was  used  to  identify  significant  differences  in 
multiple  comparisons.  Spearman  test  was  used  to  measure 
statistical  dependence  between  EZH2  mRNA  levels  and 
ALDI 1 1A1  mRNA  levels.  For  all  statistical  analyses,  the  level 
of  significance  was  set  at  0.05. 

Results  and  Discussion 

Genome-wide  mapping  of  EZH2/H3K27Me3  direct 
target  genes  in  human  EOC  cells 

To  identify  genes  whose  expression  was  suppressed  by 
EZFI2,  we  conducted  gene  expression  microarray  analysis  in 
control  and  EZH2-knockdown  SKOV3  human  EOC  cells. 
Two  individual  shRNAs  to  the  human  EZH2  gene  (shEZH2) 
were  used  to  limit  potential  off-target  effects.  Knockdown 
efficacy  was  confirmed  by  immunoblotting  analysis 
(Fig.  1A).  EZH2  knockdown  notably  decreased  the  levels 
of  FI3K27Me3,  which  is  consistent  with  the  idea  that  EZH2 
plays  a  major  role  in  regulating  the  levels  of  H3K27Me3  in 
human  EOC  cells  (Fig.  1  A).  A  list  of  148  genes  and  501  genes 
were  upregulated  more  than  1.5-fold  by  shEZH2#l  and 
shEZH2#2,  respectively,  whereas  128  genes  overlapped 
between  the  2  different  shEZFI2s  (Fig.  IB).  Further  data 
can  be  found  at  GEO  database  upon  publication  (GEO 
access  number:  GSE3 1433).  The  difference  in  the  number  of 
the  genes  altered  by  2  individual  shEZIT2s  may  be  due  to 
different  degrees  of  EZIT2  knockdown.  Consistent  with  this 
possibility,  shEZLI2#2,  which  decreases  EZH2  levels  with 
better  efficacy  than  shEZLI2#l,  resulted  in  a  greater  number 
of  upregulated  genes  (Fig.  1A  and  B).  Alternatively,  the 
differences  in  gene  upregulation  observed  with  individual 
EZLI2  shRNAs  may  be  due  to  off- target  effects.  To  avoid  this 
potential  issue,  we  chose  to  analyze  the  genes  that  are 
upregulated  by  both  shEZH2s.  Of  note,  some  of  the  known 
EZH2  target  genes  were  approaching  the  1.5-fold  upregula- 


A 
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Figure  1 .  Identification  of  EZH2  target  genes  in  human  SKOV3  EOC  cells. 
A,  SKOV3  cells  were  infected  with  indicated  lentivirus  encoding  shEZH2 
or  control.  Drug-selected  cells  were  examined  for  expression  of  EZH2, 
H3K27Me3,  and  GAPDH  by  immunoblotting  analysis  using  indicated 
antibodies.  B,  schematic  of  experimental  strategies  used  to  identify  EZH2 
target  genes.  Genes  whose  expression  was  upregulated  more  than  1 .5- 
fold  upon  EZH2  knockdown  by  2  individual  shEZH2  were  identified  by 
global  gene  expression  microarray  analysis.  Genomic  loci  occupied  by 
EZH2/H3K27Me3  were  profiled  by  ChIP-seq  analysis.  C,  cross- 
examination  of  gene  expression  profiling  and  ChIP-seq  analysis  as 
illustrated  in  (B)  revealed  a  list  of  60  putative  EZH2  target  genes  in  human 
SKOV3  EOC  cells.  GAPDH,  glyceraldehyde-3-phosphate 
dehydrogenase. 

tion  cutoff  point  but  were  not  included  in  further  analysis, 
including  VASH1  (5)  and  E-cadherin  (ref.  15;  data  not 
shown).  Although  the  conservative  approach  we  imple¬ 
mented  may  lead  to  missing  certain  EZH2  target  genes,  we 
felt  these  rigorous  methods  allowed  us  to  minimize  false¬ 
positive  EZFI2  target  genes  in  human  EOCs. 

We  next  sought  to  identify  genomic  loci  that  are  directly 
bound  by  EZLI2/H3K27Me3.  Toward  this  goal,  we  con¬ 
ducted  ChIP-seq  analysis  in  SKOV3  human  EOC  cells  using 
antibodies  specific  to  EZH2  or  H3K27Me3.  EZLI2  and 
I13K27Me3  occupancy  was  mapped  to  the  genomic  loci  of 
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3,223  and  5,922  genes,  respectively,  and  2,678  genes  were 
associated  with  both  EZH2  and  H3K27Me3  (Fig.  IB  and 
Supplementary  Table  SI).  The  difference  between  the  num¬ 
ber  of  genes  whose  locus  was  occupied  by  EZH2  and 
H3K27Me3  may  reflect  the  difference  in  affinity  of  the 
antibodies  used  for  ChlP.  Alternatively,  for  methylated  sites 
not  bound  by  EZH2,  it  is  possible  that  other  epigenetic 
regulators  in  addition  to  EZH2  can  also  generate 
H3K27Me3.  Consistently,  EZH1,  a  homolog  of  EZFI2  in 
human  cells,  is  also  capable  of  catalyzing  H3K27Me3  epi¬ 
genetic  modifications,  albeit  at  a  lower  rate  than  EZH2  (16). 
Furthermore,  genes  bound  by  EZH2,  but  not  H3K27Me3, 
may  reflect  an  H3K27Me3-independent  function  for  EZFI2 
as  previous  reports  have  suggested  (for  example,  see  ref.  17). 

To  identify  the  genes  that  are  directly  silenced  by  EZH2, 
we  cross-examined  the  gene  expression  and  ChIP-seq  data. 
As  a  result,  we  identified  a  list  of  60  EZH2/FI3K27Me3  target 
genes  whose  expression  was  upregulated  more  than  1.5-fold 
upon  EZH2  knockdown  in  SKOV3  human  EOC  cells  (Fig. 
1 C) .  Further  confirming  our  approach,  4  of  the  60  identified 
genes  have  previously  been  shown  as  EZH2/H3K27Me3 
target  genes,  namely,  SNAC  (18),  H19  (19),  DIRAS3  (20), 
and  DACT3  (21).  Notably,  Ingenuity  networks  analysis 
revealed  that  the  networks  enriched  by  the  identified  genes 
included  (i)  cell  death,  growth,  and  proliferation  (e.g.,  BMP, 
DAPK2,  NLRP1,  and  D1RAS3)  and  (ii)  reproductive  system 
development  and  cancer  (e.g.,  EAP2,  ALDH1A1 ,  SSTR1,  and 
MAGEH1;  data  not  shown).  This  is  consistent  with  the 
proliferation-promoting  and  apoptosis-suppressing  func¬ 
tion  of  EZFI2,  which  we  have  previously  reported  in  human 
EOC  cells  (4). 

Interestingly,  the  number  of  genes  upregulated  more  than 
1.5 -fold  upon  EZH2  knockdown  is  notably  lower  than  the 
number  of  genes  whose  genomic  loci  are  directly  occupied 
by  EZH2/H3K27Me3  (Fig.  IB).  This  result  suggests  that 
additional  mechanisms  may  cooperate  with  EZIT2/ 
H3K27Me3  in  silencing  or  reactivating  EZH2  target  genes. 
Consistent  with  this  idea,  previous  reports  have  shown  that 
EZH2  target  genes  are  also  subject  to  epigenetic  silencing  by 
H3K9Me3  (22)  or  histone  deacetylase  (23).  This  implies 
that  to  achieve  maximum  reactivation  of  EZIT2/ 
H3K27Me3-silenced  target  genes  in  human  EOC  cells, 
additional  epigenetic  gene  silencing  mechanisms  may  be 
considered  for  simultaneous  targeting  together  with  EZIT2 


inhibition.  Alternatively,  this  result  may  be  due  to  the 
bivalent  modification  (i.e.,  H3K27Me3  and  H3K4Me3)  at 
the  genomic  loci  of  those  upregulated  genes,  which  primes 
those  genes  for  activation  while  keeping  them  silenced  (24) . 
Further  studies  are  warranted  to  differentiate  these 
possibilities. 

Validation  of  the  selected  EZH2  target  genes  in  human 
EOC  cells 

The  list  of  upregulated  genes  was  prioritized  for  valida¬ 
tion  by  examining  their  expression  in  the  newly  released  the 
Cancer  Genomics  Atlas  (TCGA)  ovarian  database  (25).  We 
first  chose  those  genes  whose  expression  was  downregu- 
lated  more  than  2-fold  in  more  than  75%  of  EOC  cases  in 
TCGA  ovarian  database.  In  addition,  known  imprinted 
genes  such  as  H19  (19)  and  DIRAS3  (20)  or  poorly  anno¬ 
tated  genes  were  excluded.  Given  that  EZH2  promotes 
proliferation  and  invasion,  suppresses  apoptosis,  and  reg¬ 
ulates  stem  cell-like  population  in  human  EOCs  (4,  6),  we 
selected  3  identified  EZF12/H3K27Me3  target  genes  with 
one  or  more  of  these  roles  for  validation  studies.  Those 
genes  are  ALDH1A1  (11),  SSTR1  (26),  and  DACT3 
(ref.  21;  Table  1). 

We  first  validated  the  upregulation  of  the  selected  3  genes 
in  EZFI2-knockdown  SKOV3  human  EOC  cells  by  qRT- 
PCR.  Indeed,  all  3  selected  genes  were  significantly  upre¬ 
gulated  in  shEZH2-expressing  SKOV3  cells  compared  with 
controls  (Fig.  2,  P  <  0.05  vs.  controls.).  In  addition,  all  3 
genes  were  upregulated  by  both  shEZFI2s,  and  there  was  a 
correlation  between  the  degree  of  EZIT2  knockdown  and 
the  levels  of  upregulation  of  these  genes  (Figs .  1 A  and  2 ) .  We 
conclude  that  EZH2  knockdown  upregulates  the  expression 
of  ALDH1A1,  SSTR1,  and  DACT3  in  SKOV3  human  EOC 
cells. 

Next,  we  sought  to  validate  the  binding  of  EZFI2/ 
I!3K27Me3  to  the  genomic  loci  of  the  selected  genes. 
Toward  this  goal,  we  conducted  ChlP  analysis  using  anti¬ 
bodies  specific  to  EZH2  or  IT3K27Me3,  respectively.  An 
isotype-matched  IgG  was  used  as  a  negative  control,  and 
an  antibody  to  the  core  histone  H3  was  used  as  a  positive 
control  for  ChlP  analysis.  Indeed,  we  observed  the  binding 
of  both  EZH2  and  H3K27Me3  to  the  genomic  loci  of  the 
selected  EZH2  target  genes  in  SKOV3  human  EOC  cells  as 
determined  by  ChlP  analysis  (Fig.  3). 


Table  1.  Three  putative  EZH2  target  genes  identified  by  genome-wide  approaches  selected  for  further 
validation 


Gene  name 

Location 

Function 

%  TCGA  cases 
downregulated  >2-fold 

NM_000689  ALDH1A1 

9q21 .13 

Cancer  stem  cell  marker 

96 

NM„001049  SSTR1 

14q13 

Proliferation  and  invasion  inhibitor, 
cell  signaling 

95 

NM_1 45056  DACT3 

19q13.32 

Apoptosis  inducer,  Wnt  signaling 
antagonist 

76 
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Figure  2,  Validation  of  upregulation 
of  the  selected  EZH2-silenced 
target  genes  in  SKOV3  human  EOC 
cells  upon  EZH2  knockdown  by 
qRT-PCR.  SKOV3  cells  were 
infected  with  lentivirus  encoding 
the  indicated  shEZH2s  or  control. 
After  drug  selection,  mRNA  was 
extracted  and  examined  for 
expression  of  (A)  ALDH1A1  mRNA, 
(B)  SSTR1  mRNA,  and  (C)  DACT3 
mRNA  by  qRT-PCR.  Expression  of 
fi-2-microglobulin  was  used  to 
normalize  the  expression  of 
ALDH1A1,  SSTR1  and  DACT3 
mRNA.  *,  P  <  0.05  compared  with 
controls. 


We  next  sought  to  determine  whether  the  occupancy  of 
H3K27Me3  on  the  genomic  loci  of  EZH2  target  genes 
depends  upon  EZH2.  Toward  this  goal,  we  conducted  the 
ChIP  analysis  in  EZH2-knockdown  SKOV3  EOC  cells. 
Indeed,  knockdown  of  EZH2  severely  weakened  the  asso¬ 
ciation  of  both  EZH2  and  H3K27Me3  to  the  genomic  loci  of 
the  selected  EZH2  target  genes  (Fig.  3).  This  result  suggests 
that  EZH2  plays  a  major  role  in  regulating  H3K27Me3 
modification  on  the  genomic  loci  of  these  genes  in  human 
EOC  cells.  This  result  also  implies  that  the  binding  of  EZH2/ 
I  I3K27Me3  to  the  genomic  loci  of  these  genes  we  observed 
here  is  specific. 

Expression  of  EZH2  inversely  correlates  with 
expression  of  ALDH1A1 

We  next  sought  to  determine  whether  there  is  an  inverse 
correlation  between  expression  of  EZIT2  and  expression  of 
the  EZH2  target  genes  that  we  have  identified  and  validated 
in  this  study.  In  addition  to  EOC  cells,  EZH2  is  upregulated 


Figure  3.  Validation  of  occupancy  of  the  genomic  loci  of  the  selected  EZH2 
target  genes  by  EZH2  and  H3K27Me3  in  SKOV3  human  EOC  cells  using 
ChIP  analysis.  Control  and  shEZH2#2-expressing  SKOV3  cells  were 
subjected  to  ChIP  analysis  using  antibodies  specific  to  EZH2  or 
H3K27Me3,  respectively.  An  isotype-matched  IgG  was  used  as  a 
negative  control,  and  an  antibody  specific  to  core  histone  H3  was  used  as 
a  positive  control.  After  ChIP  analysis,  the  genomic  loci  of  the  indicated 
genes  were  subjected  to  PCR  amplification  using  primers  detailed  in 
Materials  and  Methods.  Please  see  Fig.  1 A  for  shEZH2  knockdown 
efficacy.  Shown  are  representative  images  of  3  independent 
experiments.  Note  that  for  H3K27Me3  ChIP,  a  low  number  of  PCR  cycles 
were  used  compared  with  EZH2  or  histone  H3  ChIP  to  avoid 
oversaturation  of  signals. 


in  ovarian  tumor- associated  endothelial  cells  (5).  To  limit 
the  confounding  effects  of  EOC-associated  stromal  cells 
(including  EOC-associated  endothelial  cells),  we  chose  to 
analyze  the  correlation  between  expression  of  EZH2  and  its 
target  genes  in  specimens  from  LCM  high-grade  serous 
subtype  EOC,  which  accounts  for  the  majority  of  EOC- 
associated  mortalities  (27). 

EZH2  is  expressed  at  higher  levels  in  human  EOC  cells 
than  in  normal  HOSE  cells  (4).  Therefore,  we  hypothesized 
that  EZH2  target  genes  that  are  silenced  by  EZH2  would  be 
expressed  at  lower  levels  in  human  EOC  cells.  Toward 
testing  this  hypothesis,  we  examined  the  expression  of 
EZH2  and  the  3  validated  EZH2  target  genes  in  a  published 
microarray  database,  which  compares  the  gene  expression 
profile  in  53  cases  of  LCM  high-grade  serous  EOCs  and  10 
individual  isolations  of  primary  HOSE  cells  (28).  Consis¬ 
tent  with  our  previous  report  (4),  EZH2  was  expressed  at 
significantly  higher  levels  in  human  EOCs  than  in  primary 
HOSE  cells  (P  <  0.001;  Fig.  4A  and  B).  Notably,  the  EZLI2 
target  gene  ALDH1A1  was  expressed  at  significantly  lower 
levels  in  human  EOCs  than  in  normal  HOSE  cells  (P  < 
0.001;  Fig.  4C  and  D).  Indeed,  there  was  a  negative  corre¬ 
lation  between  expression  of  EZH2  and  its  target  gene 
ALDH1A1  in  a  Spearman  statistical  analysis  of  the  cases 
including  both  EOC  and  primary  HOSE  cells  (P  <  0.00 1  and 
r  =  —0.41;  Fig.  4E,  including  both  open  circles  and  solid 
dots).  However,  the  coefficiency  Spearman  r  is  0.41.  This 
result  indicates  that  other  factors  may  also  play  a  role  in  the 
expression  relation.  Consistently,  there  is  evidence  to  sug¬ 
gest  that  Notch  signaling  also  regulates  ALDH1A1  expres¬ 
sion  (29).  In  addition,  the  correlation  between  expression 
of  EZH2  and  ALDI 1 1  Al  is  not  significant  among  EOC  cases 
(P  =  0.81;  Fig.  4E,  solid  dots  only).  This  may  be  due  to  the 
fact  the  ALDH 1  Al  is  expressed  at  very  low  levels  in  the  vast 
majority  of  EOC  cases,  and  thus,  the  variation  in  expression 
may  simply  be  a  reflection  of  experimental  variations. 

Comparing  EOCs  with  normal  HOSE  cells,  ALDH1A1 
showed  a  high  fold  change  in  expression  (>8-fold),  whereas 
SSTR1  (~3. 5-fold)  or  DACT3  (<1. 5-fold)  only  showed  a 
moderate  to  minimal  fold  change  in  expression  (Fig.  4D 
and  Supplementary  Fig.  SI).  It  is  possible  that  EZH2  is  the 
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2  4  6 

Relative  EZH2  expression 


Figure  4.  EZH2  targets  ALDH 1A1  in  high-grade  serous  subtype  EOC.  A,  relative  expression  of  EZH2  mRNA  in  53  cases  of  LCM  high-grade  serous  EOC  and  1 0 
individual  isolations  of  normal  HOSE  cells.  B,  quantitation  of  (A).  *,P<  0.001  compared  with  high-grade  serous  EOC.  C,  same  as  (A),  but  for  relative  expression 
of  ALDH1A1  mRNA.  D,  quantitation  of  (C).  *,  P<  0.001  compared  with  high-grade  serous  EOC.  E,  correlation  between  expression  of  ALDH1A1  ar\6EZH2  as 
determined  by  Spearman  statistical  analysis  using  GraphPad  Prism  version  5.0  software. 


major  epigenetic  regulator  of  ALDH1A1,  whereas  addition¬ 
al  epigenetic  silencing  mechanisms  may  contribute  to  sup¬ 
pression  of  the  other  2  validated  EZH2  target  genes.  In 
support  of  this  possibility,  EZH2  knockdown  induced 
much  greater  levels  of  upregulation  for  ALDH1A1  (up  to 
22-fold)  than  the  other  2  EZH2  target  genes  (Fig.  2). 

Next,  we  sought  to  examine  the  correlation  between 
expression  of  EZH2  or  ALDH1A1  and  survival  of  patients 
with  high-grade  serous  EOC.  Consistent  with  our  previous 
report  (4),  the  difference  in  overall  survival  between  high 
and  low  EZH2  expression  in  patients  with  high-grade  serous 
EOC  was  not  significant  (P  =  0.1684;  Supplementary  Fig. 
S2).  In  addition,  the  difference  in  survival  between  low 
ALDH1A1  expression  group  and  highALDEIlAl  group  was 
not  significant  (P  =  0.7789;  Supplementary  Fig.  S2).  In 
contrast  to  our  results,  high  ALDH  1 A  has  previously  been 
reported  to  be  associated  with  poor  prognosis  in  patients 
with  EOC  ( 1 0,  30) .  The  basis  for  this  discrepancy  remains  to 
be  determined.  However,  the  discrepancy  could  be  due  to 
different  methods  that  were  used  in  this  study  (microarray) 
versus  the  other  2  studies  (immunohistochemical  staining) . 

ALDH1A1  has  been  reported  as  a  marker  of  cancer  stem 
cells  in  certain  types  of  cancers  including  breast  and  ovarian 
(10,  29,  31-33).  Likewise,  EZH2  plays  an  important  role  in 
stem  cell  biology  (34).  Our  data  suggest  that  EZH2  directly 
regulates  the  levels  ofALDHIAl  in  EOC  cells,  implying  that 
EZH2  may  regulate  the  EOC  stem  cell  population  by  con¬ 
trolling  the  levels  of  ALDH1A1  expression.  Similarly,  it  has 
been  shown  that  EZEI2  directly  regulates  the  epigenetic 
status  of  Nanog,  an  important  factor  for  both  embryonic 


stem  cell  and  induced  pluripotent  stem  cells,  to  balance  the 
equilibrium  between  self-renewal  and  differentiation  of 
stem  cells  (35).  Similarly,  other  putative  markers  of  cancer 
stem  cells  [such  as  CD133  (ref.  36)  andTACSTD2  (ref.  37)] 
have  been  also  reported  to  be  hypermethylated  in  cancerous 
cells.  While  the  expression  in  stem  cells  will  be  masked  by 
the  vastly  more  abundant  non-stem  cell  population  in  our 
or  any  similar  analysis,  it  is  nevertheless  intriguing  that 
differential  expression  ofALDHIAl  has  been  described  as  a 
marker  of  sternness  in  cancer,  including  in  EOC  (10,  33). 
Together,  our  data  suggest  that  EZIT2  may  regulate  the  EOC 
stem  cell  population  by  controlling  the  levels  of  ALDI 1 1A1 
expression.  Our  future  work  will  test  this  hypothesis. 

In  the  present  study,  using  a  combination  of  global 
gene  expression  profiling  and  genome-wide  ChlP-seq 
analysis,  we  identified  a  list  of  60  EZ1I2  direct  target 
genes,  whose  expression  was  upregulated  more  than 
1.5-fold  upon  EZEI2  knockdown  in  human  EOC  cells. 
These  genes  include  56  novel  putative  EZEI2  target  genes 
and  4  known  EZH2  target  genes.  We  validated  3  selected 
EZEI2  target  genes  that  are  implicated  in  regulating  cancer 
stem  cells,  cell  proliferation,  apoptosis,  and  cell  invasion. 
We  showed  that  ALDH1A1,  a  putative  marker  of  EOC 
stem  cells  (11),  was  expressed  at  lower  levels  in  high- 
grade  serous  EOCs  than  in  normal  HOSE  cells,  and  there 
was  a  negative  correlation  between  expression  of  EZH2 
and  expression  of  ALDII1A1.  Further  studies  are  war¬ 
ranted  to  mine  the  data  presented  here  as  well  as  func¬ 
tional  characterization  of  the  identified  EZH2  target 
genes.  These  studies  should  provide  important  insights 
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into  the  biology  of  EOC  development  and  the  identifi¬ 
cation  of  potential  candidate  targets  for  prevention  and 
intervention  of  EOC. 
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Abstract 

Epithelial  ovarian  cancer  (EOC)  remains  the  most  lethal  gynecologic  malignancy  in  the  United  States.  Thus, 
there  is  an  urgent  need  to  develop  novel  therapeutics  for  this  disease.  Cellular  senescence  is  an  important  tumor 
suppression  mechanism  that  has  recently  been  suggested  as  a  novel  mechanism  to  target  for  developing  cancer 
therapeutics.  Wnt5a  is  a  noncanonical  Wnt  ligand  that  plays  a  context-dependent  role  in  human  cancers.  Here,  we 
investigate  the  role  of  Wnt5a  in  regulating  senescence  of  EOC  cells.  We  show  that  Wnt5a  is  expressed  at 
significantly  lower  levels  in  human  EOC  cell  lines  and  in  primary  human  EOCs  (n  =  130)  compared  with  either 
normal  ovarian  surface  epithelium  (n  =  31-,P  =  0.039)  or  fallopian  tube  epithelium  (n  =  28 \P<  0.001).  Notably,  a 
lower  level  of  Wnt5a  expression  correlates  with  tumor  stage  (P  =  0.003)  and  predicts  shorter  overall  survival  in 
EOC  patients  ( P=  0.003).  Significantly,  restoration  of  Wnt5a  expression  inhibits  the  proliferation  of  human  EOC 
cells  both  in  vitro  and  in  vivo  in  an  orthotopic  EOC  mouse  model.  Mechanistically,  Wnt5a  antagonizes  canonical 
Wnt/|3-catenin  signaling  and  induces  cellular  senescence  by  activating  the  histone  repressor  A/promyelocytic 
leukemia  senescence  pathway.  In  summary,  we  show  that  loss  of  Wnt5a  predicts  poor  outcome  in  EOC  patients 
and  Wnt5a  suppresses  the  growth  of  EOC  cells  by  triggering  cellular  senescence.  We  suggest  that  strategies  to 
drive  senescence  in  EOC  cells  by  reconstituting  Wnt5a  signaling  may  offer  an  effective  new  strategy  for  EOC 
therapy.  Cancer  Res;  71(19);  6184-94.  ©2011  AACR. 


Introduction 

Cellular  senescence  is  an  important  tumor  suppression 
mechanism  in  vivo  (1).  In  primary  mammalian  cells,  cellular 
senescence  can  be  triggered  by  various  inducers  including 
critically  shortened  telomeres  and  activated  oncogenes  (such 
as  oncogenic  RAS-,  ref.  1).  Senescent  cells  are  viable  but  non¬ 
dividing  (2).  Senescent  cells  also  exhibit  several  distinctive 
morphologic  characteristics  and  molecular  markers,  including 
a  large  flat  cellular  morphology  and  expression  of  senescence- 
associated  p-galactosidase  (SA-|3-gal)  activity  (3).  In  murine 
liver  carcinoma  and  sarcoma  models,  reactivation  of  the  tumor 
suppressor  p53  induces  senescence  and  is  associated  with 
tumor  regression  (4,  5).  Hence,  driving  cancer  cells  to  undergo 
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cellular  senescence  represents  a  novel  mechanism  for  devel¬ 
oping  cancer  therapeutics  (6,  7). 

More  than  85%  of  ovarian  cancers  are  of  epithelial  origin  (8). 
Epithelial  ovarian  cancers  (EOC)  are  classified  into  distinct 
histologic  types  including  serous,  mucinous,  endometrioid, 
and  clear  cell  (9).  The  most  common  histology  of  EOC  is  serous 
(~60%  of  all  cancers)  and  less  common  histologies  include 
endometrioid,  clear  cell,  and  mucinous  (9).  Recently,  an  alter¬ 
native  classification  has  been  proposed,  in  which  EOC  is 
broadly  divided  into  2  types  (10).  Type  I  EOC  includes  endo¬ 
metrioid,  mucinous,  low-grade  serous,  and  clear-cell  carcino¬ 
mas,  and  type  II  EOC  includes  high-grade  serous  carcinomas 
(10).  EOC  remains  the  most  lethal  gynecologic  malignancy  in 
the  United  States  (8).  Thus,  there  is  an  urgent  need  to  better 
understand  the  etiology  of  EOC  to  develop  novel  therapeutics 
for  this  devastating  disease. 

Wnt  signaling  is  initiated  by  binding  of  the  Wnt  ligand  to  its 
cognate  frizzled  receptor  (11).  Canonical  Wnt  signaling  results 
in  stabilization  of  the  key  transcription  factor  [l-catenin,  which 
then  translocates  into  the  nucleus  and  drives  expression  of  its 
target  genes,  such  as  CCND1  (cyclin  Dl),  FOSL1,  and  c-MYC 
(12,  13).  Canonical  Wnt  signaling  is  active  in  the  putative 
somatic  stem/progenitor  cells  of  the  coelomic  epithelium  of 
the  mouse  ovary  (14).  Underscoring  the  importance  of  Wnt 
signaling  in  EOC,  in  a  murine  ovarian  cancer  model,  activation 
of  canonical  Wnt  signaling  cooperates  with  inactivation  of  the 
tumor  suppressor  PTEN  in  driving  ovarian  carcinogenesis  (15). 
However,  the  role  of  Wnt  signaling  in  EOC  is  not  fully 
understood. 
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Wnt5a  is  a  noncanonical  Wnt  ligand  that  plays  opposing 
roles  in  different  types  of  cancer  and  has  variable  expression 
dependent  on  the  cancer  context  (16).  Specifically,  in  EOC  the 
role  of  Wnt5a  remains  unclear.  Thus,  in  this  study,  we  inves¬ 
tigated  Wnt5a  expression  and  its  potential  function  in  human 
EOC  cells.  We  discovered  that  Wnt5a  was  expressed  at  signif¬ 
icantly  lower  levels  in  primary  human  EOC  compared  with 
either  primary  human  ovarian  surface  epithelium  or  fallopian 
tube  epithelium.  Notably,  loss  of  Wnt5a  expression  was  asso¬ 
ciated  with  tumor  stage  and  predicted  shorter  overall  survival 
in  EOC  patients.  Significantly,  Wnt5a  reconstitution  inhibited 
the  growth  of  EOC  cells  both  in  vitro  and  in  vivo  in  an  orthotopic 
EOC  mouse  model  by  promoting  cellular  senescence.  These 
studies  show,  for  the  first  time,  a  functional  role  of  the 
noncanonical  Wnt  ligand,  Wnt5a,  in  promoting  senescence. 
Importantly,  they  also  suggest  that  promoting  EOC  cells  to 
undergo  senescence  represents  a  potential  novel  strategy  for 
developing  urgently  needed  EOC  therapeutics. 

Materials  and  Methods 

Cells  and  culture  conditions 

Primary  human  ovarian  surface  epithelial  (HOSE)  cells  were 
isolated  and  cultured  as  previously  described  (17).  Human  EOC 
cell  lines  were  obtained  from  American  Type  Culture  Collec¬ 
tion  (ATCC)  and  were  passaged  for  less  than  6  months.  EOC 
cell  line  identification  was  further  confirmed  by  DNA  Diag¬ 
nostic  Center  (www.dnacenter.com).  EOC  cell  fines  were  cul¬ 
tured  according  to  ATCC  in  RPMI-1640  medium  supplemented 
with  10%  FBS.  5-Aza-cytidine  (Aza-C;  Sigma)  was  used  at 
working  concentration  of  5  Umol/L  (18). 

Human  ovarian  specimens  and  immunohistochemistry 

The  protocol  to  evaluate  deidentified  human  tissue  speci¬ 
mens  was  approved  by  Fox  Chase  Cancer  Center  (FCCC) 
institutional  review  board.  Ovarian  tumor  microarray  and 
normal  human  ovary  and  fallopian  tube  specimens  were 
obtained  from  the  FCCC  Biosample  Repository  Core  Facility 
(BRCF).  Histopathology  of  the  selected  specimens  on  the 
tumor  microarrays  was  provided  by  BRCF.  Immunohis¬ 
tochemistry  (IHC)  was  conducted  by  using  goat  anti-Wnt5a 
polyclonal  antibody  (R&D  Systems)  and  mouse  anti-Ki-67 
(Dako)  with  a  DAKO  EnVision  System  and  the  Peroxidase 
(DAB)  kit  (DAKO  Corporation)  following  the  manufacturer's 
instructions  and  as  previously  described  (17).  Wnt5a  staining 
was  scored  in  a  double-blinded  manner  by  Dr.  Qi  Cai  at  the 
BRCF,  and  a  proportion  of  the  cases  were  independently 
confirmed  by  Dr.  Hong  Wu,  a  board-certified  pathologist,  at 
the  FCCC  Department  of  Pathology. 

Anchorage-independent  soft  agar  colony  formation 
assay 

Soft  agar  assay  were  carried  out  as  previously  described  (17). 
Briefly,  3,500  cells  were  resuspended  in  0.35%  low  melt  agarose 
dissolved  in  RPMI-1640  medium  supplemented  with  10%  FBS 
and  inoculated  on  top  of  0.6%  low  melt  agarose  base  in  6-well 
plates.  After  2  weeks  in  culture,  the  plates  were  stained  with 


0.005%  crystal  violet,  and  the  number  of  colonies  was  counted 
by  using  a  dissecting  microscope. 

Retrovirus  production,  infection,  and  drug  selection 

The  following  retrovirus  constructs  were  used:  pBABE-puro 
was  obtained  from  Addgene,  hygro-pWZL-luciferase  was  a 
kind  gift  of  Dr.  Denise  Connolly,  and  pBABE-Wnt5a  was 
generated  by  using  standard  cloning  protocol.  Retrovirus 
packaging  was  done  as  previously  described  by  using  Phoenix 
packaging  cells  (19,  20).  To  increase  infection  efficacy,  double 
virus  infection  was  carried  out.  For  drug  selection,  3  |tg/mL  of 
puromycin  was  used  for  the  OVCAR5  human  EOC  cell  line. 

Reverse  transcriptase  PCR,  quantitative  reverse 
transcriptase  PCR,  and  immunoblotting 

RNA  from  cultured  primary  HOSE  cells  or  human  EOC  cell 
lines  was  isolated  by  using  TRIzol  (Invitrogen)  according  to 
manufacturer's  instruction.  For  quantitative  reverse  transcrip¬ 
tase  PCR  (qRT-PCR),  TRIzol-isolated  RNA  was  further  purified 
by  using  an  RNeasy  kit  (QIAGEN)  following  manufacture's 
instruction.  The  Wnt5a,  CCND1,  FOSL1,  and  c-MYC  primers 
used  for  qRT-PCR  were  purchased  from  SABiosciences.  mRNA 
expression  of  the  housekeeping  gene  P-2-microglobulin  (B2M) 
was  used  to  normalize  mRNA  expression.  Soluble  p-calenin 
was  extracted,  using  a  buffer  that  consisting  of  10  mmoI/L  Tris- 
HC1  (pH  7.5),  0.05%  NP-40, 10  mmol/L  NaCl,  3  mmol/L  MgCk,  1 
mmoI/L  EDTA,  and  protease  inhibitors  (Roche)  as  previously 
described  (21,  22).  The  following  antibodies  were  used  for 
immunoblotting  from  the  indicated  suppliers,  goat  anti-Wnt5a 
(R&D  Systems),  mouse  p-actin  (Sigma),  mouse  anti-GAPDH 
(Millipore),  mouse  anti-p-catenin,  mouse  anti-Rb  (BD  Bios¬ 
ciences),  and  rabbit  anti-pRBpS780  (Cell  Signaling). 

Immunofluorescence  and  SA-fl-gal  staining 

Indirect  immunofluorescence  staining  was  carried  out  as 
previously  described  (19, 20, 22).  The  following  antibodies  were 
used  for  immunofluorescence:  a  cocktail  of  mouse  anti-HIRA 
monoclonal  antibodies  (WC19,  WC117,  and  WC119;  1:10; 
ref.  20)  and  a  rabbit  anti-PML  antibody  (Chemicon,  1:5,000). 
Images  were  captured  by  a  DS-QilMc  camera  on  a  Nikon 
Eclipse  80i  microscope  and  processed  by  NIS-Elements  BR3.0 
software  (Nikon).  SA-P-gal  staining  was  carried  out  as 
described  previously  (3,  23).  For  SA-P-gal  staining  in  sections 
from  xenografted  tumors,  8  separate  fields  were  examined 
from  2  individual  tumors  for  each  of  the  groups. 

In  vivo  orthotopic  xenograft  tumorigenesis  study 

The  protocol  was  approved  by  the  FCCC  Institutional  Ani¬ 
mal  Care  and  Use  Committee.  OVCAR5  cells  were  infected  with 
a  luciferase-encoding  retrovirus  (hygro-pWZL-luciferase)  and 
infected  cells  were  selected  with  50  lig/ml,  hygromycin.  Drug- 
selected  cells  were  then  infected  with  control  or  Wnt5a-encod- 
ing  retrovirus  and  subsequently  selected  with  3  lig/ml,  puro¬ 
mycin  and  50  pg/mL  hygromycin.  A  total  of  3  x  106  drug- 
selected  cells  were  unilaterally  injected  into  the  ovarian  bursa 
sac  of  immunocompromised  mice  (6  mice  per  group;  ref.  24). 
From  day  10  postinfection,  tumors  were  visualized  by  injecting 
luciferin  (intraperitoneal,  4  mg/mice)  resuspended  in  PBS  and 
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imaged  with  an  IVIS  Spectrum  imaging  system  every  5  days 
until  day  30.  Images  were  analyzed  by  Live  Imaging  4.0  software. 
At  day  30,  tumors  were  surgically  dissected  and  either  fixed  in 
10%  formalin  or  fresh-frozen  in  Optimal  Cutting  Temperature 
compound  (Tissue-Tek).  Sections  of  the  dissected  tumors  were 
processed  by  the  FCCC  Histopathology  Core  Facility. 

Statistical  analysis 

Quantitative  data  are  expressed  as  mean  ±  SD,  unless 
otherwise  indicated.  ANOVA  with  Student's  t  test  was  used 
to  identify  significant  differences  in  multiple  comparisons.  The 
Pearson  test  was  used  to  analyze  the  relationship  between 
categorical  variables.  Overall  survival  was  defined  as  the  time 
elapsed  from  the  date  of  diagnosis  and  the  date  of  death  from 
any  cause  or  the  date  of  last  follow-up.  Kaplan-Meier  survival 
plots  were  generated  and  comparisons  were  made  by  using  the 
log-rank  sum  statistic.  For  all  statistical  analyses,  the  level  of 
significance  was  set  at  0.05. 

Results 

Wnt5a  is  expressed  at  significantly  lower  levels  in  human 
EOC  cell  lines  and  primary  human  EOCs  compared  with 
normal  human  ovarian  surface  epithelium  or  fallopian 
tube  epithelium 

To  determine  Wnt5a  expression  in  human  EOC  cell  lines  and 
primary  HOSE  cells,  we  examined  the  relative  Wnt5a  mRNA 
levels  by  carrying  out  semiquantitative  RT-PCR.  We  observed 
that  Wnt5a  mRNA  levels  were  greatly  diminished  in  human 
EOC  cell  lines  compared  with  primary  HOSE  cells  (Fig.  1A). 
This  finding  was  further  confirmed  through  qRT-PCR  analysis 
of  Wnt5a  mRNA  in  multiple  isolations  of  primary  HOSE  cells 
and  human  EOC  cell  lines,  showing  that  the  levels  of  Wnt5a 
mRNA  were  significantly  lower  in  human  EOC  cell  lines 
compared  with  primary  HOSE  cells  (Fig.  IB;  P  =  0.008). 
Consistently,  we  observed  that  Wnt5a  protein  levels  were  also 
lower  in  human  EOC  cell  lines  compared  with  primary  HOSE 
cells  as  determined  by  immunoblotting  (Fig.  1C).  On  the  basis 
of  these  results,  we  conclude  that  Wnt5a  is  expressed  at  lower 
levels  in  human  EOC  cell  lines  compared  with  primary  HOSE 
cells. 

We  next  determined  whether  the  loss  of  Wnt5a  expression 
found  in  human  EOC  cell  lines  was  also  observed  in  primary 
human  EOCs.  We  examined  Wnt5a  expression  in  130  cases  of 
primary  human  EOC  specimens  and  31  cases  of  normal  human 
ovary  with  surface  epithelium  by  IHC,  using  an  antibody 
against  Wnt5a  (Table  1).  In  addition,  there  is  recent  evidence 
to  suggest  that  a  proportion  of  high-grade  serous  EOC  may 
arise  from  distant  fallopian  tube  epithelium  (25).  Thus,  we  also 
included  28  cases  of  normal  human  fallopian  tube  specimens  in 
our  IHC  analysis  (Table  1).  The  specificity  of  the  anti-Wnt5a 
antibody  was  confirmed  in  our  study  (Supplementary  Fig.  SI). 
A  single  band  at  predicted  molecular  weight  (~42  kDa)  was 
detected  in  OVCAR5  cells  with  ectopically  expressed  Wnt5a 
and  was  absent  after  expression  of  a  short  hairpin  RNA  to  the 
human  Wnt5a  gene  (shWnt5a),  which  effectively  knocked 
down  Wnt5a  mRNA  expression  (Supplementary  Fig.  S1A  and 
data  not  shown).  In  addition,  Wnt5a  staining  was  lost  when 


primary  anti-Wnt5a  antibody  was  replaced  with  an  isotype- 
matched  IgG  control  (Supplementary  Fig.  SIB). 

As  shown  in  Figure  ID,  in  normal  human  ovarian  surface 
epithelial  cells  and  fallopian  tube  epithelial  cells,  both  cyto¬ 
plasm  and  cell  membrane  were  positive  for  Wnt5a  IHC  staining 
(black  arrows,  Fig.  ID).  In  contrast,  Wnt5a  staining  in  EOC  cells 
was  dramatically  decreased  (Fig.  ID).  We  scored  expression  of 
Wnt5a  as  high  (//-score  >  30)  or  low  (//  score  <  30)  on  the  basis 
of  a  histological  score  ( H  score;  26),  which  considers  both 
intensity  of  staining  and  percentage  of  positively  stained  cells, 
as  previously  described  (17).  Wnt5a  expression  was  scored  as 
high  in  58.1%  (18/31)  cases  of  normal  human  ovarian  surface 
epithelium  and  82.1%  (23/28)  cases  of  normal  human  fallopian 
tube  epithelium  (Table  1).  In  contrast,  Wnt5a  expression  was 
scored  as  high  in  37.7%  (49/ 130)  cases  of  primary  human  EOCs 
(Table  1).  Statistical  analysis  revealed  that  Wnt5a  was 
expressed  at  significantly  lower  levels  in  primary  human  EOCs 
compared  with  either  normal  human  ovarian  surface  epithe¬ 
lium  (P=  0.039)  or  normal  human  fallopian  tube  epithelium  (P 
<  0.001;  Table  1).  On  the  basis  of  these  studies,  we  conclude  that 
Wnt5a  is  expressed  at  significantly  lower  levels  in  primary 
human  EOCs  compared  with  either  normal  human  ovarian 
surface  epithelium  or  fallopian  tube  epithelium. 

Wnt5a  expression  negatively  correlates  with  tumor  stage 
and  lower  VVnt5a  expression  predicts  shorter  overall 
survival 

We  next  examined  the  correlation  between  Wnt5a  expres¬ 
sion  and  clinical  and  pathologic  features  of  human  EOCs. 
Significantly,  there  was  a  negative  correlation  between  Wnt5a 
expression  and  tumor  stage  (P  =  0.003;  Table  1).  Notably,  the 
majority  of  examined  cases  are  high-grade  serous  subtypes 
that  are  usually  of  stage  3/4.  In  addition,  we  examined  the 
correlation  between  expression  of  Wnt5a  and  a  marker  of  cell 
proliferation,  Ki-67  (ref.  27;  Fig.  IE).  There  was  a  significant 
negative  correlation  between  Wnt5a  expression  and  Ki-67  (P  = 
0.038;  Table  1).  We  next  assessed  whether  Wnt5a  expression 
based  on  H  score  might  predict  prognosis  of  EOC  patients 
(High,  //score  >  30;  Low,  //score  <30  ;n  =  123),  for  which  long¬ 
term  follow-up  data  were  available.  Significantly,  lower  Wnt5a 
expression  correlated  with  shorter  overall  survival  in  the 
examined  EOC  patients  ( P  =  0.003;  Fig.  IF).  Together,  we 
conclude  that  a  lower  level  of  Wnt5a  expression  correlates 
with  tumor  stage  and  predicts  shorter  overall  survival  in 
human  EOC  patients. 

Wnt5a  gene  promoter  hypermethylation  contributes  to 
its  downregulation  in  human  EOC  cells 

Wnt5a  gene  promoter  hypermethylation  has  been  implicat¬ 
ed  as  a  mechanism  underlying  its  silencing  in  several  types  of 
human  cancers  (16).  Consistently,  we  also  observed  Wnt5a 
gene  promoter  hypermethylation  in  a  number  of  human  EOC 
cell  lines  (Fig.  2A;  Supplementary  Table  SI).  Further  supporting 
a  role  of  promoter  hypermethylation  in  suppression  of  Wnt5a 
expression,  treatment  with  a  DNA  demethylation  drug,  Aza-C 
(28),  in  PEOl  EOC  cells  resulted  in  a  significant  increase  in 
levels  of  both  Wnt5a  mRNA  and  protein  (Fig.  2B  and  C).  We 
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Figure  1.  Wnt5a  is  expressed  at  significantly  lower  levels  in  human  EOC  cells  compared  with  normal  human  ovarian  surface  or  fallopian  tube  epithelial 
cells,  and  a  lower  level  of  Wnt5a  expression  predicts  shorter  overall  survival  in  human  EOC  patients.  A,  expression  of  Wnt5a  mRNA  in  primary  HOSE  cells  and 
the  indicated  human  EOC  cell  lines  was  determined  by  semiquantitative  RT-PCR.  Expression  of  glyceraldehyde  3-phosphate  dehydrogenase 
(GAPDH)  mRNA  was  used  as  a  loading  control.  B,  Wnt5a  mRNA  levels  were  quantified  by  qRT-PCR  in  6  individual  isolations  of  primary  HOSE  cells  and  7 
different  EOC  cell  lines.  Expression  of  p-2-microglobulin  was  used  to  normalize  Wnt5a  mRNA  expression.  * ,  P  =  0.008  compared  with  human  EOC  cells.  C, 
same  as  (A),  but  examined  for  Wnt5a  and  GAPDH  protein  expression  by  immunoblotting.  D,  examples  of  Wnt5a  IHC  staining  in  normal  human  ovarian  surface 
epithelium,  fallopian  tube  epithelium,  and  EOC  of  indicated  histologic  subtypes.  Bar,  50  pm.  Arrows  point  to  examples  of  positively  stained  human  ovarian 
surface  epithelial  cells  and  fallopian  tube  epithelial  cells.  E,  representative  images  from  tissue  microarray  depicting  low  Wnt5a  expression  correlated 
with  high  Ki-67,  a  cell  proliferation  marker.  F,  loss  of  Wnt5a  expression  is  an  independent  poor  prognosis  marker  in  human  EOC  patients.  A  lower  level  of  Wnt5a 
expression  correlates  with  shorter  overall  survival  in  human  EOC  patients.  The  univariate  overall  survival  curve  (Kaplan-Meier  method)  for  EOC  patients 
(n  =  123)  with  high-  or  low-Wnt5a  expression  as  determined  by  immunohistochemical  analysis. 


conclude  that  WntSa  gene  promoter  hypermethylation  con¬ 
tributes  to  its  downregulation  in  human  EOC  cells. 

Wnt5a  restoration  inhibits  the  growth  of  human  EOC 
cells  by  antagonizing  the  canonical  Wilt/ p-catenin 
signaling 

We  next  sought  to  determine  the  effects  of  Wnt5a  recon¬ 
stitution  in  human  EOC  cells.  Wnt5a  expression  was  recon¬ 
stituted  in  the  OVCAR5  EOC  cell  line  via  retroviral  trans¬ 
duction.  Ectopically  expressed  Wnt5a  was  confirmed  by 


both  qRT-PCR  and  immunoblotting  in  OVCAR5  cells  stably 
expressing  Wnt5a  or  a  vector  control  (Fig.  3A  and  B).  Of  note, 
the  levels  of  ectopically  expressed  Wnt5a  in  OVCAR5  cells 
are  comparable  with  the  levels  observed  in  primary  HOSE 
cells  (Fig.  3B).  Interestingly,  Wnt5a  reconstitution  in 
OVCAR5  human  EOC  cells  significantly  inhibited  both 
anchorage-dependent  and  anchorage-independent  growth 
in  soft  agar  compared  with  vector  controls  (Fig.  3C  and 
D).  In  addition,  similar  growth  inhibition  by  Wnt5a  recon¬ 
stitution  was  also  observed  in  the  PEOl  human  EOC  cell  line 
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Table  1.  Wnt5a  expression  in  primary  human  EOCsand  correlation  of  its  expression  with  clinicopathologic 


variables 

Patient  characteristics 

Wnt5a  protein  expression 

Low  (n) 

High  (n) 

Total  (n) 

High  (%) 

P 

Age  ( 23-85  y;  mean  59.2  y) 

<55 

24 

16 

40 

40.0 

>55 

52 

33 

85 

38.8 

0.900 

Unknown 

5 

0 

5 

Laterality 

Left 

22 

14 

36 

38.9 

Right 

12 

9 

21 

42.9 

0.957 

Bilaterality 

35 

24 

59 

40.7 

Undetermined 

12 

2 

14 

Histotype 

EOC 

81 

49 

130 

37.7 

Type  1 

16 

21 

37 

56.8 

Low-grade  serous 

1 

1 

2 

50.0 

Endometrioid 

4 

9 

13 

69.2 

Mucinous 

2 

3 

5 

60.0 

Clear  cell 

5 

4 

9 

44.4 

Others 

4 

4 

8 

50.0 

Type  II 

High-grade  serous 

65 

28 

93 

30.1 

0.005a 

Normal  epithelium 

Ovarian  surface 

13 

18 

31 

58.1 

0.039b 

Fallopian  tube 

5 

23 

28 

82.1 

<0.001 b 

Ki-67 

Low 

22 

23 

44 

52.3 

High 

51 

24 

75 

32.0 

0.038 

Undetermined 

7 

3 

11 

Tumor  grade 

1 

3 

7 

10 

70.0 

2 

12 

8 

20 

40.0 

3 

64 

31 

95 

32.6 

Undetermined 

2 

3 

5 

Tumor  stage 

Stage  1/2 

12 

18 

30 

60.0 

Stage  3/4 

67 

29 

96 

30.2 

0.003c 

Undetermined 

2 

2 

4 

aCompared  with  type  I  EOC. 
bCompared  with  EOC. 
cCompared  with  stage  1/2. 


(Supplementary  Fig.  S2A-C)  suggesting  that  this  effect  is  not 
cell  line  specific.  On  the  basis  of  these  results,  we  conclude 
that  Wnt5a  reconstitution  inhibits  the  growth  of  human 
EOC  cells  in  vitro. 

Canonical  Wnt  signaling  promotes  cell  proliferation  and 
Wnt5a  has  been  shown  to  antagonize  the  canonical  Wnt/ 
p-catenin  signaling  in  certain  cell  contexts  (16, 29-31).  Because 
Wnt5a  expression  inversely  correlated  with  expression  of  Ki-67 
(Fig.  IE;  Table  1),  a  cell  proliferation  marker,  we  hypothesized 
that  Wnt5a  would  suppress  the  growth  of  human  EOC  cells  by 


antagonizing  canonical  Wnt/p-catenin  signaling.  To  test  our 
hypothesis,  we  examined  the  effect  of  Wnt5a  reconstitution  on 
expression  of  markers  of  active  Wnt/p-catenin  signaling 
in  human  EOC  cells,  namely  the  levels  of  "active"  soluble 
p-catenin  (21,  22,  32)  and  expression  of  p-catenin  target  genes 
such  as  CCND1,  c-MYC,  and  FOSL1  (12, 13).  Indeed,  we  observed 
a  decrease  in  soluble  p-catenin  in  Wnt5a-reconstituted 
OVCAR5  cells  compared  with  vector  controls  (Fig.  3E).  Con¬ 
sistently,  we  also  observed  a  significant  decrease  in  the  levels  of 
p-catenin  target  genes  in  these  cells,  namely  CCND1 
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Figure  2.  Promoter  DNA  CpG  island  hypermethylation  contributes  to 
Wnt5a  downregulation  in  human  EOC  cells.  A,  schematic  structure  of  the 
human  Wnt5a  gene  transcript  and  its  promoter  CpG  islands.  Locations  of 
exon  1  (open  rectangle),  CpG  sites  (vertical  lines)  and  coding  exons  (filled 
rectangle),  and  the  transcription  start  site  (curved  arrow)  are  indicated. 
Flat  arrows  indicate  the  positions  of  primers  used  for  PCR  amplification, 
and  the  size  of  PCR  product  is  also  indicated.  B,  PEOI  cells  were  treated 
with  5  |imol/L  Aza-C  for  4  days,  and  mRNA  was  isolated  from  control-  and 
Aza-C-treated  cells  and  examined  for  Wnt5a  mRNA  expression  by  qRT- 
PCR.  Mean  of  3  independent  experiments  with  SD.  C,  same  as  (B)  but 
examined  for  Wnt5a  protein  expression  by  immunoblotting. 

(. P  =  0.0095),  FOSL1  (P  =  0.0012),  and  c-MYC  (P  =  0.0286;  Fig. 
3F).  Similar  effects  of  Wnt5a  reconstitution  on  expression  of 
markers  of  active  Wnt/p-catenin  signaling  (such  as  decreased 
levels  of  soluble  (3-catenin)  were  also  observed  in  PEOI  human 
EOC  cells  (Supplementary  Fig.  S2D),  suggesting  that  this  is  not 
cell  line  specific.  On  the  basis  of  these  results,  we  conclude  that 
Wnt5a  suppresses  the  growth  of  human  EOC  cells  by  antag¬ 
onizing  canonical  Wnt/p-catenin  signaling  in  human  EOC 
cells. 

Wnt5a  reconstitution  drives  cellular  senescence  in 
human  EOC  cells 

Next,  we  sought  to  determine  the  cellular  mechanism 
whereby  Wnt5a  inhibits  the  growth  of  human  EOC  cells.  We 
have  previously  shown  that  suppression  of  canonical  Wnt 
signaling  promotes  cellular  senescence  in  primary  human 
fibroblasts  by  activating  the  senescence-promoting  histone 
repressor  A  (HIRA)/promyelocytic  leukemia  (PML)  pathway 
(22).  PML  bodies  are  20  to  30  dot-like  structures  in  the  nucleus 
of  virtually  all  human  cells.  PML  bodies  are  sites  of  poorly 
defined  tumor  suppressor  activity  and  are  disrupted  in  acute 
PML  (33).  PML  has  been  implicated  in  regulating  cellular 
senescence.  For  example,  the  foci  number  and  size  of  PML 
bodies  increase  during  senescence  (33,  34)  and  inactivation  of 
PML  suppresses  senescence  (35).  Activation  of  the  HIRA/PML 
pathway  is  reflected  by  the  recruitment  of  HIRA  into  PML 
bodies  (36). 

To  determine  whether  Wnt5a  reconstitution  activates  the 
HIRA/PML  senescence  pathway  and  induces  senescence  in 
EOC  cells,  we  first  sought  to  determine  whether  the  HIRA/PML 
pathway  is  conserved  in  human  ovarian  epithelial  cells.  Ectop- 


ically  expressing  activated  oncogenes  (such  as  oncogenic  BAS) 
is  a  standard  approach  for  inducing  senescence  in  a  synchro¬ 
nized  manner  in  primary  human  cells  (1,  2,  19,  20).  Indeed, 
ectopic  expression  of  oncogenic  H-RA!f’llv  induced  senes¬ 
cence  of  primary  HOSE  cells,  as  evident  by  an  increase  in 
SA- (3-gal  activity,  a  universal  marker  of  cellular  senescence 
(Supplementary  Fig.  S3A  and  B).  Notably,  the  HIRA/PML 
pathway  was  activated  during  senescence  of  primary  HOSE 
cells  induced  by  oncogenic  RAS,  as  evident  by  the  relocalization 
of  HIRA  into  PML  bodies  (Supplementary  Fig.  S3C  and  D).  This 
result  shows  that  the  senescence-promoting  HIRA/PML  path¬ 
way  is  conserved  in  human  ovarian  epithelial  cells.  In  addition, 
primary  HOSE  cells  with  HIRA  foci  displayed  a  marked 
decrease  in  BrdU  incorporation,  a  marker  of  cell  proliferation, 
compared  with  HIRA  foci-negative  cells  (Supplementary  Fig. 
S3E  and  F).  This  result  is  consistent  with  the  idea  that  acti¬ 
vation  of  the  HIRA/PML  pathway  is  directly  correlated  with 
senescence-associated  cell  growth  arrest  (37). 

We  next  asked  whether  Wnt5a  expression  is  regulated 
during  natural  senescence  of  primary  HOSE  cells.  Indeed,  we 
observed  an  increase  in  the  levels  of  Wnt5a  mRNA  in  senescent 
primary  HOSE  cells  compared  with  young  cells  (Fig.  4A-C).  In 
addition,  we  found  that  ectopic  Wnt5a  induces  senescence  of 
primary  HOSE  cells  (Fig.  4D-F).  Together,  we  conclude  that 
Wnt5a  plays  a  role  in  regulating  senescence  of  primary  HOSE 
cells. 

As  Wnt5a  antagonizes  canonical  Wnt  signaling  in  human 
EOC  cells  (Fig.  3E  and  F),  we  sought  to  determine  whether 
Wnt5a  restoration  might  activate  the  senescence-promoting 
HIRA/PML  pathway  and  induce  senescence  in  human  EOC 
cells.  Toward  this  goal,  we  examined  the  localization  of  HIRA  in 
OVCAR5  EOC  cells  reconstituted  with  Wnt5a  or  vector  control. 
Notably,  there  was  a  significant  increase  in  the  percentage  of 
cells  with  HIRA  localized  to  PML  bodies  in  Wnt5a  restored 
human  EOC  cells  compared  with  controls  (Fig.  5A  and  B;  P  = 
0.004).  In  addition,  we  also  observed  an  increase  in  the  number 
and  size  of  PML  bodies  in  the  Wnt5a  restored  OVCAR5  EOC 
cells  (Fig.  5A),  which  are  also  established  markers  of  cellular 
senescence  (35,  38).  Similarly,  we  observed  activation  of  the 
HIRA/PML  pathway  by  Wnt5a  restoration  in  PEOI  human 
EOC  cells  (Supplementary  Fig.  S4A  and  B),  suggesting  that  the 
observed  effects  are  not  cell  fine  specific.  Together,  we  con¬ 
clude  that  Wnt5a  reconstitution  activates  the  HIRA/PML 
senescence  pathway. 

The  p53  and  pRB  tumor  suppressor  pathways  play  a  key  role 
in  regulating  senescence  (1).  Thus,  we  sought  to  determine 
whether  activation  of  the  HIRA/PML  pathway  depends  on  the 
p53  and  pRB  pathways.  Interestingly,  pl6INK4a,  the  upstream 
repressor  of  pRB,  is  deleted  in  OVCAR5  human  EOC  cell  line 
(39).  In  addition,  the  levels  of  total  phosphorylated  pRB  were 
not  decreased  by  Wnt5a,  whereas  the  levels  of  cyclin  D1/CKD4- 
mediated  Serine  780  phosphorylation  on  pRB  (pRBpS780)  were 
decreased  by  Wnt5a  (ref.  40;  Fig.  5C  and  D).  Furthermore,  p53  is 
null  in  OVCAR5  cells  (41).  We  conclude  that  activation  of  the 
HIRA/PML  pathway  is  independent  of  the  p53  and  pl6INK4a. 

We  next  sought  to  determine  whether  Wnt5a  restoration 
induces  SA-|3-gal  activity,  a  universal  marker  of  cellular  senes¬ 
cence  (1).  Indeed,  SA-)3-gal  activity  was  notably  induced  by 
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Figure  3.  Wnt5a  restoration  inhibits  the  growth  of  human  EOC  cells  by  antagonizing  canonical  Wnt/p-catenin  signaling.  A,  OVCAR5  cells  were 
transduced  with  a  control  or  Wnt5a-encoding  puromycin-resistant  retrovirus.  The  infected  cells  were  drug-selected  with  3  gg/ml_  puromycin.  Expression  of 
Wnt5a  mRNA  in  drug-selected  cells  was  determined  by  qRT-PCR.  B,  same  as  (A),  but  examined  for  expression  of  Wnt5a  and  fi-actin  in  control 
or  Wnt5a-infected  OVCAR5  and  primary  HOSE  cells  by  immunoblotting.  Relative  levels  of  Wnt5a  expression  was  indicated  on  the  basis  of  the  densitometric 
analysis,  using  NIH  ImageJ  software.  C,  same  as  (A),  but  equal  number  (5,000)  of  drug-selected  control  (open  triangles  and  dotted  Hne)orWnt5a-infected  cells 
(open  circles  and  solid  line)  were  cultured  on  plastic  plates  for  4  days,  and  the  number  of  cells  was  counted  [control  ±  SDorWnt5a±  SD(n  =  3);  Student's?  test 
was  used  for  calculating  Rvalue]  at  day  1  (6,666  ±  1,258  vs.  5,000  ±  1 ,000;  P  =  0.1469),  day  2  (14,583  ±  954  vs.  9,583  ±  3,463;  P  —  0.084),  day  3  (41 ,250  ± 
6,538  vs.  1 4, 750  ±  2,787;  *,  P  =  0.0038),  and  day  4  (83,055  ±  8,978  vs.  35,41 6  ±  2,055;  **,  P  =  0.001).  Mean  of  3  independent  experiments  with  SDand  linear 
regression.  D,  same  as  (C),  but  grown  under  anchorage-independent  condition  in  soft  agar.  The  number  of  colonies  was  counted  2  weeks  after  initial 
inoculation.  Mean  of  3  independent  experiments  with  SD.  E,  same  as  (A),  but  examined  for  the  levels  of  soluble  (3-catenin  and  p-actin  expression  by 


immunoblotting.  NT,  nontreated.  F,  same  as  (A),  but  examined  for  expression  c 
microglobulin  was  used  to  normalize  the  expression  of  indicated  genes.  *,  P  = 

Wnt5a  reconstitution  in  both  OVCAR5  and  PEOl  human  EOC 
cells  compared  with  controls  (Fig.  5E  and  F;  Supplementary 
Fig.  S4C  and  D,  respectively).  On  the  basis  of  these  results,  we 
concluded  that  Wnt5a  restoration  induced  senescence  of 
human  EOC  cells  by  activating  the  HIRA/PML  senescence 
pathway. 

Wnt5a  inhibits  the  growth  of  human  EOC  cells  in  vivo  by 
inducing  cellular  senescence 

We  next  sought  to  determine  whether  Wnt5a  would  mediate 
growth  inhibition  and  induce  senescence  in  vivo  in  an  ortho¬ 
topic  EOC  model  in  immunocompromised  mice.  A  luciferase 
gene  was  retrovirally  transduced  into  control  or  Wnt5a-recon- 
stituted  OVCAR5  cells  to  monitor  the  cell  growth  in  vivo  via 
noninvasive  imaging.  These  cells  were  injected  unilaterally  into 
the  bursa  sac  covering  the  ovary  in  female  immunocompro¬ 
mised  mice  (n  =  6  for  each  of  the  groups;  Supplementary  Fig. 

S5).  Tumor  growth  was  monitored  every  5  days  starting  at  day 
10  postinjection  by  measuring  luciferase  activity,  and  the 
growth  of  the  tumor  was  followed  for  a  total  of  30  days  (Fig. 

6A).  Wnt5a  significantly  suppressed  the  growth  ofxenografted 
OVCAR5  human  EOC  cells  compared  with  controls  (Fig.  6B;  P< 

0.03).  Consistently,  following  general  pathologic  examination 
during  surgical  dissection  at  day  30,  we  observed  that  tumor 


f  indicated  p-catenin  target  genes  by  qRT-PCR.  Expression  of  p-2- 
0.0095;  **,  P  =  0.0012;  and  ***,  P  =  0.0286  compared  with  controls. 

sizes  were  notably  smaller  from  mice  injected  with  Wnt5a- 
reconstituted  OVCAR5  cells  compared  with  controls  (data  not 
shown).  The  expression  of  ectopic  Wnt5a  was  confirmed  by 
IHC  staining  in  sections  from  dissected  tumors  (Fig.  6C). 

We  next  sought  to  determine  whether  cell  proliferation  was 
suppressed  by  Wnt5a  reconstitution  in  dissected  tumors. 
Toward  this  goal,  we  examined  the  expression  of  Ki-67  by 
IHC.  We  observed,  there  was  a  significant  decrease  in  the 
number  of  Ki-67-positive  cells  in  tumors  formed  by  Wnt5a- 
reconstituted  OVCAR5  cells  compared  with  controls  (Fig.  6D 
and  E).  In  addition,  intensity  of  Ki-67  staining  was  also  notably 
weaker  in  Ki-67-positive  Wnt5A-reconstituted  OVCAR5  cells 
than  in  control  Ki-67-positive  cells  (Fig.  6D).  On  the  basis  of 
these  results,  we  conclude  that  Wnt5a  reconstitution  inhibits 
the  proliferation  of  human  EOC  cells  in  vivo  in  an  orthotopic 
xenograft  EOC  model. 

We  next  investigated  whether  the  growth  inhibition 
observed  by  Wnt5a  reconstitution  in  vivo  was  due  to  induction 
of  cellular  senescence.  Toward  this  goal,  we  examined  the 
expression  of  SA-p-gal  activity  in  fresh  sections  of  dissected 
tumors  formed  by  OVCAR5  cells  reconstituted  with  Wnt5a  or 
control  cells.  Indeed,  we  observed  a  significant  increase  in  the 
number  of  cells  positive  for  SA-p-gal  activity  in  OVCAR5 
cells  reconstituted  with  Wnt5a  compared  with  control  tumors 
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Figure  4.  Wnt5a  promotes  senescence  of  primary  HOSE  cells.  A,  young  proliferating  primary  HOSE  cells  were  passaged  to  senescence 
(after  7  population  doublings).  Expression  of  SA-p-gal  activity  was  measured  in  young  and  naturally  senescent  primary  HOSE  cells.  B,  same  as 
(A).  Quantitation  of  SA-p-gal-positive  cells.  **,P<  0.001 .  C,  same  as  (A),  but  mRNA  was  isolated  and  examined  for  Wnt5a  expression  by  qRT-PCR.  Expression 
of  B2M  was  used  as  a  control.  *,  P  =  0.003.  D,  young  primary  HOSE  cells  were  transduced  with  retrovirus  encoding  human  Wnt5a  gene  or  a  control. 
Expression  of  Wnt5a  in  indicated  cells  was  determined  by  qRT-PCR.  Expression  of  B2M  was  used  as  a  control.  E,  same  as  (D),  but  stained  for  expression  of 
SA-p-gal  activity  in  drug-selected  cells.  F,  quantitation  of  (E).  Mean  of  3  independent  experiments  with  SD.  *,  P  <  0.05. 


(Fig.  6F  and  G;  P  =  0.003).  Together,  we  conclude  that  Wnt5a 
reconstitution  inhibits  the  growth  of  human  EOC  cells  in  vivo 
by  inducing  cellular  senescence. 

Discussion 

Driving  cancer  cells  to  undergo  cellular  senescence  has 
recently  been  proposed  to  be  a  novel  mechanism  to  target  for 
developing  cancer  therapeutics  (1,  6).  For  example,  pharma¬ 
cologic  inhibitor  of  PTEN  drives  senescence  and,  consequently, 
inhibits  tumorigenesis  in  vivo  in  xenograft  models  of  PTEN 


heterozygous  prostate  cancer  cells  (42,  43).  Compared  with 
apoptosis,  therapeutics  that  drive  cellular  senescence  are 
proposed  to  have  less  cytotoxic  side  effects  (6),  which  makes 
prosenescence  therapy  attractive.  Herein,  we  describe  that 
restoration  of  Wnt5a  signaling  drives  senescence  of  human 
EOC  cells  both  in  vitro  and  in  vivo  in  an  orthotopic  mouse 
model  of  EOC  (Figs.  5  and  6).  Restoring  gene  expression  by  gene 
therapy  has  had  limited  success.  Therefore,  restoring  Wnt5a 
signaling  via  exogenous  ligand  could  prove  to  be  an  alternative 
approach.  Interestingly,  it  has  been  previously  reported  that  a 
Wnt5a-derived  hexapeptide  is  sufficient  to  restore  Wnt5a 


Figure  5.  Wnt5a  restoration 
triggers  cellular  senescence  in 
human  EOC  cells.  A,  control  and 
Wnt5a-expressing  OVCAR5  EOC 
cells  were  stained  with  antibodies 
to  HIRA  and  PML.  Arrows  point  to 
examples  of  colocalized  HIRA  and 
PML  bodies.  Bar,  10  |im.  B, 
quantitation  of  (A).  A  total  of  200 
cells  from  control  and  Wnt5a- 
expressing  cells  were  examined  for 
HIRA  and  PML  colocalization. 
Mean  of  3  independent 
experiments  with  SD.  C,  same  as 
(A),  but  examined  for  pRB  and 
GAPDH  expression.  D,  same  as  (C), 
but  examined  for  pRBpS780  and 
GAPDH  expression.  E,  same  as  (A), 
but  examined  for  SA-p-gal  activity. 
F,  quantitation  of  (E).  Mean  of  3 
independent  experiments  with  SD. 
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Figure  6.  Wnt5a  restoration  inhibits  tumor  growth  and  promotes  senescence  of  human  EOC  cells  in  vivo.  A,  OVCAR5  cells  were  transduced 
with  luciferase-encoding  hygromycin-resistant  retrovirus  together  with  a  control  orWnt5a-encoding  puromycin-resistant  retrovirus.  Drug-selected  cells  were 
unilaterally  injected  into  the  periovarian  bursa  sac  of  the  female  immunocompromised  mice  (n  =  6  for  each  of  the  groups).  The  radiance  of  luciferase 
bioluminescence,  an  indicator  of  the  rate  for  tumor  growth,  was  measured  every  5  days  from  day  10  until  day  30  by  using  the  IVIS  imaging  system.  Shown  are 
images  taken  at  day  1 0  and  day  30,  respectively.  B,  quantitation  of  tumor  growth  from  injected  OVCAR5  cells  expressing  Wnt5a  or  control  at  indicated  time 
points.  *,  P  =  0.038  compared  with  controls.  C,  following  tumor  dissection,  expression  of  Wnt5a  in  tumors  formed  by  control  or  Wnt5a-expressing 
OVCAR5  EOC  cells  was  determined  by  immunohistochemical  staining  against  Wnt5a  (magnification,  40x).  Bar,  50  |im.  D,  same  as  (C),  but  examined  for 
expression  of  Ki-67,  a  marker  of  cell  proliferation  (magnification,  40x).  Bar,  50  urn.  E,  quantitation  of  (D).  *,  P  =  0.008  compared  with  controls.  F,  expression  of 
SA-fi-gal  activity  was  examined  on  sections  of  fresh-frozen  tumors  formed  by  OVCAR5  cells  expressing  control  or  Wnt5a  (magnification,  40x).  Bar,  1 00  (jm.  G, 
quantitation  of  (F).  *,  P  =  0.003  compared  with  controls.  Arrow  points  to  an  example  of  SA-fi-gal  positive  cells. 


signaling  both  in  vitro  and  in  vivo  in  xenograft  models  of  breast 
cancer  (44).  It  would  be  interesting  to  test  whether  the  Wnt5a- 
derived  hexapeptide  will  be  sufficient  to  reconstitute  Wnt5a 
signaling  and  drive  senescence  of  EOC  cells.  Our  data  suggest 
that  cellular  senescence  is  a  potential  target  for  developing 
EOC  therapeutics.  In  addition,  these  data  imply  that  restora¬ 
tion  of  Wnt5a  signaling  represents  a  potential  novel  strategy  to 
drive  senescence  of  EOC  cells. 

This  study  is  the  first  to  show  a  role  for  Wnt5a  in  regulating 
senescence.  We  showed  that  Wnt5a  activated  the  senescence- 
promoting  HIRA/PML  pathway  in  human  EOC  cells  (Fig.  5A; 
Supplementary  Fig.  S4A).  In  primary  human  cells,  activation  of 
HIRA/PML  pathway  is  sufficient  to  drive  senescence  by  facil¬ 
itating  epigenetic  silencing  of  proliferation-promoting  genes 
(such  as  E2F  target  genes;  ref.  19).  Herein,  we  reported  for  the 
first  time  that  the  key  HIRA/PML  senescence  pathway  can  be 
reactivated  to  drive  senescence  of  human  cancer  cells.  Further 
studies  are  warranted  to  elucidate  the  molecular  basis  by 
which  Wnt5a  restoration  and  activation  of  HIRA/PML  path¬ 
way  drive  cellular  senescence  in  human  EOC  cells. 

Interestingly,  senescence  induced  by  Wnt5a  restoration  in 
human  EOC  cells  was  independent  of  both  the  p53  and 
pl6INKIa  tumor  suppressors,  which  implies  that  EOC  cells  that 
lack  p53  and  pl6INK4a  retain  the  capacity  to  undergo  senes¬ 
cence  via  HIRA/PML  pathway  through  suppressing  the  canon¬ 
ical  Wnt  signaling.  This  is  consistent  with  previous  reports 
showing  that  cancer  cells  that  lack  p53  and  pRB  retain  the 


capacity  to  undergo  senescence  when  treated  with  anticancer 
agents  or  ionizing  radiation  (6).  Notably,  although  the  levels  of 
total  phosphorylated  pRB  were  not  decreased  by  Wnt5a,  we 
observed  a  decrease  in  the  levels  of  pRBpS780  that  is  mediated 
by  cyclin  D1/CDK4  (Fig.  5C  and  D).  Future  studies  will  deter¬ 
mine  whether  the  decrease  in  pRBpS780  levels  plays  a  role  in 
regulating  senescence  of  human  EOC  cells. 

Expression  of  Wnt5a  is  altered  in  many  types  of  cancers  (45). 
For  example,  in  melanoma,  Wnt5a  overexpression  correlates 
with  cancer  progression  and  a  higher  tumor  stage  (16).  How¬ 
ever,  in  colorectal  and  esophageal  squamous  cell  carcinomas, 
Wnt5a  has  been  described  to  be  a  tumor  suppressor  and  was 
frequently  silenced  by  promoter  hypermethylation  (16,  46). 
Consistently,  we  also  observed  Wnt5a  promoter  hypermethy¬ 
lation  in  a  number  of  human  EOC  cell  lines  in  which  Wnt5a  is 
downregulated  (Fig.  2;  Supplementary  Table  SI).  This  result  is 
consistent  with  the  idea  that  Wnt5a  promoter  hypermethyla¬ 
tion  contributes  to  Wnt5a  downregulation  in  human  EOC  cells. 

Wnt5a  function  is  highly  dependent  on  cellular  context  (45). 
For  example,  the  cellular  Wnt  receptor/coreceptor  context 
dictates  the  downstream  signaling  pathways  upon  the  binding 
of  Wnt5a,  which  include  activating  noncanonical  Wnt  signal¬ 
ing  or  antagonizing  canonical  Wnt/(3-catenin  signaling  (47). 
These  reports  illustrate  that  Wnt5a  expression  and  its  resulting 
activity  are  cell  type  and  context  dependent.  The  Wnt  receptor/ 
coreceptor  profile  in  EOC  cells  is  currently  unknown,  and  our 
future  studies  will  elucidate  the  mechanism  by  which  Wnt5a 
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antagonizes  Wnt/p-catenin  signaling  in  human  EOC  cells. 
Regardless,  our  data  show  that  Wnt5a  downregulation  is  an 
independent  predictor  for  overall  survival  in  EOC  patients.  In 
contrast,  2  other  studies  showed  that  higher  Wnt5a  expression 
predicts  poor  survival  in  EOC  patients  (48,  49).  The  basis  for 
this  discrepancy  remains  to  be  elucidated.  An  explanation  may 
be  that  our  study  included  more  cases  than  the  other  2  studies 
(130  EOC  cases  in  our  study  vs.  38  cases  in  the  study  by 
Badiglian  and  colleagues  or  63  cases  in  the  study  by  Peng  and 
colleagues).  It  may  also  be  due  to  the  difference  in  the  com¬ 
position  of  type  I  and  type  II  cases  in  this  study  compared  with 
the  other  2  studies.  The  vast  majority  of  EOC  cases  in  this  study 
are  of  type  II  high-grade  serous  subtypes.  Consistently,  our  data 
showed  that  there  is  a  difference  in  Wnt5a  expression  between 
type  I  and  type  II  EOC  (P  =  0.005;  Table  1).  Furthermore,  it  has 
been  shown  in  microarray  analysis  that  Wnt5a  is  expressed  at 
lower  levels  in  laser  capture  and  microdissected  high-grade 
serous  EOC  compared  with  normal  primary  HOSE  cells  (50). 

In  summary,  the  data  reported  here  show  that  Wnt5a  is  often 
expressed  at  lower  levels  in  human  EOCs  compared  with  either 
normal  human  ovarian  surface  epithelium  or  fallopian  tube 
epithelium.  A  lower  level  of  Wnt5a  expression  correlates  with 
tumor  stage  and  predicts  shorter  overall  survival  in  EOC 
patients.  Reconstitution  of  Wnt5a  signaling  inhibits  the  growth 
of  human  EOC  cells  both  in  vitro  and  in  vivo.  In  addition,  Wnt5a 
reconstitution  suppresses  the  proliferation-promoting  canon¬ 
ical  Wnt/p-catenin  signaling  in  human  EOC  cells.  Significantly, 
Wnt5a  reconstitution  drives  cellular  senescence  in  human  EOC 
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cells  and  this  correlates  with  activation  of  the  senescence- 
promoting  HIRA/PML  pathway.  Together,  our  data  imply  that 
reconstitution  of  Wnt5a  signaling  to  drive  senescence  of 
human  EOC  cells  is  a  potential  novel  strategy  for  developing 
EOC  therapeutics. 
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Cell  Cycle,  Cell  Death,  and  Senescence 

Enhancer  of  Zeste  Homolog  2  Promotes  the  Proliferation 
and  Invasion  of  Epithelial  Ovarian  Cancer  Cells 

Hua  Li1,  Qi  Cai2,  Andrew  K.  Godwin1,  and  Rugang  Zhang1,3 


Abstract 

Enhancer  of  zeste  homolog  2  (EZH2)  is  the  catalytic  subunit  of  the  polycomb  repressive  complex  2  (PRC2)  that 
includes  noncatalytic  subunits  suppressor  of  zeste  12  (SUZ12)  and  embryonic  ectoderm  development  (EED). 
When  present  in  PRC2,  EZH2  catalyzes  trimethylation  on  lysine  27  residue  of  histone  H3  (H3K27Me3), 
resulting  in  epigenetic  silencing  of  gene  expression.  Here,  we  investigated  the  expression  and  function  of  EZH2  in 
epithelial  ovarian  cancer  (EOC).  When  compared  with  primary  human  ovarian  surface  epithelial  (pHOSE)  cells, 
EZH2,  SUZ12,  and  EED  were  expressed  at  higher  levels  in  all  8  human  EOC  cell  lines  tested.  Consistently, 
H3K27Me3  was  also  overexpressed  in  human  EOC  cell  lines  compared  with  pHOSE  cells.  EZH2  was 
significantly  overexpressed  in  primary  human  EOCs  ( n  =  134)  when  compared  with  normal  ovarian  surface 
epithelium  ( n  =  46;  P  <  0.001).  EZH2  expression  positively  correlated  with  expression  of  Ki67  (P  <  0.001;  a 
marker  of  cell  proliferation)  and  tumor  grade  {P  =  0.034)  but  not  tumor  stage  (P  =  0.908)  in  EOC.  There  was  no 
correlation  of  EZH2  expression  with  overall  (P  =  0.3)  or  disease-free  survival  (P  =  0.2)  in  high-grade  serous 
histotype  EOC  patients  ( n  =  98).  Knockdown  of  EZH2  expression  reduced  the  level  of  H3K27Me3  and 
suppressed  the  growth  of  human  EOC  cells  both  in  vitro  and  in  vivo  in  xenograft  models.  EZH2  knockdown 
induced  apoptosis  of  human  EOC  cells.  Finally,  we  showed  that  EZH2  knockdown  suppressed  the  invasion  of 
human  EOC  cells.  Together,  these  data  demonstrate  that  EZH2  is  frequently  overexpressed  in  human  EOC  cells 
and  its  overexpression  promotes  the  proliferation  and  invasion  of  human  EOC  cells,  suggesting  that  EZH2  is  a 
potential  target  for  developing  EOC  therapeutics.  Mol  Cancer  Res;  8(12);  1610—8.  ©2010  AACR. 


Introduction 

Enhancer  of  zeste  homolog  2  (EZH2)  is  the  catalytic 
subunit  of  polycomb  repressive  complex  2  (PRC2;  refs. 
1—4).  In  addition  to  EZH2,  PRC2  also  contains  the  non¬ 
catalytic  subunits  embryonic  ectoderm  development  (EED) 
and  suppressor  of  zeste  12  (SUZ12;  ref.  5).  PRC2  plays  an 
important  role  in  epigenetic  gene  silencing  via  methylation 
of  lysine  27  residue  of  histone  H3  (H3K27)  and  can  add  up 
to  3  methyl  groups  to  the  lysine  side  chain.  EZH2  lacks 
enzyme  activity  on  its  own,  and  has  to  complex  with  EED 
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and  SUZ12  to  attain  robust  histone  methyltransferase 
activity  (5,  6).  The  trimethylated  form  of  H3K27 
(H3K27Me3)  is  thought  to  be  the  main  form  that  confers 
transcriptional  silencing  function  (7-10). 

EZH2  is  overexpressed  in  several  types  of  cancers  (11-15) 
and  is  correlated  with  aggressiveness  and  poor  prognosis  in 
breast  and  prostate  cancers  (11-13).  In  breast  epithelial 
cells,  EZH2  overexpression  causes  anchorage-independent 
growth  and  increases  cell  invasiveness  in  vitro  (11).  In 
prostate  cancer  cells,  inhibition  of  EZH2  blocked  the 
growth  of  prostate  cancer  cells  (13,  15).  In  addition, 
SUZ12  is  also  upregulated  in  certain  types  of  cancer, 
including  colon,  breast,  and  liver  (16-18). 

More  than  85%  of  ovarian  cancers  are  of  epithelial  origin 
(19).  Epithelial  ovarian  cancers  (EOC)  are  classified  into 
distinct  histologic  subtypes  including  serous,  mucinous, 
endometrioid,  and  clear  cell  (19).  The  most  common 
histology  of  EOC  is  serous  (50%-60%  of  all  EOCs), 
approximately,  75%  of  which  is  high-grade  and  25%  is 
low-grade  (20-22).  Less  common  histologies  include  endo¬ 
metrioid  (25%),  clear  cell  (4%),  and  mucinous  (4%;  ref.  20, 
21).  Recently,  an  alternative  classification  has  gained  trac¬ 
tion,  in  which  EOC  is  broadly  divided  into  2  types  (22). 
Type  I  EOC  includes  endometrioid,  mucinous,  low-grade 
serous,  and  clear  cell  carcinomas,  and  type  II  EOC  includes 
high-grade  serous  carcinomas  (22).  EOC  remains  the  most 
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lethal  gynecologic  malignancy  in  the  Western  world  (19). 
Thus,  there  is  an  urgent  need  to  identify  new  targets  for 
developing  novel  therapeutics  for  EOC.  Although  EZH2  is 
overexpressed  in  tumor-associated  endothelial  cells  in  inva¬ 
sive  EOC  (23)  and  regulates  tumor  angiogenesis  in  EOC 
(24),  its  role  in  pathogenesis  of  EOC  remains  poorly 
understood.  Here,  we  examined  the  expression  of  the 
subunits  of  PRC2  and  H3K27Me3  in  human  EOC  cell 
lines.  In  addition,  we  determined  EZH2  expression  in 
primary  human  EOCs  of  different  histologic  subtypes  by 
immunohistochemistry  (IHC).  Further,  we  investigated  the 
effects  of  EZH2  knockdown  by  short  hairpin  RNA 
(shRNA)  on  H3K27Me3  expression,  cell  growth,  and 
invasion  of  human  EOC  cells. 

Material  and  Methods 

Cell  culture 

Primary  human  ovarian  surface  epithelial  (pHOSE)  cells 
were  isolated  and  cultured  as  previously  described  (25).  The 
protocol  was  approved  by  Fox  Chase  Cancer  Center 
(FCCC)  institutional  review  board.  Human  EOC  cell  lines 
A1847,  A2780,  OVCAR3,  OVCAR5,  OVCARIO,  PEOl, 
SKOV3,  and  UPN289  were  kindly  provided  by  Drs. 
Thomas  Hamilton  and  Steve  Williams  at  FCCC  and  were 
maintained  in  1640  medium,  supplemented  with  10% 
FBS,  2  mmol/L  of  L-glutamine,  penicillin  (100  units/ 
mL),  and  streptomycin  (100  pg/mL). 

shRNA,  lentivirus  packaging,  and  infection 

The  sense  sequences  of  2  individual  shRNA  EZH2 
are:  5'-CCAACACAAGTCATCCCATTA-3'  and  5'- 
CGGAAATCTTAAACCAAGAAT-3',  respectively.  Lenti¬ 
virus  packaging  was  performed  using  virapower  system 
(Invitrogen)  according  to  manufacturer's  instruction. 
PEOl  and  SKOV3  at  40%  to  50%  confluence  were 
infected  with  lentivirus  expressing  shRNA  to  the  human 
EZH2  gene  or  vector  control.  The  infected  cells  were  drug- 
selected  with  1  pg/mL  (for  PEOl)  or  3  pg/mL  (for  SKOV3) 
of  puromycin,  respectively. 

Human  ovarian  tissue  microarrays 

Tissue  microarrays,  including  core  samples  from  134 
primary  human  EOCs  and  46  cases  of  normal  ovary  tissues 
were  obtained  from  FCCC  Biosample  Repository  Core 
Facility.  Use  of  these  human  specimens  was  approved  by 
the  Institutional  Review  Board. 

Immunohistochemical  staining  and  scoring 

The  expression  ofEZH2  and  Ki67  proteins  was  detected 
using  avidin-biotin-peroxidase  methods.  Briefly,  tissue 
sections  were  subjected  to  antigen  retrieval  by  steaming 
in  0.01  mol/L  of  sodium  citrate  buffer  (pH  6.0)  for  30 
minutes.  After  quenching  endogenous  peroxidase  activity 
with  3%  hydrogen  peroxide  and  blocking  nonspecific 
protein  binding  with  1%  bovine  serum  albumin,  sections 
were  incubated  overnight  with  primary  monoclonal  anti¬ 
body  (anti-EZH2:  Millipore,  1:100;  anti-Ki67:  DAKO, 


1:100)  at  4°C,  followed  by  biotinylated  goat  anti-mouse 
IgG  (DAKO,  1:400)  for  1  hour,  detecting  the  antibody 
complexes  with  the  labeled  streptavidin-biotin  system 
(DAKO),  and  visualizing  them  with  the  chromogen  3,3'- 
diaminobenzidine.  Sections  were  lightly  counterstained 
with  hematoxylin.  Tissues  in  which  nuclei  were  stained 
for  EZH2  or  Ki67  protein  were  considered  positive.  Two  1- 
mm  cores  were  examined  in  each  specimen  on  the  tissue 
microarray  and  cells  were  counted  in  at  least  5  high-power 
fields,  with  approximately  200  cells  analyzed  per  high- 
power  field. 

FACS,  immunoflurescence  staining,  and  Western  blot 
analysis 

FACS  and  indirect  immunoflurescence  (IF)  staining  were 
performed  as  described  previously  (26-28).  The  following 
antibodies  were  used  for  IF:  rabbit  anti-H3K27Me3  (Cell 
Signaling,  1:1,000),  and  rabbit  anti-H3K9Me3  (Abeam, 
1:500).  The  antibodies  used  for  Western  blotting  were  from 
indicated  suppliers:  mouse  anti-EZH2  (Millipore;  1:2,500), 
rabbit  anti-H3K27Me3  (Cell  signaling,  1:1,000),  rabbit 
anti-H3K9Me3  (Abeam,  1:2,000),  mouse  anti-histone 
H3  (Millipore,  1:10,000),  mouse  anti-GAPDH  (Millipore, 
1:10,000),  rabbit  anti-PARP  p85  fragment  (Promega, 
1:1,000),  rabbit  anti-cleaved  caspase  3  (Cell  Signaling, 
1:1,000),  and  rabbit  anti-cleaved  Lamin  A  (Cell  signaling, 
1:1,000). 

Soft  agar  colony  formation  assay 

A  total  of  1  x  1 04  cells  per  well  were  inoculated  in  a  6-well 
plate  in  1.5  mL  of  RPMI  1640  medium  supplemented  with 
10%  FBS  and  0.35%  agar  on  a  base  layer  of  1.5  mL  of  the 
same  medium  containing  0.6%  agar.  Three  weeks  after 
plating,  the  cells  were  stained  with  1%  crystal  violet  (Sigma) 
in  PBS  to  visualize  the  colonies.  Number  of  colonies  that 
were  larger  than  50  pm  (approximately  100  cells)  in 
diameter  in  each  well  was  counted. 

Matrigel  invasion  assay 

BD  BioCoat  Matrigel  Invasion  Chamber  was  used  to 
measure  cell  invasion  according  to  manufacturer's  instruc¬ 
tion.  Cells  (1  x  105  cells  per  well)  suspended  in  0.5  ml 
RPMI  1640  medium  were  added  to  the  upper  compartment 
of  24-well  matrigel-coated  or  noncoated  8  |im  membrane, 
and  RPMI  1640  medium  supplemented  with  10%  FBS  was 
applied  to  the  lower  compartment.  After  incubating  22 
hours  at  37°C,  5%  C02,  the  cells  were  fixed  with  4% 
formaldehyde  and  stained  with  1%  crystal  violet  in  PBS. 
The  number  of  cells  that  migrated  across  control  membrane 
or  invaded  through  Matrigel-coated  membrane  was  deter¬ 
mined  in  9  fields  across  the  center  and  the  periphery  of  the 
membrane. 

Annexin  V  staining  for  detecting  apoptotic  cells 

Phosphatidylserine  externalization  was  detected  using  an 
Annexin  V  staining  kit  (Millipore)  following  manufacture's 
instruction.  Annexin  V  positive  cells  were  detected  by 
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Guava  system  and  analyzed  with  Guava  Nexin  software 
Module  (Millipore). 

In  vivo  tumorigenicity  assay 

A  total  of  5x  106  cells  in  PBS  (pH  7.3)  per  mouse  were 
injected  subcutaneously  into  the  flank  of  6-week-old  female 
nude  athymic  mice.  The  mice  were  sacrificed  4  weeks  post¬ 
inoculation.  Width  andlengthoftumorsizeweremeasuredand 
the  tumor  volume  (mm3)  was  calculated  using  the  following 
formula:  tumor  volume  (in  mm3)  =  length  X  width2  x  0.52. 

Statistical  analysis 

Quantitative  data  were  expressed  as  mean  ±  SD,  unless 
otherwise  stated,  analysis  of  variance  (AN OVA)  with 
Student's  t  test  was  used  to  identify  significant  differences 
in  multiple  comparisons.  The  X~  test  was  used  to  analyze 
the  relationship  between  categorical  variables.  Overall 
survival  was  defined  as  the  time  elapsed  from  the  date 
of  diagnosis  and  the  date  of  death  from  any  cause  or  the 
date  of  last  follow-up.  Disease-free  survival  was  defined  as 
the  time  elapsed  from  the  date  of  surgery  and  the  date  of 
the  first  recurrence.  Kaplan-Meier  survival  plots  were 
generated  and  comparisons  made  using  the  log-rank 
statistic.  For  all  statistical  analyses,  the  level  of  significance 
was  set  at  0.05. 

Results 

The  catalytic  and  noncatalytic  subunits  of  PRC2 
complex  and  H3K27Me3  are  expressed  at  higher  levels 
in  human  EOC  cell  lines  compared  with  pHOSE  cells 

Expression  of  EZH2,  EED,  and  SUZ12  was  examined  by 
Western  blotting  in  cultures  of  pHOSE  cells  isolated  from 
5  different  individuals  and  8  human  EOC  cell  lines.  When 
compared  with  pHOSE  cells,  EZH2,  EED,  and  SUZ12 
were  expressed  at  higher  levels  in  all  human  EOC  cell  lines 
tested  (Fig.  1A).  Consistently,  the  levels  of  H3K27Me3,  the 
product  of  EZH2  histone  H3  lysine  27  methyltransferase 
activity,  was  also  increased  in  human  EOC  cell  lines  com¬ 
pared  with  pHOSE  cells  (Fig.  1A).  On  the  basis  of  these 
results,  we  conclude  that  the  catalytic  and  noncatalytic 
subunits  of  PRC2  and  H3K27Me3  are  expressed  at  higher 
levels  in  human  EOC  cell  lines  compared  with  pHOSE  cells. 

EZEI2  is  overexpressed  in  primary  human  EOCs  and  its 
expression  positively  correlates  with  expression  of  Ki67, 
a  cell  proliferation  marker 

We  next  sought  to  examine  the  expression  of  EZH2,  the 
catalytic  subunit  of  PRC  complex,  in  134  primary  human 
EOCs  and  46  normal  ovary  tissue  specimens  by  IHC 
staining  (Table  1;  Fig.  IB;  Supplementary  Fig.  SI).  The 
specificity  of  the  EZH2  antibody  used  for  IHC  staining  was 
confirmed  by  the  following  (Supplementary  Fig.  S1A  and 
B).  First,  a  single  band  at  right  molecule  weight  (~95KD) 
was  obtained  in  Western  blotting  of  human  EOC  cell  line 
SKOV3  using  the  EZH2  antibody,  and  this  band  was 
absent  after  expression  of  shRNA  to  the  human  EZH2 
gene  (shEZH2)  that  effectively  knocked  down  EZH2 


Figure  1.  EZH2  is  expressed  at  higher  levels  in  human  EOC  cell  lines 
and  primary  human  EOCs  compared  with  normal  ovarian  surface 
epithelium.  A,  expression  of  EZH2,  EED,  SUZ12,  H3K27Me3,  and  GAPDH 
in  5  individual  batches  of  pHOSE  cell  cultures  and  indicated  human  EOC 
cell  lines  was  determined  by  Western  blotting.  B,  examples  of  EZH2 
immunohistochemical  staining  in  normal  ovary  and  primary  human 
EOC  (shown  is  an  example  of  high-grade  serous  histotype  EOC). 

Bar  =  50  pm. 

mRNA  expression  (Supplementary  Fig.  S1A  and  data  not 
shown).  In  addition,  EZH2  staining  signal  was  lost  when 
primary  anti-EZH2  was  replaced  with  an  isotype-matched 
IgG  control  (Supplementary  Fig.  SIB).  Importantly,  the 
nuclei  of  human  EOC  cells  were  strongly  stained  by  the 
anti-EZH2  antibody  (Fig.  IB;  Supplementary  Fig.  SIC). 
By  contrast,  ovarian  surface  epithelial  cells  were  negative  for 
EZH2  staining  (Fig.  IB). 

We  scored  expression  ofEZH2  as  low  (H-score  <  100)  or 
high  (H-score  >  100)  based  on  the  histochemical  score  (29, 
30),  which  considers  both  the  intensity  of  staining  and  the 
percentage  of  positively  stained  cells.  EZH2  expression  in 
the  surface  epithelium  of  all  46  normal  ovaries  was  scored  as 
low  (in  fact,  negative  EZH2  staining).  EZH2  was  scored  as 
low  in  34%  (46  of  134)  and  high  in  66%  (88  of  134)  of 
primary  EOCs  tested,  respectively.  When  compared 
with  normal  ovarian  surface  epithelium,  EZH2  was 
expressed  at  significantly  higher  levels  in  primary 
human  EOCs  ( P  <  0.001).  Because  EZH2  has  been 
implicated  in  promoting  cell  proliferation  (12),  we  stained 
the  same  set  of  primary  human  EOC  specimens  with  Ki67, 
a  cell  proliferation  marker,  and  compared  the  expression  of 
EZH2  and  Ki67  expression  in  consecutive  sections.  There 
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Table  1.  Correlation  between  EZH2  expression  and  tumor  cell  proliferation  (Ki67)  or  clinicopathologic 
variables 


EXH2  protein  expression 


Patient  characteristics 

Low  (n) 

High  (n) 

Total  (n) 

High  (%) 

P 

Age  (23-85  y,  mean  59.6  y) 

<55 

17 

31 

48 

64.58 

>55 

29 

57 

86 

66.28 

0.843 

Laterality 

Left 

12 

25 

37 

67.57 

Right 

7 

18 

25 

72.00 

Bilaterality 

22 

39 

61 

63.93 

0.764 

Undetermined 

5 

6 

11 

54.55 

Histotype 

Epithelial  ovarian  cancer 

46 

88 

134 

65.67 

Type  1 

11 

25 

36 

69.44 

Low-grade  serous 

2 

0 

2 

Endometrioid 

5 

8 

13 

Mucinous 

1 

4 

5 

Clear  cell 

2 

6 

8 

Others 

1 

7 

8 

Type  II 

35 

63 

98 

64.29 

0.577a 

High-grade  serous 

35 

63 

98 

64.29 

Normal  ovarian  epithelium 

Ki67 

46 

0 

46 

0.00 

<0.001 b 

0%-10% 

27 

2 

29 

6.90 

10%-40% 

10 

17 

27 

62.96 

40%-100% 

9 

69 

78 

88.46 

<0.001 

Tumor  grade 

Well  differentiation 

6 

3 

9 

33.33 

Moderate  differentiation 

10 

12 

22 

54.55 

Poor  differentiation 

29 

70 

99 

70.71 

0.034 

Undetermined 

1 

3 

4 

75.00 

Tumor  stage 

Stage  1/2 

12 

21 

33 

63.64 

Stage  3/4 

33 

55 

88 

62.50 

0.908 

Undetermined 

1 

12 

13 

92.31 

aCompared  with  type  I,  P  =  0.577. 

bCompared  with  epithelial  ovarian  cancer,  P  <  0.001 . 


was  a  significant  correlation  between  EZH2  expression  and 
Ki67  expression  (P  <  0.001;  Table  1).  Together,  we  con¬ 
clude  that  EZH2  is  significantly  overexpressed  in  primary 
human  EOCs  compared  with  normal  ovarian  surface 
epithelium  and  its  expression  correlates  with  a  high  pro¬ 
liferation  index  revealed  by  Ki67  staining. 

EZH2  expression  is  positively  correlated  with  tumor 
grade  but  not  tumor  stage,  or  overall  or  disease-free 
survival 

We  next  sought  to  determine  the  correlation  between 
EZEI2  expression  and  clinical  and  pathologic  features  of 
human  EOCs.  There  was  a  significantly  positive  correlation 
between  EZH2  expression  and  tumor  grade  ( P  =  0.034; 


T able  1 ) .  Elowever,  EZH2  expression  was  not  associated  with 
tumor  stage  (P  =  0.908;  Table  1).  Next,  we  sought  to 
determine  whether  EZH2  expression  correlates  with  prog¬ 
nosis  of  type  II  high-grade  serous  histotype  EOC  patients  for 
which  long-term  follow-up  data  were  available  {n  —  98).  The 
difference  in  overall  ( P=  0.3)  or  disease-free  {P=  0.2)  survival 
between  low  EZEI2  expression  group  ( n  =  35)  and  high 
EZH2  expression  group  ( n  =  63)  was  not  significant  (Fig.  2). 

EZH2  knockdown  inhibits  the  growth  of  human  EOC 
cells  in  vitro  and  in  vivo 

Because  EZH2  expression  positively  correlates  with  Ki67 
expression  (Table  1),  we  sought  to  determine  the  effects  of 
EZH2  knockdown  on  proliferation  ofhuman  EOC  cells.  To 
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A 


Figure  2.  EZH2  expression  does  not  correlate  with  disease-free  or  overall 
survival  in  high-grade  serous  histotype  EOC  patients.  The  univariate 
disease-free  (A)  and  overall  survival  (B)  curves  (Kaplan-Meier  method)  for 
high-grade  serous  histotype  EOC  patients  with  low  or  high  EZH2  protein 
levels  as  assessed  by  immunohistochemistry  analysis. 

knockdown  EZH2  expression  in  SKOV3  cells,  we  devel¬ 
oped  2  individual  lentivirus  encoded  shEZH2.  The  knock¬ 
down  efficacy  of  shEZH2  in  SKOV3  cells  was  confirmed  by 
Western  blotting  (Fig.  3A).  Consistently,  the  level  of 
H3K27Me3  level  was  significantly  reduced  by  shEZH2 
expression  in  SKOV3  cells  as  determined  by  both  Western 
blotting  and  IF  staining  (Fig.  3B  and  C).  As  a  negative 
control,  shEZFI2  expression  has  no  effects  on  the  level  of 
trimethylated  lysine  9  histone  FI3  (H3K9Me3)  that  is 
generated  by  histone  methyltransferase  Suv39H  (ref.  31; 
Fig.  3B  and  C).  Compared  with  controls,  EZH2  knockdown 
significantly  reduced  both  anchorage-dependent  and  -inde¬ 
pendent  growth  in  soft  agar  in  SKOV3  cells  ( P  <  0.001; 
Fig.  3D  and  E) .  The  degree  of  growth  inhibition  by  shEZH2 
correlated  with  the  level  of  EZFI2  knockdown  in  SKOV3 
cells  by  2  individual  shEZH2  (Fig.  3),  suggesting  that  the 
observed  growth  inhibition  by  shEZH2  was  not  because  of 
off-target  effects.  In  addition,  EZFI2  knockdown  in  PEOl 
cells  has  same  effects  on  the  expression  of  H3K27Me3  and 
also  suppressed  both  anchorage-dependent  and  -indepen¬ 
dent  cell  growth  (Supplementary  Fig.  S2),  suggesting  that 
the  observed  growth  inhibition  is  not  cell  line  specific. 


We  next  sought  to  determine  the  effects  of  EZH2  knock¬ 
down  on  the  growth  of  SKOV3  cells  in  vivo  in  immuno¬ 
compromised  nude  mice.  Control  and  shEZH2  expressing 
SKOV3  cells  were  injected  subcutaneously  into  nude  mice 
with  5  x  106  cells  per  mice  and  5  mice  in  each  group.  Four 
weeks  after  injection,  the  sizes  of  xenografted  tumors  were 
compared  between  control  and  shEZH2  expressing  cells 
(Fig.  4A  and  B).  EZH2  knockdown  by  shEZFI2  in  the 
xenograft  tumors  was  confirmed  by  IHC  staining  (Fig.  4C). 
shEZFI2  expression  significantly  inhibited  the  growth  of 
xenografted  SKOV3  cells  (Fig.  4A  and  B). 

EZH2  knockdown  inhibits  the  invasion  of  human  EOC 
cells 

EZH2  has  been  implicated  in  regulating  cell  invasion  in 
several  types  of  cancer  cells  (11,  15,  32,  33).  Thus,  we 
sought  to  determine  the  effects  of  EZH2  knockdown  on 
invasion  of  human  EOC  cells.  Toward  this  goal,  control  and 
shEZH2  expressing  SKOV3  cells  were  tested  for  their 
ability  to  migrate  through  uncoated  control  membrane 
or  invade  through  matrigel-coated  membrane.  Compared 
with  controls,  EZH2  knockdown  significantly  inhibited  the 
invasion  of  SKOV3  cells  as  revealed  by  a  decreased  invasion 
index  that  is  calculated  as  the  ratio  between  the  number  of 
cells  invaded  through  matrigel-coated  membrane  and  the 
number  of  cells  migrated  through  control  membrane 
(Fig.  5).  Inhibition  of  invasion  was  observed  by  2  individual 
shEZH2  in  SKOV3  cells  (Fig.  5).  In  addition,  the  degree  of 
invasion  inhibition  correlated  with  the  degree  of  EZH2 
knockdown  (Figs.  3  and  5),  suggesting  that  this  is  not 
because  of  off-targets  effects.  On  the  basis  of  these  results, 
we  conclude  that  EZH2  knockdown  inhibits  the  invasion  of 
human  EOC  cells. 

EZH2  knockdown  triggers  apoptosis  in  human  EOC 
cells 

We  next  sought  to  determine  the  mechanisms  by  which 
EZH2  knockdown  inhibits  the  growth  of  human  EOC 
cells.  DNA  content  analysis  determined  by  FACS  showed 
that  there  was  no  statistical  difference  in  cell-cycle  distribu¬ 
tion  between  control  and  shEZH2  expressing  cells  (Sup¬ 
plementary  Fig.  S3).  We  next  examined  the  markers  of 
apoptosis  in  control  and  shEZH2  expressing  SKOV3  cells. 
When  compared  with  controls,  markers  of  apoptosis  were 
significantly  induced  by  shEZH2  expression  (Fig.  6).  Those 
apoptotic  markers  include  increased  percentage  of  cells  at 
sub-Gl  phase  as  measured  by  FACS  analysis  (Fig.  6A), 
increased  percentage  of  Annexin  V  positively  stained  cells  as 
measured  by  Guava  Nexin  assay  (Fig.  6B),  upregulation  of 
cleaved  Lamin  A,  PARP  p85,  and  caspase  3  (Fig.  6C;  ref. 
34).  Together,  we  conclude  that  EZH2  knockdown  induces 
apoptosis  of  human  EOC  cells. 

Discussion 

Consistent  with  our  findings,  EZH2  mRNA  expression 
was  upregulated  2-fold  or  more  in  more  than  80%  of  high- 
grade  serous  human  EOC  specimens  in  the  newly  released 
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Figure  3.  EZH2  knockdown  inhibits  the  growth  of  SKOV3  cells  in  vitro.  A,  SKOV3  cells  were  infected  with  2  individual  lentivirus  encoded  shEZH2  or 
control.  Drug-selected  cells  were  examined  for  expression  of  EZH2  and  GAPDH  by  Western  blotting.  B,  same  as  A,  but  examined  for  expression  of  H3K27Me3 
and  H3K9Me3  by  Western  blotting.  C,  same  as  A,  but  examined  for  expression  of  H3K27Me3  and  H3K9Me3  by  immunofluorescence  staining.  DAPI 
counterstaining  was  used  to  visualize  the  cell  nuclei.  D,  same  as  A,  but  equal  number  of  drug-selected  cells  were  seeded  and  counted  at  indicated  time 
points.  Mean  of  3  independent  experiments  with  SD.  *,  P  <  0.001.  E,  same  as  A,  1  x  104  cells  were  seeded  in  soft  agar  and  the  number  of  colonies 
were  counted  after  3  weeks  of  culture.  Mean  of  3  independent  experiments  with  SD.  *,  P  <  0.001 . 


the  Cancer  Genomics  Atlas  (TCGA)  serous  ovarian 
cystoadenocarcinoma  gene  expression  database  (http// 
cancergenome.nih.gov/).  EZH2  gene  is  located  at  chromo¬ 
some  7q36.1.  Gene  amplification  contributes  to  EZH2 
overexpression  in  several  types  of  cancer  (14,  35).  However, 
TCGA  gene  copy-number  analysis  indicates  that  EZH2 
gene  amplification  occurs  only  in  a  very  small  percentage  of 
EOCs  (<10%  specimens  show  >4  copy  of  EZH2  gene; 
http//cancergenome. nih.gov/),  suggesting  that  gene  ampli¬ 
fication  is  not  a  major  mechanism  that  leads  to  EZH2 


upregulation  in  human  EOCs.  EZH2  is  an  E2F  target  gene 
(35).  A  very  recent  study  showed  that  VEGF  stimulates 
EZH2  expression  in  human  EOC  cells  via  E2F  family 
members,  E2F1  and  E2F3  (24).  However,  VEGF  only 
stimulates  the  expression  of  EZH2  mRNA  up  to  3-fold 
(24),  which  is  far  below  the  level  of  increase  in  EZH2 
mRNA  or  protein  in  human  EOC  cells  compared  with 
cultured  pHOSE  cells  (Fig.  1  and  data  not  shown),  suggest¬ 
ing  additional  mechanisms  contribute  to  EZH2  upregula¬ 
tion  in  human  EOC  cells.  In  the  future,  we  will  elucidate 
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Figure  4.  EZH2  knockdown  suppresses  the  growth  of  SKOV3  cells  in  vivo 
in  immunocompromised  mice.  A,  5  x  106  control  or  shEZH2  expressing 
SKOV3  cells  were  injected  subcutaneously  into  immunocompromised 
nude  mice  (n  =  5).  Four  weeks  postinjection,  tumors  were  removed  from 
mice.  B,  quantitation  of  A,  the  size  of  tumors  was  measured.  Mean  of 
tumor  sizes  with  SEM.  C,  xenografted  tumors  formed  by  control  or 
shEZH2  expressing  SKOV3  cells  were  sectioned  and  stained  for  EZH2 
expression.  Bar  =  50  pm. 

additional  mechanisms  that  contribute  to  EZH2  upregula- 
tion  in  human  EOCs. 

EZH2  has  been  demonstrated  as  a  prognostic  marker  for 
breast  and  prostate  cancers  and  positively  correlates  with 
disease-free  survival  and  overall  survival  in  those  patient 
populations  (11-13).  In  addition,  EZH2  overexpression 
correlates  with  more  advanced  disease  stages  of  breast  and 
prostate  cancers  (11,  13).  However,  EZH2  expression  was 
not  a  prognostic  marker  in  other  types  of  cancers  including 
renal  clear  cell  carcinoma  and  hepatocellular  carcinomas 
(33,  36).  We  showed  that  there  was  no  significant  correla¬ 
tion  between  EZH2  expression  and  disease-free  or  overall 
survival  in  high-grade  serous  EOC  patients  (Fig.  2).  Con¬ 
sistent  with  our  findings,  low  expression  of  H3K27Me3  has 
been  demonstrated  to  be  a  poor  prognostic  marker  in  EOC 
(37).  In  contrast,  a  very  recent  study  showed  that  EZH2 
expression  in  either  EOC  cells  or  ovarian  tumor  vasculature 
is  predictive  of  poor  clinical  outcome  (24).  The  basis  for  the 
discrepancy  between  our  study  and  that  of  Lu  et  al.'s  is 
unclear.  A  possible  reason  may  be  that  we  correlated  EZH2 
expression  with  overall  or  disease-free  survival  only  in  high- 
grade  serous  histotype  EOCs  (Fig.  2),  whereas  the  study  by 
Lu  et  al.  includes  additional  histotypes  of  EOCs  (24). 

We  showed  that  EZH2  expression  positively  correlated 
with  Ki67  expression  in  EOCs  (Table  1).  There  are  con¬ 
flicting  results  regarding  the  prognostic  value  of  Ki67  in 
ovarian  carcinoma  (38-42).  A  recent  study  by  Kobel  et  al. 


Figure  5.  EZH2  knockdown  suppresses  the  invasion  of  SKOV3  cells.  A, 
equal  number  of  control  and  shEZH2  expressing  SKOV3  cells  were 
assayed  for  migration  through  uncoated  control  membrane  or  invasion 
through  matrigel-coated  membrane.  The  cells  migrated  through  control 
membrane  or  invaded  through  matrigel-coated  membrane  were  stained 
with  1%  crystal  violet  in  PBS.  B,  quantitation  of  A.  Relative  percentage  of 
shEZH2  expressing  cells  migrated  through  control  membrane  or  invaded 
through  matrigel-coated  membrane  compared  with  controls  was 
indicated.  Mean  of  3  independent  experiments  with  SD.  *,  P  <  0.05. 

C,  invasion  index  of  shEZH2  expressing  SKOV3  cells  compared  with 
controls.  Invasion  index  is  the  ratio  between  cells  invaded  through 
matrigel-coated  membrane  and  cells  migrated  through  control  membrane. 
Mean  of  3  independent  experiments  with  SD.  *,  P  <  0.05. 

(43)  suggests  that  differences  in  Ki67  index  among  different 
subtypes  of  EOCs  confound  the  Ki67  survival  analysis 
because  nearly  all  high-grade  serous  EOCs  have  a  high 
Ki67  index.  In  analysis  of  each  individual  subtypes,  Ki67  is 
no  longer  a  prognostic  marker.  Consistent  with  this, 
although  EZH2  correlates  with  Ki67  expression  (Table  1), 
EZH2  expression  was  not  a  prognostic  indicator  for  either 
overall  or  disease-free  survival  in  the  tested  high-grade 
serous  histotype  EOC  patients  (Fig.  2). 

Interestingly,  when  compared  with  normal  ovarian  surface 
epithelium,  EZH2  expression  is  significantly  upregulated 
(up  to  23-fold)  in  ovarian  epithelial  inclusion  cysts  (44), 
which  are  thought  to  be  the  precursor  lesion  of  a  subset  of 
EOC  (45).  This  suggests  that  EZH2  overexpression  is  an 
early  event  during  EOC  development.  Although  ovarian 
surface  epithelium  is  thought  to  be  the  cell  origin  of  EOC 
(46),  there  are  still  several  histopathology-based  theories  that 
differ  in  their  explanations  about  the  origins  of  EOC  (46- 
48).  Notably,  recent  evidence  suggests  that  a  proportion  of 
high-grade  serous  EOC  may  arise  from  distal  fallopian  tube 
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A 


Figure  6.  EZH2  knockdown  induces  apoptosis  of  SKOV3  cells.  A,  control 
and  shEZH2  expressing  SKOV3  cells  were  examined  for  cell-cycle 
distribution  by  FACS.  The  percentage  of  sub-GI  cells  was  indicated.  B, 
control  and  shEZH2  expressing  cells  were  stained  for  Annexin  V,  a  cell 
surface  marker  of  apoptosis.  Annexin  V-positive  cells  were  measured  by 
Guava  Nexin  assay.  Mean  of  3  independent  experiments  with  SD.  *,  P, 
0.019  and  **,  P  <  0.001  compared  with  controls.  C,  same  as  B,  but 
examined  for  expression  of  cleaved  Lamin  A,  PARP  p85,  and  caspase  3,  all 
markers  of  apoptosis  in  control  and  shEZH2  expressing  cells. 

(48) .  Therefore,  it  will  be  interesting  to  examine  EZH2 
expression  in  normal  fallopian  tube  epithelium. 

Multiple  genes  have  been  implicated  in  EZH2  inhibi¬ 
tion-induced  apoptosis.  For  example,  FBX032  contributes 
to  EZH2  inhibition  induced-apoptosis  in  breast  cancer  cells 
(45),  and  Bim  expression  has  been  demonstrated  to  mediate 
EZH2  inhibition-induced  apoptosis  in  prostate  cancer  cells 

(49) .  Likewise,  E-cadherin,  DAB2IP,  and  SLIT2  have  all 
been  implicated  in  mediating  increased  invasiveness  con¬ 
ferred  by  high  levels  of  EZH2  expression  (32,  50-52). 
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Wnt5a-dependent  induction  of  senescence  suppresses  epithelial  ovarian  cancer. 
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Epithelial  Ovarian  Cancer  (EOC)  remains  the  most  lethal  gynecological  malignancy  in 
US  and  is  the  fifth  leading  cause  of  cancer  deaths  among  American  women.  Thus,  there  is  an 
urgent  need  to  understand  the  etiology  of  EOC  to  develop  novel  therapies  for  this  disease.  Here, 
we  demonstrated  that  a  non-canonical  Wnt  ligand,  Wnt5a,  is  expressed  at  significantly  lower 
levels  in  human  EOC  cell  lines  and  in  primary  human  EOC  compared  with  either  normal  ovarian 
surface  epithelial  cells  or  fallopian  tube  epithelial  cells.  Importantly,  expression  of  Wnt5a  in 
primary  human  EOC  inversely  correlates  with  tumor  stage  but  not  tumor  grade.  Interestingly, 
Wnt5a  expression  is  significantly  lower  in  Type  II  high-grade  serous  EOC  compared  to  Type  I 
EOC  that  includes  low-grade  serous,  mucinous,  clear  cell  and  endometrioid  subtypes  of  EOC.  In 
addition,  we  discovered  that  hypermethylation  of  promoter  CpG  island  contributes  to  Wnt5a 
downregulation  in  human  EOC  cells.  Significantly,  restoration  of  Wnt5a  expression  in  human 
EOC  cells  promoted  senescence  of  EOC  cells  and  resulted  in  a  dramatic  decrease  in  cell 
proliferation  both  in  vitro  and  in  vivo  in  an  orthotopic  model  of  EOC  in  the  ovary  of  SCID  mice. 
Mechanistically,  Wnt5a  inhibited  canonical  Wnt/p-catenin  signaling  and  resulted  in  the 
activation  of  senescence-promoting  histone  repressor  A/PML  pathway.  In  summary,  we  show 
that  Wnt5a  is  often  expressed  at  lower  levels  in  primary  human  EOC  and  Wnt5a  expression 
suppresses  the  growth  of  EOC  cells  by  triggering  senescence  through  antagonizing  canonical 
Wnt  signaling.  These  results  also  suggest  that  loss  of  Wnt5a  expression  is  a  putative  marker  of 
EOC  and  that  non-canonical  Wnt  signaling  is  a  potent  target  for  developing  novel  EOC 
therapeutics. 
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Epithelial  ovarian  cancer  (EOC)  remains  the  most  lethal  gynecological 
malignancy  in  the  US.  Thus,  there  is  an  urgent  need  to  develop  novel  therapeutics  for 
this  disease.  Cellular  senescence  is  a  tumor  suppression  mechanism  that  has  been 
suggested  as  a  novel  mechanism  for  developing  cancer  therapeutics.  Wnt5a  is  a  non- 
canonical  Wnt  ligand  that  plays  a  context-dependent  role  in  human  cancers.  Recently, 
we  discovered  a  lower  level  of  Wnt5a  expression  predicted  shorter  overall  survival  in 
EOC  patients.  Significantly,  Wnt5a  restoration  inhibited  proliferation  of  EOC  cells  both  in 
vitro  and  in  vivo  in  an  orthotopic  mouse  model.  Mechanistically,  Wnt5a  antagonizes 
Wnt/p-catenin  signaling  and  induces  cellular  senescence  by  activating  the  epigenetic 
histone  repressor  A  (HIRA)/promyelocytic  leukemia  (PML)  senescence  pathway.  These 
data  suggest  that  promoting  cells  to  undergo  senescence  by  reconstituting  Wnt5a 
signaling  represents  a  novel  strategy  for  developing  urgently  needed  EOC  therapeutics. 
Here  we  performed  preliminary  experiments  with  a  Wnt5a-mimetic  hexapeptide,  Foxy5. 
In  EOC  cell  lines,  Foxy5  significantly  inhibits  the  proliferation  of  EOC  cells  by 
downregulating  p-catenin  and  its  transcriptional  target,  CCND1  (cyclin  D1).  In  addition, 
Foxy5  inhibits  EOC  cell  migration  in  a  classical  wound-healing  assay.  In  summary,  we 
show  that  loss  of  Wnt5a  predicts  poor  outcome  in  EOC  patients  and  Wnt5a  suppresses 
the  growth  of  EOC  cells  by  triggering  cellular  senescence.  We  suggest  that  strategies  to 
drive  senescence  in  EOC  cells  by  reconstituting  Wnt5a  signaling  via  using  Foxy5  may 
offer  an  effective  new  strategy  for  EOC  therapy. 
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